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An echosounder is used to probe various atmospheric 
parameters. An acoustic wave is transmitted into the 
atmosphere and information deduced from the backscattered 
energy . 

This thesis seeks tp understand the range limitations of 
the echosounder and to explore methods to quantify 
atmospheric turbulence parameters at a given range. The 
propagation of the acoustic energy, including the effects of 
excess attenuation, are modeled to predict the performance 
of an echosounder when various parameters are changed. The 
electronics of an existing echosounder are investigated to 
understand inherent or design limitations. 
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I . INTRODUCTION 



Many atmospheric turbulence-dependent processes take 
place in the troposphere, the lowest ten to twenty 
kilometers above the earth's surface. Acoustic energy 
interacts with irregularities in the atmosphere more 
strongly than electromagnetic waves and is potentially a 
better probe for these irregularities. Acoustic echoes have 
been used to detect: 

1) wind speed and direction profiles, 

2) humidity profiles, 

3) strength and location of temperature inversions, 

4) temperature inhomogeneities, 

5) mechanical turbulence. 

The acoustic sounder, also known by the names 
echosounder, echosonde, sodar, and acoustic radar, transmits 
a pulse of acoustic energy into the atmosphere. The various 
atmospheric parameters can be determined based on the 
intensity and frequency of the scattered energy. The range 
to the scattering center is determined from the elapsed time 
between transmission and return of the scattered energy. 

In this thesis we are attempting to improve our 
understanding of the fundamental range limitations of the 
acoustic sounder. In general, the acoustic sounder suffers 



from multiple scattering. Turbulence at shorter ranges 
alters the phase front of the propagating pulse and reduces 
the magnitude of the return from longer ranges. 

We are also exploring the use of return from shorter 
ranges, along with theory for the degradation of the energy 
along the path, to compensate for short range degradation in 
a boot-strap fashion and quantify the return from longer 
ranges. Acoustic sounders at present indicate the presence 
of inhomogeneities at a specific altitude but the magnitude 
of the inhomogeneities may have errors of a factor of four 
or more. 

The following factors affect the range of the 
echosounder : 

1) Atmospheric: 

a) Pressure, 

b) Temperature, 

2 

c) Temperature Structure Parameter C T , 

2 

d) Velocity Structure Parameter Cv , 

e) Water-vapor Pressure, and 

f) Ambient Noise. 

2) Echosounder: 

a) Antenna Aperture Factor, 

b) Antenna Diameter, 

c) Efficiency, 

d) Frequency, 



7 



e) Power Transmitted, and 



f) Pulse Length. 
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I I . BACKGROUND 



Acoustic energy propagates in the atmosphere as a 
longitudinal wave of pressure variations. The energy is 
scattered, attenuated, and refracted. Energy which is 
scattered constitutes the returned signal. The attenuation 
decreases the energy ensonifying a given volume and 
decreases the returned signal. Energy that is scattered 
through small angles or that is refracted results in further 
degradation or excess attenuation of the signal. 

A. SCATTERING 

The echosounder transmits a pulse of acoustic energy 
that is scattered by temperature and velocity 
inhomogeneities. Information about the inhomogeneities is 
then based on elapsed time between transmission and receipt 
of the return signal, the strength of the returned signal, 
the equations for scattering, and the Doppler shift. 
Variations in the propagation velocity of the wave, which 
are a function of the temperature and velocity variations in 
the air mass, produce the scattering. In addition to 
scattering, the temperature variations cause refraction 
while the velocity turbulence causes both a shift in phase 
and direction of propagation [Ref. l:p. 60]. 



Random temperature or wind structure caused by 
turbulence and uniform gradients of temperature or velocity 
contribute to the scattering. The gradient must change in 
a scale size that is comparable to or smaller than the 
acoustic wavelength in order to contribute to the 
scattering. It appears that the effects of uniform 
gradients are limited to beam bending and possible specular 
reflection for the acoustic frequencies of interest here. 
This refractive structure of the troposphere and 
stratosphere would be better probed with what has come to be 
called infrasound, sound of frequencies below twenty Hertz. 
[Ref. 1 : p . 52] 

The turbulent atmospheric temperature and velocity 

fluctuations, at high wave number, follow the Kolmogorov 
K 3//3 power spectral density law (in one dimension). The 
fluctuations are correlated spatially on the order of one 
centimeter to tens of meters. The nonzero correlation 
lengths and the declining power spectral density imply a 
different scattering process than for random point 
scatterers, even though the turbulence is a stochastic 
process [Ref. l:p 61], 

An expression for the power scattered from a unit volume 
per unit incident flux into a unit solid angle is [Ref. 2:p 
84] 

(J =.Q,55 cos 2 © (£y 2 cos 2 (i) + .13 C± 2 ) ( sin(-|) ) _11/3 

c /L z T z z 



t 



where 



is the wavelength of the transmitted sound, 

T is the average atmospheric temperature, 

0 is the angle of scattering from the direction of 

propagation , 
o 

Cv is the velocity structure parameter, 

Cv 2= <V (x) j- / y(x+r )> 2 , 

and Ct * is the temperature structure parameter, 

Ct 2 = < T( x ) j^(x+r) > 2 . 

This scattering equation indicates that backscatter 
(&=rr) is only a function of the temperature structure. The 
velocity at any point is the vector sum of the phase 
velocity of the sound and the particle velocity of the 
turbulence. This vector sum is always in the forward 
direction because the particle velocity is always less than 
the speed of sound c. Therefore the backscatter is only a 
function of the temperature structure [Ref. l:p 60]. The 
vertical turbulent velocity does cause a Doppler shift of 
the frequency of the backscattered radiation. 

Scattering over a region with correlated scattering 
centers produces constructive and destructive interference. 
The backscattered waves are partially coherent. This 
results in much greater intensities than would be received 
from incoherent scattering such as Rayleigh scattering [Ref. 



l:p 61]. Also the scattering with the interference, over 
the nonzero coherence length, acts like an array of 
scattering centers in a regular crystal lattice of spacing 



L. The Bragg condition 



L= ^ /sin( f> /2) , 



is satisfied. The dominant scattering is for scale sizes of 



B. ATTENUATION 

The atmosphere absorbs some acoustic energy that 
propagates through the atmosphere and reradiates this energy 
at different frequencies. Quantitatively the power lost 
(Pl) over a path length 1 is given by 




where is the attenuation coefficient in Nepers per meter. 
Historically the attenuation coefficient has been considered 
to be the sum of several terms; 

^ 1 = classical viscous losses, 

C X 

ck r = molecular rotation losses, 

o< vib = m °lecular vibrational losses, N 2 and 0 2 • 

Classical and rotational loses are negligible below 
about three kilohertz, the region of interest for 
echosounder operation [Ref. l:p 54 and Ref. 3:p. 18-2], 
consequently vibration of N 2 and 0 2 produce most of the 
attenuation . 




[Ref. 1 :p 61] 
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Rotational and vibrational losses are referred to as 
relaxation processes. The acoustic energy excites internal 
energy modes of the N 2 and O 2 molecules. The rate of 
collisions with water vapor determines the rate and 
efficiency of conversion of the energy into translational 
energy ( heat ) . 

The phase is shifted due to the relaxation processes. 
This is one of the reasons some of the microwave radar pulse 
compression techniques cannot be used in echosounders . A 
number of the pulse compression techniques rely on the phase 
not changing during propagation. 

The dependence of the attenuation on the water vapor 
pressure is believed to be due to a resonance process 

between the lowest vibrational states of the O 2 and N 2 
molecules with the water molecules. Henderson and Hertfeld 
in Reference 4 [p.986] state that the lowest vibrational 

states of 0 2 an d water vapor are only thirty-nine 
wavenumbers (56^K) apart at 1600 cm - ^. For this reason 
oxygen was thought to be primarily responsible for the 
humidity dependent absorption of sound. Henderson and 
Herzfeld [Ref. 4], and many others, assumed nitrogen was an 
inactive dilutent having no effect on the relaxation 
processes [Ref. 4:p. 987]. Unfortunately the theory did not 
agree with the data for low frequencies and high 
humidities . 
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Theory incorporating the relaxation processes involving 



nitrogen and water at low frequencies and relative 
humidities greater than twenty-five percent bring the theory 
into agreement with the data. At high frequencies and low 
relative humidities, theory and data match well with oxygen 
making the main contribution and nitrogen acting as an inert 
dilutent. At low frequencies and high humidities, nitrogen 
seems to make the main contribution and oxygen acts as an 
inert dilutent. For the range of relative humidities in 
which many people live and over a good part of the audible 
frequency range it appears nitrogen is the major contributer 
to the relaxation processes [Ref. 5:p. 165]. This is in 
contrast to the previous assumption that oxygen was 
responsible [Ref. 6:p. 604], 

An expression for attenuation due to the molecular 
absorption is given by [Ref. 7:p. 34], 



^max < 018*f V +/ 2 * ( £/£ m ) 2 ? \ 2 1/2 , 

30 4.8^ fm) ^1 +(f/fm) / „ 

in dB/m, 

where 

f is the frequency of the transmitted pulse, 
f is the Napier frequency, 

f I- 2 \ / r\ \ / ^ ~ \ 0 . 




1- 519 
< 1.8T+492 
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where 

T is temperature in degrees Celcius, 
h is the percent mole ratio of water vapor, 
h= 100 (e/p), 

e= water vapor pressure in mb, 

P= atmospheric pressure in mb. 

To convert from dB/m to Nepers/m note that 
10 log ( I ( x ) /I ( 0 ) )= 10 log ( exp( -X x ) , 

==> 4 . 34 .>< (Nepers/m ) = X(dB/m). 

The Napier frequency is the frequency of maximum 

absorption per wavelength and ^rnax t * ie attenuat i° n at 
the Napier frequency. The Napier frequency is shifted to 
higher frequencies by even small amounts of water vapor. 
Plots of attenuation versus water-vapor pressure in 
millibars (Figure 1) and versus relative humidity (Figure 2) 
are shown for a range of frequencies. 

For a given temperature, the attenuation has a maximum 
and decreases for higher or lower humidities. With higher 
temperatures, the maximum attenuation increases and the 
relative humidity at which that maximum occurs decreases. 
Plots of attenuation versus water-vapor pressure in 
millibars (Figure 3) and versus relative humidity (Figure 4) 
are shown for a range of temperatures. 
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Fig. 1. Attenuation versus Water-vapor Pressure (mb.) 
for Frequencies around i.6 KHz. 




Fig. 2. Attenuation versus 
Relative Humidity (%) 
for Frequencies around 1.6 KHz. 
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ATTENUATION 1/m 



02 



ATMOSPHERIC PRESSURE- 1 0 1 3 mb 
FREQUENCY- 1600 HZ 




Fig. 3. Attenuation versus Water-vapor Pressure (mb.) 
for Temperatures Around 20 Celsius. 




Fig* 4. Attenuation versus Relative Humidity (%) 
for Temperatures around 20 Celsius. 



C. EXCESS ATTENUATION 



The atmosphere scatters and attenuates the acoustic 
energy. In the previous section on scattering, we were only 
concerned with the energy scattered from a given volume. 
Excess attenuation accounts for scattering and energy loss 
for the round trip in the atmosphere to and from a given 
volume of interest. Excess attenuation arises because the 
atmosphere degrades the mutual coherence of the acoustic 
wave. The divergent solid angle of the acoustic wave is 
larger than would occur for a coherent, diffraction limited 
wave . 

The excess attenuation, Ze 2 , accounts for this energy 
lost due to small-angle scattering. This turbulent beam 
broadening reduces the on axis intensity. Clifford and 
Brown in Reference 8 [p. 1972] develop the equation 
Ze = 1/ ( 1 +N ) for N<1 , 

Ze =1 . 5/ (1+N) for N>1, 
where 



N ' 'V/Se > 2 ■ 

D q — antenna diameter, 

is the atmospheric acoustic coherence length 

= 1.46*k 2 * Q f R 0 ds*Cn e 2 (s)*( ^l-s) 5 / 3 +/ s_ j 5/73 ) 

2R 0 ^2R 0 7 

The term varies from a value of one for no excess 



-3/5 
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attenuation to an asymptote of zero, implying the energy 
would be spread over a 2fTsolid angle. 

There is a step in the functional dependence of the 
excess attenuation when N=1 that can be seen in the 
equations above and in the plots. There clearly is no 
physical discontinuity but rather a transition between the 
theoretical dependence between two asymptotic regions. 

Figures 5,6, and 7 are plots of excess attenuation 
versus range for antenna sizes of .5, 1, and 1.5 meters, 
respectively, for several frequencies. As can be seen from 
these plots or the equations above, the relationship of the 
antenna size to the coherence length has a significant 
effect on the excess attenuation. A larger antenna will 
not, by itself, produce greater range. The antenna size 
must be matched to the transmission frequency in terms of 
antenna design guidelines and effects upon the excess 



attenuation . 
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Fig. 5. Range versus Excess Attenuation 
for an Antenna Diameter of .5 m. 
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EXCESS ATTENUATION 



Fig. 6. Range versus Excess Attenuation 
for an Antenna Diameter of 1.0 m. 




Fig. 7. Range versus Excess Attenuation 
for an Antenna Diameter of 1.5 m. 
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Ill EXPERIMENT 



The acoustic sounder is an excellent tool for probing 
the lower troposphere. It suffers from several 
shortcomings : 

1) The range has typically been limited to a few hundred 
meters . 

2) It has proven difficult to quantify the measurements 
accurately for a given range. 

In this thesis, we approached the problem from several 
directions. We looked at the hardware to understand the 
signal to noise limitations and to understand how the 
magnitude of the returned signal could be calibrated. We 
developed a software model which would allow one to estimate 
the power backscattered from the atmosphere at a given 
range, based on profiles of atmospheric characteristics and 
input parameters for the acoustic radar. 

A . HARDWARE 

The acoustic sounder we were working with was the 
Aerovironment Model 300. It consists of an electronic 
module which generates a 1600 Hertz electrical pulse. The 
pulse is converted to acoustic energy by a transducer which 
feeds a 1.25 meter parabolic reflector. Energy 
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backscattered by the atmosphere is then received by the 
reflector and transducer. The electrical signal is then 
filtered and amplified. In addition, a ramp amplifier 
compensates for the 1/r decreasing signal amplitude with 
range to decrease the dynamic range requirements of the 
present data display, a strip chart recorder. 

We replaced the various integrated circuit amplifiers 
and filters with more current designs with lower noise. We 
replaced the pre-amplifier with an OPA 111 and the rest with 
LF 356 BN devices. 

We traced the amplifier and filters of the receiver 
board to understand the undocumented choices the 
manufacturer had made. Figures 8 and 9 show a pre- 
amplifier, high and low pass filters and two stage 
amplification. The filters are bi-quad configured with 
notched outputs. Figure 10 shows the notching. Figure 10 
represents the frequency spectrum output of the receiver 
board as measured with the HP 3561A Frequency Spectrum 
Analyser with random noise from the HP 3561A providing the 
input signal before the filters. Figure 11 represents 
frequency spectrum of the receiver board with the random 

noise across the transducer. A 10^ ohm resistor was used to 
match impedences as shown in Figure 12. 




Fig . 8 Schematic of the Receiver Board 
of the Aerovironment Model 300. 
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. 9. Schematic of the High and Low Pass Filters of 
the Receiver Board of Aerovironment Model 300. 




Fig. 10. Frequency Spectrum of the Receiver Board of 
the Aerovironment Model 300. (Input before filters) 




FREQUENCY (Hz) 



Fig. 11. Frequency Spectrum of the Receiver Board of 
the Aerovironment Model 300. (Input across transducer) 
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Fig. 12. Connection of Noise Source to Transducer 



The manufacturer had a filter-oscillator board after the 
receiver board with a switch for three pass bands; narrow, 
medium, and wide. With the same input as in Figure 10 we 
measured the frequency output for the three settings of the 
filter-oscillator board (Figure 13). As can be seen, there 
is little difference between the three settings. 

From previous measurements, the resistor labeled in 

the schematic (Figure 8) of the receiver board was found to 
load down the input. The resistor was initially 25 kilo- 
ohms. Figures 14, 15, and 16 represent the signal from the 
input transformer with =25K, 100K, and ohms 

respectively. Random noise from the HP 3561A was input 
across the transducer as shown in Figure 12. The Q of the 
input transformer was improved by increasing the resistance. 
The lower curve in each Figure is with no noise signal input 
from the HP 3561A. 

It appeared the Q of the filters could also be increased 
by adding the resistors labeled Rq in Figure 17. The upper 
curves in Figures 18, 19, and 20 are the output of the 
receiver board with each Rq = oc> , 750, and 560 ohms 
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Fig. 13. Frequency Spectrum of the Filter Oscillator Board 
of the Aerovironment Model 300. 

(Input before the filters of the receiver board) 




Fig. 14. Frequency Spectrum of the Input Transformer 
(Input across transducer with 1*25=25^ ) 
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Fig. 15. Frequency Spectrum of the Input Transformer 
(input across transducer with 1 * 25 = 10 °^ ) 




Fig. 16. Frequency Spectrum of the Input Transformer 
(input across transducer with R25= =*?) 
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Fig. 17 . Schematics of the Receiver Board with Rq 
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Fig. 18. Frequency Spectrum of the Receiver Board with Rg= 

(input before filters) 




Fig. 19. Frequency Spectrum of the Receiver Board with 
R(T 750 -A (input before filters) 







Fig. 20. Frequency Spectrum of the Receiver Board with 
r q =560 XL (input before filters) 
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respectively. The input was random noise £rom the HP 3561A 
after the pre-amplifier but before the filters of the 
receiver board. The effect of decreasing the resistance 
was an increase in gain but the pass-band was no longer flat 
on top. However, for the range of vertical wind velocities, 
about two or three meters per second, a pass-band of fifty 
or sixty Hertz is adequate to allow for the Doppler shift. 

As can be seen, the filters could be tuned for a narrow 
pass-band which would have better shape. We would gain 
amplification and still have a pass-band that is wide enough 
to meet our needs. 

The lower curves in Figures 18, 19, and 20 are with 
notching, the upper without. The effect of the notching was 
a sharper cut off near the bandpass frequency (1600 Hz) but 
a loss at higher and lower frequencies. 

B . SOFTWARE 

The signal received by an acoustic sounder from a given 
range tells us that scattering centers exist but we can say 
little about the size and magnitude of the scattering 
centers. We don't know how much energy was incident on the 
scattering volume nor how much the signal was degraded on 
the return path. If, on the other hand, previous signals 
are used to determine the atmospheric characteristics then 



an estimate of the degradation of the incident power and the 
degradation of the return signal can be made. 

The software programs I have written here go in the 
opposite direction. Given certain atmospheric 
characteristics, the power returned is estimated. With the 
modules tested, it would then be a matter of turning it 
around to take actual data and estimate atmospheric 
conditions. In the present form they serve our purposes in 
allowing us to explore the effects of parameter changes. 

All the programs are fundamentally built around the 

echosonde equation, also sometimes referred to as the radar 

equation in meteorology [Ref. 9 and Ref. 7:p. 3]. 

Pr = Er*(Pt*Et)*(exp(-2XR) )* tf (R,f )*(clT )*(AiG)*Ze 2 , 

u 2 R 2 

where 

Pr is the power returned from a range R, 

Pt is the power transmitted at frequency f, 

Er is the efficiency of conversion of acoustic power 
to electrical power by the transducer , 

Et is the efficiency of conversion of electrical power 
to acoustic by the transducer, 
c is the speed of sound in m/sec, 

T is the pulse length, 

A is the area of the antenna, 

R is the range, and 

G is the effective-aperture factor of the antenna. 
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Transducer efficiencies must be measured for each driver 



and typically range from a few percent [Ref. 7:p. 3] to 
twenty five percent [Ref. 10:p. 11-10] . 

The scattering cross section per unit volume ( ’T g ) i- s 

the fraction of incident power backscattered per unit 

distance into a unit solid angle at a given frequency. From 

Reference 7 [p. 4] and Reference 2, 

(f 0 = .0039 * k 1/3 * Ct 2 , 

T 0 

where 

k is the wavenumber = 2 ^/wavelength, 

Tq is the local mean temperature in Kelvin, 

Ct 2 is the temperature structure parameter. 

The power scattered from a scattering volume is 
Pb( I ) = ( Pt*E t-Pb ( 1-1 ) )*exp(-2^R)*(c//2)*Ze 2 * X , 

The power returned to the antenna is 
Pr=Pb*A*G*Er/R 2 , 

where the return path attenuation was already included in 
Ps . The excess attenuation Ze 2 was discussed in the 
background section. 

The dependence with height of the velocity structure 
parameter Cv 2 was needed for the calculation of the excess 
attenuation. Reference 2 and Reference 11 [p.149] give 
Cv 2 = 2 2//3 , 

where £ is the dissipation rate of turbulent kinetic energy. 




Reference 11 [p. 154] and Reference (p 194) give 

<£ = (.33m/s) 3 * (1 + .07 * R 3/5 ) 3/2 , 

k * R 

for a stable surface layer. Figure 21 is a plot of range 
versus this Cv . 

The programs have four temperature structure parameter 
( Cip 2 ) profiles. The operator must choose one. The first is 
based on data as presented in Reference 13 [p . 398] for 
midday clear weather above the Tularosa Basin desert in New 
Mexico. The second choice is for the same data multiplied 
by a factor two to approximate looking up a convective 
plume. The third case for a nocturnal atmosphere assumes a 
dependence proportional to the negative exponential of the 
range as presented in Reference 13 [p. 399] with tower data 

from the same reference used for the first sixty-five 
meters. The fourth case assumes a dependence proportional 
to height to the -4/3 and a surface vertical heat flux of 
.095 [Ref. 7:p. 7]. Also case four allows for the operator 
to input the height of an inversion layer with C T being 
proportional to height to the -4/3 above the inversion 
layer. Figures 22, 23, 24, and 25 are plots of these four 
profiles . 

^ 2 9 

Cv and Ct z were then used to calculate the acoustic 
refractive index factor [Ref. 14:p. 119] 
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Fig. 21. Velocity Structure Parameter (C v 2 ) Profile. 
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Fig. 22. Range versus Temperature Structure Parameter (C T 2) 

for Profile 1 
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Fig. 23. Range versus Temperature Structure Parameter (C^) 

for Profile 2 
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Fig . 
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Range versus Temperature Structure Parameter 
for Profile 3 



(c T 2 ) 




TEMPERATURE STRUCTURE PARAMETER 



Fig. 25. Range versus Temperature Structure Parameter (C T 2 ) 

for Profile 4 
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Cne 2 = Ct 2 + Cv 2 

4T 0 2 c5 

Figure 26 is a plot of Cne^ using the third profile for Ct . 

Figure 27 represents an input flow chart. It summarizes 
the effect of each variable or atmospheric parameter. 

Figures 28 and 29 are plots of the range versus the power 
returned to the echosounder for profiles 3 and 4. 

The segments of the programs were tested against 
existing data to verify proper operation. Based on input 
temperature, atmospheric pressure, water-vapor pressure, and 
the frequency of operation of the acoustic sounder the 
programs calculate the attenuation. If the operator 
desires, the first program will plot the attenuation as a 
function of absolute water-vapor pressure and/or relative 
humidity for frequencies at one-third octaves around trie 
input frequency and then again for temperatures at ten 
degree Celcius intervals around the input temperature. 

These plots were used to check the attenuation against data 
[Ref. 15] and [Ref. 5]. 

Ambient noise levels of acoustic sounders are found to 
be about ten to forty dB above the theoretical Johnson noise 
limit [Ref. 16:p. 19-4]. This noise level determines the 
maximum range. This maximum range is comparable to range 
capabilities of the Aerovironment System 300 when the 
operating parameters of the Aerovironment are used in the 
program . 
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Fig. 27. Input Flow Chart for Computer Model 
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Fig. 28. Range versus Power Returned to the Antenna 
for Or 2 Profile 3, Parameters of the Aerovironment 
1 Model 300, and S.T.P. 




Fig. 29. Range versus Power Returned to the Antenna 
for C T ^ Profile 4, Parameters of the Aerovironment 
Model 300, and S.T.P. 
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Trends, as input parameters where changed, were also 
used to check the program output. As the frequency 
increased, the backscatter of the turbulence in the 
atmosphere decreased slowly and the absorption of the 
atmosphere increased rapidly in agreement with Reference 17 
tp. 1-10]. 

The calculation for excess attenuation was difficult to 
check. At ranges of about 400 to 500 meters the excess 
attenuation of a typical acoustic radar has been found to be 
about .25, in agreement with the program calculation. 

PROGRAM 1 

The first program takes atmospheric parameters and the 

eehosounder parameters as input and outputs nine plots. The 

inputs, by the operator at the keyboard, are: 

Atmospheric parameters; 

Atmospheric pressure in millibars, 

2 

The profile of Ct from four options. 

Temperature in degrees Celsius, 

Water-vapor pressure in millibars, 

Eehosounder data; 

Antenna diameter in meters, 

Frequency of the eehosounder in Hertz, 

Power transmitted by the eehosounder in Watts, and 
Pulse length of the transmitted acoustic eneegy. 
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The program outputs the following plots: 

1. Attenuation (1/m) versus water-vapor pressure (mb) for 
five frequencies at one-third octaves around the input 
frequency. 

2. Attenuation (1/m) versus relative humidity (%) for fiv 
frequencies at one-third octaves around the input 
frequency . 

3. Attenuation (1/m) versus water-vapor pressure (mb) for 
five temperatures at ten degree intervals around the 
input temperature. 

4. Attenuation (1/m) versus relative humidity (%) for fiv 
temperatures at ten degree intervals around the input 
temperature . 

5. Range (m) versus excess attenuation. 

2 

6. Range (m) versus the temperature structure factor Ct • 

2 

7. Range (m) versus the velocity structure factor Cv . 

8. Range (m) versus the acoustic refractive index 
structure factor Cne . 

9. Range (m) versus the power backscattered to the 
echosounder . 

All of the programs prompt the operator for inputs and 
with a series of yes/no questions allows the operator to re 
run with the same inputs or change the inputs. Figure 30 i 
a flow chart of program one and appendix 1 is a copy. 
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Figures 1, 2, 3, 4, 21, 22 , 23 , 24 , 25, 26, 28, and 29 are 
output plots of program one. 

The second program takes the same atmospheric parameters 
used in the first program. With the exception of the 
frequency, it takes the same echosounder parameters also. 
This program outputs a plot of range versus frequency and a 
plot of range versus excess attenuation for various 
frequencies. Figure 31 is a flow chart of the program and 
appendix 2 is a copy. Some output plots of the program are 
included in the conclusions section. 

The third program takes the same atmospheric and 
echosounder parameters as the first. The program outputs a 
plot of range as a function of efficiency of the transducer, 
assuming the same efficiency for transmit and receive. All 
the other programs and plots in this thesis assume 
efficiencies of 25%, which is on the high side of typical 
performance. Appendix 3 is a copy and an output plot is 
presented in the conclusions section. 

The fourth program has the same inputs as the first with 
the exception of antenna area, which is the dependent 
variable for the output plot. It plots the range as a 
function of antenna area for several frequencies and is 
presented and discussed in the conclusions. The model does 
not include different efficiencies based on optimum antenna 
design for different frequencies. The output reflects the 
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Fig. 31. Flowchart of Computer Model Two. 
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effect of different antenna diameters on the excess 
attenuation and the echosonde equation. I am uncertain ho// 
much these drive the considerations for optimum antenna 
design. Appendix 4 is a copy of the program. 

The fifth program has the power to the transducer as the 
dependent variable plotted against range for several 
frequencies. Appendix 5 is a copy and the output plot is 
presented and discussed in the conclusions section. 

The sixth program has the ambient, atmospheric 
background noise as the dependent variable. Appendix 6 
contains a copy of the program and the output plot is 
presented and discussed in the conclusions section. 

These programs, taken together, allow a parametric 
analysis of the effect of different parameter changes. 



IV CONCLUSIONS 



Reference 10 [p. II-l] indicates that as the frequency 
increases : 

1) Background and wind noises decrease except for marked 
peaks due to fans, etc. This relationship between 
noise and frequency is not included in the this model. 

2) The reflectivity of the turbulence in the atmosphere 
decreases slowly and the absorption of the atmosphere 
increases rapidly. The model agrees with this for low 
values of the water-vapor pressure. The effects of 
water-vapor pressure on the attenuation are shown in 
Figures 1 to 4 . 

3) Wildlife sounds tend to increase with increasing 
frequency and dominate the background noise at about 
3000 Hertz. This effect is not included in the model. 

4) The Doppler shift requires the receiver bandwidth to be 
increased. The model doesn't take this into 
account. The bandwidth is used to calculate the 
Johnson noise but Johnson noise is not significant. 

The noise really should be scaled with frequency to 
depict the gain in range accurately for lower 
frequencies . 
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Transducer efficiencies vary with frequency. This is 
not included in the model. An increase in efficiency will 
have a much greater impact on the potential range than 
increases in, say, power. An increase in efficiency will 
help both with transmission and return, increasing the 
transmitted power and the returned electrical signal 
strength. Figures 32 is a plot of the maximum range for 
efficiencies of .05 to .5. This figure is the output of 
program 3. The range increases quickly with improvements in 
efficiency. The efficiency could be improved by using 
better designed horns, such as catenoidal or exponential. 

Optimum antenna diameters vary with frequency. The 
output of program 4, Figure 33, demonstrates this effect. 

Not included are the effects on the antenna effective 
aperture factor G, except to the extent it may be effected 
by the excess attenuation. The discontinuities in the 
curves are due to the discontinuities in the equations for 
the excess attenuation. For a given frequency, the range 
increases and then, as antenna diameters increase beyond an 
optimum, the range decreases. The decrease for larger 
antenna diameters is the result of the excess attenuation. 

One other consideration as to choice of frequency is 
that the resolution is increased with increasing frequency. 

A specific need for the echosounder may drive this 
constraint and is not considered in this model. 
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Fig. 32. Range versus Transducer Efficiency 

Profile 3, Parameters of the Aerovironment 
T Model 300, and S.T.P. 
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Fig. 33. Range versus Antenna Diameter 
for C 2 p r ofile 3, Parameters of the Aerovironment 
Model 300, and S.T.P. 
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Figure 34 is a plot of range as a function of frequency 
for an antenna size of 1.5 meters. The water-vapor pressure 
input was ten millibars. Figure 35 is for the same antenna 
size but with a water-vapor pressure of 2 millibars. The 
step in each curve is not a physical effect but rather due 
to the step of 1.5 in the theoretical equation for the 
excess attenuation. Figure 1 showed the greater value of 
the attenuation at around 2 mb. water-vapor pressure for the 
range of frequencies we are dealing with. These previous 
two plots demonstrate the effect. Figures 34 and 35 are 
some of the output plots of program 2. Both figures 
demonstrate the greater attenuation of the acoustic wave as 
the frequency increases. 

Significant increases in range can be achieved by 
judicious choice of frequency. Considerations are the 
water-vapor pressure (Figures 1 to 4), the excess 
attenuation (Figures 5 to 7), which depends on the antenna 
diameter (Figure 33), the frequency (Figure 33 to 35), and 
the anticipated frequency spectrum of the background noise. 

Figure 35 models the effect of increasing the 
transmitted power to increase the range. This is the output 
plot of program 5. The slope of the curve is not very steep 
and becomes less so for higher frequencies. A considerable 
increase in power is required to double the range. Also the 
model does not include nonlinear transducer effects which 
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Fig. 34. Range versus Echosounder Frequency 
for c«4 Profile 3, Parameters of the Aerovironment 
t Model 300, and S.T.P. 
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Fig. 35. Range versus Echosounder Frequency 
for Profile 3, Parameters of the Aerovironment 

Model 300, and S.T.P. except a Water-vapor Pressure of 2 mb. 
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arise when the power is increased. Reference 1 [p. 57] 
points out that at some point the acoustic wave becomes 
distorted, which implies, from Fourier analysis, there is a 
flow of energy out of the fundamental frequency into the 
higher harmonics. Since the higher frequencies are 
attenuated more quickly we will soon reach a point of 
saturation, which Brown [Ref. 1 :p 57] called nonlinear 
saturation to distinguish it from another saturation effect 
he discusses. 

Figure 37 is the output of program 5. This shows the 
effect of decreasing the ambient background noise level. 

This decrease in noise might be achieved by using digital 
processing and fast Fourier transforms to achieve a narrower 
bandwidth. The bandwidth, and therefore the noise, can be 
much smaller. 

These plots demonstrate which parameter changes might 
best improve the range of an echosounder, which was one of 
our goals in this thesis. The modeled performance changes 
should allow for intelligent decisions of the necessary 
parameters for the expected uses of the echosounder. 

We also sought to explore the ability to use the 
returned signal to quantify atmospheric parameters 
accurately at a given range. As can be seen from Figure 27 
(the input flow chart), for a single acoustic radar return 
you could calculate C T based on an assumed profile of Cv^. 

V. 
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Fig. 36. Range versus Power Transmitted 
for Citi^ Profile 3, Parameters of the Aerovironment 
T Model 300 f and S.T.P. 
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Fig. 37. Range versus Noise Level 
for C,j,2 Profile 3, Parameters of the Aerovironment 
Model 300, and S.T.P. 
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However this need to assume a Cv^ profile could be 
eliminated with the use of two or more echosounders . 
Reference 11 and many others describe techniques. 

Another approach would be to measure the Doppler width. 
The Doppler width is the spread of frequency around the 
Doppler shifted frequency. Epsilon and therefore Cv^ can be 
related to the Doppler width. With the use of fast Fourier 

y 2 

transforms and digital processing Cv and C^, could be 
measured simultaneously with one echosounder. 

Using the returned signal from lesser ranges the energy 
Incident on a given volume could be estimated and the 
degradation of the return signal could be estimated. In 
this way the computer program would allow one to essentially 
boot-strap up to a given range and more accurately depict 
the atmospheric parameters based on the returned signal. 
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APPENDIX A 

COMPUTER PROGRAM 1 
FULLER, ROBERT PROG 1 

1® SEP 85 



*******PURPOSE******* 



I THIS PROGRAM WILL TAKE PARAMETERS OF AN ECHOSOUNDER AND 

! ESTIMATE THE RANGE. THE FOLLOWING INPUTS ARE REQUIRED: 

! ATMOSPHERIC DATA 

! 1 ) ATOMOSPHERIC PRESSURE in millibars 

! 2 )SELECT ONE OF FOUR ATMOSPHERIC PROFILES FOR THE TEMPERATURE 

! STRUCTURE PROFILE. 

! a)FOR ONE PROFILE THE HEIGHT OF THE INVERSION LAYER IS ALSO INPUT 

! 3) TEMPERATURE IN DEGREES CELSIUS 

! 4)UATER VAPOR PRESSURE IN millibars 

! ECHOSOUNDER DATA 

! 5) ANTENNA DIAMETER IN METERS 

! 6 )FREQUENCY OF ECHOSOUNDER IN Hz 

! 7)POWER TRANSMITTED BY ECHOSOUNDER IN WATTS 
! 8 )PULSE LENGTH OF THE TRANSMITTED ACOUSTIC ENERGY 

! THE PROGRAM OUTPUTS THE FOLLOWING PLOTS: 

! 1)ATTENUATION( I/m. ) VERSUS WATER-VAPOR PRESSURE! mb ) FOR FIVE 
! FREQUENCIES AT ONE-THIRD OCTIVES AROUND THE INPUT FREQUENCY. 

! 2 ) ATTENUATION! I /m. > VERSUS RELATIVE HUMIDITY(X) FOR FIVE FREQUENCIES 

! AT ONE-THIRD OCTIVES AROUND THE INPUT FREQUENCY. 

I 3)ATTENUATI0N! 1/m. ) VERSUS WATER-VAPOR PRCSSURE(mb) FOR FIVE 
! TEMPERATURES AT TEN DEGREE INTERVALS AROUND THE INPUT TEMPERATURE. 

! 4)ATTENUATI0N< 1/m. ) VERSUS RELATIVE HUMIDITY! 2) FOR FIVE TEMPERATURES 
! AT TEN DEGREE INTERVALS AROUND THE INPUT TEMPERATURE. 

! 5) RANGE! m. ) VERSUS EXCESS ATTENUATION 

! G) RANGE! m. ) VERSUS TEMPERATURE STRUCTURE FACTOR. 

! 7) RANGE! m. > VERSUS VELOCITY STRUCTURE FACTOR. 

! 8 ) RANGE! m. > VERSUS ACOUSTIC STRUCTURE FACTOR. 

! 9 ) RANGE! m. > VERSUS BACKSCATTERED TO ECHOSOUNDER. 



***** VARIABLES***** 



Again TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 

WISHES TO MAKE ANOUTHER RUN. 

Ano_change TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 

WISHES TO CHANGE THE VALUE OF INPUT DATA BEFORE MAKING 
ANOUTHER RUN. 

Ant_area ECHOSOUNDER ANTENNA AREA CALCULATED FROM INPUT OF 
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Ant_diam 
Atom_jDres 
Att en 

Att_f req 



Att max 



I Att_temp 



Bn 



C 

Cl & CZ 
C3 

CneZ ( * ) 



CteZ( *> 
CveZ ( * ) 



Dx 

Dy 

Epsilon 

Et 

Er 

Es 

Exc_att( * ) 
F 

Fnax 

Freq 

Freq_con 

G 



ANTENNA DIAMETER. 

INPUT OF ECHOSOUNDER ANTENNA DIAMETER IN METERS. 

INPUT OF ATMOSPHERIC PRESSURE IN mb. 

ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 
SUBPROGRAM ATTENUATION. 

VARIABLE USED TO DETERMINE IF OPERATOR 
WANTS TO PLOT ATTENUATION VERSUS UATER_VAPOR 
PRESSURE FOR VARIOUS FREQUENCIES AROUND THE INPUT 
FREQUENCY. IF SO THEN PLOT IS DONE IN SUBPROGRAM 
Att_f req. 

VARIABLE IN SUBPROGRAM ATTENUATION. IT IS THE 
ATTENUATION AT THE FREQUENCY OF THE MAXIMUM 
ATTENUATION FOR THE INPUT CONDITIONS. 

VARIABLE USED TO DETERMINE IF OPERATOR WANTS 
TO PLOT ATTENUATION VERSUS WATER_VAPOR PRESSURE 
FOR VARIOUS TEMPERATURE AROUND THE INPUT TEMPERATURE. 
THE PLOT IS DONE IN THE SUBPROGRAM Att_temp. 

BANDWIDTH OF RECEIVER NEEDED TO RECEIVE SIGNALS 
DOPPLER SHIFFED BY THREE METER PER SECOND 
VERTICAL VELOCITIES. 

VELOCITY OF SOUND CALCULATION AT ZERO CELSIUS. 

VARIABLES USED TO SCALE INVERSION HEIGHT IN 
TEMPERATURE STRUCTURE PROFILE THREE. 

CONSTANT USED IN CALCULATION FOR CteZ PROFILE 4. 

ARRAY OF VALUES OF THE ACOUSTIC REFRACTIVE INDEX 
PARAMETER. CALCULATED BASED ON SELECTION OF 
PROFILE FOR TEMPERATURE STRUCTURE PARAMETER 
AND ASSUMED PROFILE FOR THE VELOCITY STRUCTURE 
PARAMETER. 

ARRAY OF VALUES OF THE TEMPERATURE STRUCTURE PARAMETER. 
THE OPERATOR SELECTS A PROFILE FROM SEVERAL PROVIDED. 

ARRAY OF VALUES OF THE VELOCITY STRUCTURE PARAMETER. 
VALUES BASED ON CALCULATION USING ASSUMED DISSAPATION 
RATE. 

STEP SIZE FOR X AXIS FOR VARIOUS PLOTS. 

STEP SIZE FOR Y AXIS FOR VARIOUS PLOTS. 

DISSAPATION RATE USED IN CALCULATION OF CveZ. 

TRANSMISSION EFFICIENCY OF ECHOSOUNDER. 

EFFICIENCY OF ECHOSOUNDER WHEN RECEIVING ACOUSTIC 
BACKSCATTER. 

SATURATION VAPOR PRESSURE AT GIVEN TEMPERATURE. 

EXCESS " ATTENUATION" AT GIVEN RANGE. 

VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 

FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
ATTENUATION. 

INPUT FREQUENCY OF ECHOSOUNDER. 

USED TO INSURE VALUE OF FREQUENCY PASSED TO SUBPROGRAM 
At t_f req WAS NOT CHANGED. 

ANTENNA EFFECTIVE APERATURE FACTOR. 
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H 

I 


VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 
MAIN INDEX FOR LOOP TO CHECK RETURNED POWER UNTIL 
POWER WAS LESS THAN MINIMUM DETECTABLE SIGNAL. 


Interval 
! Inver 
J 
K 
L 

Mess_up 


MAXIMUM EFFECTIVE-SCATTERING VOLUME THICKNESS. 

HEIGHT OF INVERSION LAYER. 

FIRST ORDER INDEX FOR ASSORTED LOOPS. 

WAVENUMBER 

SECOND ORDER INDEX FOR ASSORTED LOOPS. 

TEST VALUE TO CHECK RESPONSES WHEN OPERATOR IS INPUTTING 
RESPONSES. 


N 


TEST VALUE IN CALCULATION OF COHERENCE LENTH IN CALCULATION 
OF EXCESS ATTENUATION. 


Neu_va 

1 


TEST VALUE IN CHECKING IF OPERATOR WISHES TO CHANGE 
A VARIABLE BEFORE A NEW RUN. 


Noise 


ASSUMED MINIMUM DETECTABLE SIGNAL. 



Pou_back<*) POWER BflCKSCRTTERED FROM GIVEN RANGE . 

Pow_ret < * ) POWER BACKSCATTERED TO ECHOSOUNDER FROM GIVEN RRNGE. 



Pou_trans 

Profile 

Pstar 


INPUT OF POWER SUPPLIED TO TRANSDUCER OF ECHOSOUNDER. 
OPERATOR INPUT OF CteZ PROFILE FROM RVRILRBLE PROFILES. 
VRRIRBLE USED IN SUBPROGRRM RTTENURTION FOR INTERMEDIRTE 
CRLCULRTION. 


Pulse 


TRRNSMITTED PULSE LENGTH IN MILLISECONDS. INPUT BY 
OPERRTOR. 


R 

R1 

Range< * ) 
Remainder 


RRNGE FOR INPUT TO INTEGRRL OF EXCESS RTTENURTION. 

VRRIRBLE USED IN SCRLING INVERSION HEIGHT FOR PROFILE THREE 
RRRRY OF RRNGE VRLUES. 

REMRINDER OF MODULO FUNCTION USED TO DECREES THE 
NUMBER OF PRSSES THROUGH THE INTEGRRTI ON FOR EXCESS 
RTTENURTION. 


[ Rge 
Rho 


VRLUE OF RRNGE IN INTEGRRL FOR EXCESS RTTENURTION. 
CORRELRTI ON LENGTH USED IN CRLCULRTION OF EXCESS 
RTTENURTION. 



Sigma( * ) FRRCTION OF ENERGY BRCK5CRTTERED FROM GIVEN RRNGE. 
SpeecLsound SPEED OF SOUND IN RIR RT INPUT TEMPERATURE. 
Sumpou_back SUM OF THE BACKSCATTERED ENERGY 



T 

Temp 

Temp_con 


INPUT TEMPERATURE IN DEGREES KELVIN. 

INPUT TEMPERATURE IN DEGREES CELSIUS. 

USED TO INSURE VALUE OF TEMPERATURE PASSED TO SUBPROGRAM 
Att_tei*ip WAS NOT CHANGED. 


Temp$ 


VARIABLE STRING USED IN FUNCTION YES. 


Titles 
Tst ar 


STRING PASSED TO SUBPROGRAM PT FOR TITLE OF PLOT. 
INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 
IN SUBPROGRAM ATTENUATION. 


Var 


USED TO SELECT THE VARIABLE THE OPERATOR WISHED TO 
CHANGE BEFORE MAKING ANOTHER RUN. 


Wat_pres 


ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 
OPERATOR. 


w 

X X 


THIRD ORDER INDEX USED IN VARIOUS LOOPS. 

STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 
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! 

! 



! 



! 



! 

! DECLARE VARIABLES 
INTEGER I 

DIM P ou_back ( 1 500 ) f Pou_ret ( 1 500 ) , Range ( I 500 ) , Ct eZ ( 1 500 ) , S igma( 1 500 ) 
DIM CveZ( 1500) ,CneZ< 1 500) , Exc_att< 1500) ,Xvar( 1500) 

DIM T it 1 e$[ 50] ,XlabelSC50] ,Ylabel$[ 15] 
l 

PLOTTER IS 705,"HPGL" 

LINE TYPE 1 
I 

! INPUT ATMOSPHERIC DATA 

INPUT "ENTER TEMPERATURE IN DEGREES CELCI US", Temp 
INPUT "ENTER ATMOSPHERIC PRESSURE IN MILLIBARS" , At om_pres 
INPUT "ENTER WATER VAPOR PRESSURE IN MILLIBARS" , Wat_pres 



RESPONSE TO YES OR NO QUESTION. 

Xlabel $ LABEL ON X AXIS PASSED TO SUBPROGRAM Pt FOR PLOTTING. 

Xmax VALUE OF LARGEST VALUE OF X FOR VARIOUS PLOTS. 

Xmin VALUE OF SMALLEST VALUE OF X FOR VARIOUS PLOTS. 

Xrange VALUE OF RANGE OF X VALUES FOR VARIOUS PLOTS. 

Xvar<*> VALUE OF ARRAY OF X VALUES TO BE PLOTTED IN SUBPROGRAM Pt . 

YlabelS LABEL ON Y AXIS PASSED TO SUBPROGRAM Pt FOR PLOTTING. 

Ymax LARGEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

Ymin SMALLEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

Yrange RANGE OF Y VARIABLES FOR VARIOUS PLOTS. 

Ze INTERMEDIATE VALUE IN CALCULATION OF EXCESS ATTENUATION. 



! 

! INPUT ECHOSOUNDER DATA 

INPUT "ENTER FREQUENCY OF ECHOSOUNDER IN HERTZ" ,Freq 
INPUT "ENTER ANTENNA DIAM IN METERS" , Ant_diam 
INPUT "ENTER ECHOSOUNDER PULSE LENGTH IN MILLI SECONDS=" , Pul se 
INPUT "ENTER POWER TO TRANSMITTER IN WATTS" ,Pou_trans 
Et= . Z5 ! TRANSMIT EFFICIENCY 

Er= . Z5 ! RECEIVER EFFICIENCY 

G s . 40 ! ANTENNA EFFECTIVE APERATURE FACTOR 

! 

! 



! SELECT THE PROFILE FOR THE TEMPERATURE STRUCTURE 
PRINT "YOU HAVE A CHOICE OF TEMPERATURE STRUCTURE PROFILES* 
‘PROFILE TO BE USED." 

*YOUR SELECTIONS ARE" 



PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 



1 A TEMPERARURE STRUCTURE PROFILE BASED ON DATA AS 

PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 WHICH " 
GIVES A HEIGHT TO THE -1.15 PROFILE" 

Z THE TEMPERATURE STRUCTURE PROFILE ABOVE BUT" 

WITH A FACTOR OF TWO TO APPROXIMATE LOOKING" 

UP A CONVECTIVE PLUME" 
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PRINT " 

PRINT M 
PRINT M 
PRINT M 
PRINT M 
PRINT n “ 

PRINT " 

PRINT n 

PRINT " ** 

PRINT M “ 

Mess_up= 1 
WHILE Mess 

INPUT “ENTER THE DESIRED PROFILE (1 OR 2 OR 3 OR 4>" .Profile 
IF Prof il e=4 THEN 

INPUT “ENTER THE HEIGHT OF THE INVERSION IN METERS" , Inver 
END IF 

IF Profile*! THEN 
Mess_Lip=0 

ELSE 

IF Prof i 1 e=Z THEN 
Mess_up=0 

ELSE 

IF Profile-3 THEN 
Mess_up=0 

ELSE 

IF Prof i 1 e=4 THEN 
Mess_up=0 

ELSE 

PRINT Profile, “ WAS NOT ONE OF THE OPTIONS !!!! “ 
Mess_up= 1 
END IF 
END IF 
END IF 
END IF 
END WHILE 
OUTPUT KBD; “K“ ; 

{ 

! 

Again* 1 

UNITIALIZE ARRAYS FOR SUCCESSIVE RUNS 
WHILE Again* 1 
FOR J=1 TO I 

Pou_back ( J ) *0 
Pou_ret ( J ) =0 
Range( J ) =0 
Cte2< J )=0 
S igma( J ) =0 
Cve2( J ) =0 
Cne2 ( J )=0 



3 A TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS 
PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 FOR“ 
NIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT “ 
IS EXP< - . 001 * HEIGHT ABOVE G5 METERS AND HEIGHT" 
TO THE -1.4G BELOW 65 METERS" 

4 A TEMPERATURE STRUCTURE PROFILE BASED ON A“ 

HEIGHT TO THE -4/3 “ 



up= 1 
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Exc_att ( J)=0 
NEXT J 

FOR J-1 TO 15 
PRINT " M 
NEXT J 

PRINT "DO YOU WISH TO SEE A PLOT OF ATTENUATION VERSUS HUMIDITY’ 1 
PRINT M FOR VARIOUS FREQUENCIES AROUND THE FREQUENCY YOU INPUT?" 

LINPUT "IF YES ENTER Y , IF NO ENTER N",X$ 

Att_freq=FNYes< X$) 

IF Att_f req= 1 THEN 

OUTPUT KBD; "K" ; 

CALL Att_f req( Atom_pres , Freq , Temp ) 

END IF 

OUTPUT KBD; "K" ; 

I 

| 

FOR J=1 TO 15 
PRINT " " 

NEXT J 

PRINT "DO YOU WISH TO SEE A PLOT OF ATTENUATION VERSUS HUMIDITY" 

PRINT "FOR VARIOUS TEMPERATURES AROUND THE TEMPERATURE YOU INPUT?" 

LINPUT "IF YES ENTER Y , IF NO ENTER N" ,X$ 

Att_tenp=FNYes< X$) 

IF At t_t emp= 1 THEN 
OUTPUT KBD; "K" ; 

CALL Att_temp< Atomjres , Freq , Temp ) 

END IF 

OUTPUT KBD; "K" ; 
l 

! CONVERT TEMPERATURE TO KELVIN 
T=Temp+Z73 
! 

! CALCULATE SPEED OF SOUND BASED ON INPUT TEMPERATURE( CELCIUS ) 
Speed_sound=Z0. 05* < T ) * . 5 

! CALCULATE THE SPEED OF SOUND AT 0 DEGREES CELCIUS 
C=Z0. 05* 273 A . 5 
! 

! CALCULATE THE ATTENUATION COEFFICIENT BASED ON ATMOSPHERIC DATA 
! EQUATION IS FROM NEFF 1975 

CALL Att enuat ion< Atom_pres , At ten , Freq , Temp , Uat_pres ) 

1=1 ! INDEX FOR LOOP 

Pou_back<0>=0 UNITIALIZE VARIABLE FOR POWER BACKSCATERRED 
K=Z*PI *Freq/ Speed_sound ! WAVENUMBER 
Ant_area=PI*< Ant_diam/Z) A Z ! ANTENNA AREA 
Interval^ Speed_sound*Pul se* 1 . E-3) /Z 

Bn=Z*Freq* < 1 - 1 /< G/Speed_soundH 1 ) ) ! BANDWIDTH FOR 3M/S VERTICAL VELOCITIES 
Noise=1 .38E-Z3*Bn*< T ) +Z . E- 1 4 I MINIMUM DETECTABLE SIGNAL AS 
! JOHNSON NOISE< NEGLIGABLE )+ ESTIMATED BACKGROUND 

I 

I CALCULATE RETURNED SIGNAL POWER UNTIL IT IS LESS THAN BACKGROUND. 
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Range( 0)=0 
Sumpou_back =0 
REPEAT 

RangeC I > “RangeC 1-1 ) + Z 
SELECT Profile 
CASE 1 

! THI S TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
! AS PRESENTED IN UALTERS/KUNDEL 1981 PAGE 398. 

! UAS FOR MIDDAY AND CLEAR UEATHER ABOVE THE 
! TULAROSA BASIN, NEU MEXICO. 

! AN ADDITONAL FACTOR OF 1 . 8 UAS INCLUDED AS PER DIFFERENT 
! AVERAGING TIME. 

CteZC I )=Z. 1Z*Range< I)"<-1 . IB) 

CASE Z 

! THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
! AS PRESENTED IN UALTERS/KUNDEL 1981 PAGE 398. 

! UAS FOR MIDDAY AND CLEAR UEATHER ABOVE THE 
! TULAROSA BASIN, NEU MEXICO. 

! AN ADDITONAL FACTOR OF 1 . 8 UAS INCLUDED AS PER DIFFERENT 
! AVERAGING TIME AND AN ADDITIONAL FACTOR OF TUO UAS 
! INCLUDED TO APPROXIMATE LOOKING UP A CONVECTIVE PLUME 
CteZC I )=Z*2 . 1 Z*Range( I >"<-1 . 1G) 

CASE 3 

! THI S TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS" 

! PRESENTED IN UALTERS/KUNDEL 1981 PAGE 398 FOR” 

I NIGHT CONDITIONS. THE DEPENDENCE UITH HEIGHT M 
IIS EXP( - . 001 ^HEIGHT ABOVE 65 METERS AND HEIGHT" 

! TO THE -1.46 BELOU 65 METERS" 

IF Ranged ><65 THEN 

CteZC I )=75. 5*Range( I > A <-1 .46) 

ELSE 

CteZC I ) =3 . 66E-Z*EXP( - . 00 1 * Range ( I > ) 

END IF 
CASE 4 

[CALCULATE THE TEMPERATURE STRUCTURE FACTOR BASED ON ASSUMPTION 
! THAT IT’S PROPORTIONAL TO HEIGHT TO THE -4/3 POUER 

! THAT THE SURFACE VERTICAL HEAT FLUX IS .095 C.n./sec 

ITHE DEPENDENCE UITH HEIGHT ABOVE THE INVERSION I PULLED FROM 
! MY LEFT EAR 

! EQUATION FROM NEFF, 1975 
SELECT Range( I )/Inver 
CASE <.9 

C3=< ( .0Z4)*(Tr< .667)) 

CteZC I >=C3*< RangeC I > >d-1.33) 

CASE .9 TO 1 

CteZC I )“CteZ( 1-1 > 

R1 =Range( I > 

C1“CteZ< I) 

CASE 1 TO 1.3 

CteZC I )= 10"( C LGTC C 3* RangeC I > >-LGT(C3*R1 > HLGTCCI > ) 
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R2=Range( I ) 

CZ«CteZ< I) 

CASE ELSE 

Cte2< I )=C3*( Range( I ) A ( - 1 . 33 ) - R2 A ( - 1 . 33 ) ) +C2 
END SELECT 
END SELECT 

I CALCULATE VELOCITY STRUCTURE FACTOR. FORMULA FROM GAYNOR 77 
! AND TETARSKI 

Epsilon*( . 2SG6/Range( I ) )*( 1 + . 07* ( Range( I ) ) A ( .6) ) A ( 1.5) 

CveZ( I )®Z*Epsilon A ( .GG7) 

‘CALCULATE ACOUSTIC REFRACTIVE INDEX STRUCTURE FACTOR. 

! FORMULA FROM TETARSKI 19G1 

CneZ ( I )-( CteZ ( I >/< Z . 98E+5 > > + ( CveZ( I )/( C*C ) ) 

‘CALCULATE THE FRACTION OF POWER BACKSCATTERED FROM INTERVAL. 

! EQUATION FROM NEFF. 1975 
S igma( I )*( .0039*< K A ( 1/3) )*CteZ< I ) >/( T) A Z 
! 

‘CALCULATE THE EXCESS ATTENUATION. EQUATION FROM CLIFFORD. 1980 
! THE EXCESS ATTENUATION IS Exc_att(I> 

! THE MODULO STRUCTURE IS TO SKIP SOME OF THE INTEGRALS ONCE THE 
! RESOLUTION IS LESS IMPORTANT. 

IF I <50 THEN 

Remainder 15 © 

ELSE 

IF I < 1 50 THEN 

Remainder®Range( I ) MODULO 10 

ELSE 

Remainder=Range( I ) MODULO 20 
END IF 
END IF 

IF Remainder© THEN 
Rho=0 
L=0 
H=0 

Rge=Range( I ) ! CONSTANT IN INTEGRAL 

R=0 

FOR J=0 TO Z*I 

F=CneZ ( INT( J/Z+1 ) ) 

F=F*( 1 - R/ ( Rge ) A ( 1 . G7 ) + ( R/ ( Rge ) ) A ( 1 . G7 ) ) 

IF J>0 THEN 

IF J<2 * I THEN 

IF INT( J/Z ) = J/Z THEN 
L = L + F 
F=0 

ELSE 

H=H+F 

F=0 

END IF 
END IF 
END IF 
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Rho=Rho+F 

R=R+1 



NEXT J 

Rho®Rho+4*L+Z*H 

Rho®< ( ( Rho* . 33)*K*K*1 .46) A ( - .6) ) 

N=( Ant_diam/Rho ) A Z 
IF N<= 1 THEN 

Ze=1 /( 1+N) 

ELSE 

Ze= 1 . B/ ( 1+N) 

(STEP OF 1 . 5==>SEE CLIFFORD 1380 
END IF 

Exc_att( I )=Ze*Ze 

ELSE 

Exc_att ( I )=Exc_att( 1-1 ) 

END IF 

! 

{CALCULATE THE POUER BACKSC ATTERED 

Pou__back ( I ) = ( Pou_trans*Et-Sumpou_back ) *EXP( - At ten* Range( I ) ) 

Pou_back< I )®Pou_back ( I ) * I nterval *Exc_at t ( I )*Sigma( I ) 

Sumpou_back = Sumpou_back+PoLj_back ( I ) 

{CALCULATE THE POUER BACKSCATTERED TO THE ANTENNA 
! TUO UAY PATH ATTENUATION ACCOUNTED FOR ABOVE . 

Pou_ret( I )=Pou_back( I )*EXP( - Atten*Range( I ) ) * Ant_area*G*Er/Range( I ) A Z 
PRINT "RANGE®" ,Range( I ) 

PRINT "POUER RETURNED®" ,Pou_ret< I ) 

1 = 1 + 1 

UNTIL Pou_ret< 1-1 )>0 AND Pow_ret< 1-1 )<=Noise 

i 
! 

OUTPUT KBD; "K" ; 

l 

PRINT "INPUT CONDITIONS" 

j 

PRINT " " 

PRINT USING "K"; "TEMPERATURE® 

PRINT USING "K" ; "ATMOSPHERIC PRESSURE= 

PRINT USING "K" ; "UATER VAPOR PRESSURE= 

PRINT USING "K" ; "PULSE LENGTH= 

PRINT USING "K" ; "TRANSMITTED FREQUENCY= 

PRINT USING "K"; "ANTENNA DIAMETER® 

PRINT USING "K"; "POUER TRANSMITTED® 

PRINT "TEMPERATURE STRUCTURE PROFILE USED ", Prof lie 
IF Prof il e=4 THEN 

PRINT "INVERSION HEIGHT ", Inver 
END IF 

SELECT Profile 
CASE 1 

PRINT " 1 A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 

PRINT " AS PRESENTED IN UAL TERS/K UNDE L 1381 PAGE 338 " 



" ,Temp,"CELCIUS" 

" , Atom_pres , "mb" 

" , Uat_pr es , " mb" 

" , Pul se , " ms" 

" , Freq , " Hz" 

" , Ant_diam , " m - " 

" ,Pou_trans," UATTS" 
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PRINT " 

CASE Z 

PRINT " 

PRINT M 
PRINT M 
PRINT " 

PRINT ** 

CASE 3 

PRINT " 

PRINT M 
PRINT " 

PRINT " 

PRINT M 
CASE 4 

PRINT " 

PRINT " 

END SELECT 
! 

PRINT USING "K'VTHE MINIMUM DETECTABLE SIGNAL WAS SET AT M , Noise , ,, UATTS M 
| 

PRINT " M 

PRINT "OUTPUT CONDITIONS" 

! 

PRINT USING ”K"; "RANGE®" , Range ( I-Z ) , " m." 
j 

PRINT "HIT CONTINUE TO CONTINUE" 

PAUSE 

OUTPUT KBD; "K" ; 

PRINT "THE FOLLOWING GRAPHS WILL NOW BE PLOTTED WITH RANGE VERSUS" 

PRINT "EXCESS ATTENUATION" 

PRINT "CTEZ (TEMPERATURE STRUCTURE PARAMETER)" 

PRINT "CVEZ (VELOCITY STRUCTURE PARAMETER)" 

PRINT "CNEZ (ACOUSTIC REFRACTIVE INDEX STRUCTURE PARAMETER)" 

PRINT "POWER RETURNED" 

PRINT "WHEN YOU ARE READY FOR THE FIRST GRAPH HIT CONTINUE" 

PAUSE 

OUTPUT KBD; "K" ; 

! 

GRAPHICS ON 
VIEWPORT 15, 1Z0, 10,70 
! 

! 

! PLOT RANGE VERSUS EXCESS ATTENUATION 
Xmin=® 

Xrcax = 1 
Ynin=0 

IF Range( I -Z ) <= 1 000 THEN 
Ymax=3 
ELSE 

Ymax=4 



WHICH GIVES A HEIGHT TO THE -1.16 PROFILE" 

Z A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 

AS PRESENTED IN WALTERS/KUNDEL 1381 PAGE 338 " 
WHICH GIVES A HEIGHT TO THE -1.16 PROFILE" 

BUT WITH A FACTOR OF TWO TO APPROXIMATE" 
LOOKING UP A CONVECTIVE PLUME" 

3 A TEMPERATURE STRUCTURE PROFILE BASED ON DATA" 
AS PRESENTED IN WALTERS/KUNDEL 1381 PAGE 398" 
FOR NIGHT CONDITIONS. THE DEPENDENCE WITH" 
HEIGHT IS EXP( - . 001 * HEIGHT ABOVE G5 METERS" 

AND HEIGHT TO THE -1.46 BELOW G5 METERS" 

4 A TEMPERATURE STRUCTURE PROFILE BASED ON A" 

HEIGHT TO THE -4/3 " 
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END IF 
Dx = . 1 
Dy= 1 

Xrange=Xnax-Xnin 
Y r ange= Ymax - Ym i n 
WINDOW Xmin , Xnax , Ymin , Ynax 
AXES Dx ,Dy ,Xnin, Ynin, 1 ,1 
CLIP OFF 
! 

(LABEL PLOT 
CSIZE 4 , . G 
LDIR 0 
LORG 5 

MOVE . 5 , Ymax* 1 . 1 

LABEL "RANGE VERSUS EXCESS ATTENUATION" 
(LABEL HORIZONTAL AXES 
LDIR 0 
LORG 5 

FOR J=0 TO Xnax STEP Dx 
MOVE J,-- 1 *Dy 
LABEL J 
NEXT J 

MOVE . 5*Xrange , - . 3*Dy 
CSIZE 4,-6 

LABEL "EXCESS ATTENUATION" 

(LABEL VERTICAL AXES 
LORG 8 

FOR J=0 TO Ynax STEP Dy 
CSIZE 4 , . G 
MOVE - . 3*Dx , J 
LABEL "10" 

CSIZE Z 

MOVE 1 *Dx , J+ . 05* Dy 

LABEL USING "K";J 
NEXT J 
LDIR PI/Z 
LORG G 

MOVE -Dx , . 5* Yrange 
CSIZE 4 , , G 

LABEL "RANGE (METERS)" 

CLIP ON 

FOR J=1 TO 1-1 

PLOT Exc_att( J) ,LGT( Range( J) ) 

NEXT J 

PRINT "HIT CONTINUE FOR NEXT PLOT" 



(PLOT RANGE VERSUS TEMPERATURE STRUCTURE PARAMETER 

PAUSE 

GCLEAR 
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I 

IF Range( I-Z ) <* 1 W&Q THEN 
Xmin=-4 
Ymax®3 
ELSE 

Xmin“~5 

Ymax*^ 

END IF 
Xmax*# 

Ymin=# 

Title$="RANGE vs TEMPERATURE STRUCTURE PARAMETER" 

XI abel $=" TEMPERATURE STRUCTURE PARAMETER" 

Ylabel $=" RANGE ( METERS)" 

CALL Log_log( Xmin , Xnax , Ymin, Ymax , CteZ( * ) , I- l , Range ( *),I-1,Title$,Xl abel $ 
tl$) 



! PLOT RANGE VERSUS VELOCITY STRUCTURE PARAMETER 
PAUSE 
GCLEAR 
! 

Xmi n=-Z 
Xmax=# 

Ymin=# 

IF Range( I-Z ><*1000 THEN 
Ymax=3 
ELSE 

Ymax=4 
END IF 

Tit 1 e$=" RANGE VERSUS VELOCITY STRUCTURE PARAMETER" 

XI abel $=" VELOCITY STRUCTURE PARAMETER" 

CALL Log_log< Xnin , Xmax , Ymin , Ymax , CveZ( * ) , I- l f Range ( * ) . I - l f Title$,Xl abel S 
?1 $ > 



! PLOT RANGE VERSUS ACOUSTIC REFRACTIVE INDEX STRUCTURE PARAMETER 

PAUSE 

GCLEAR 

Xmin*-? 

Xmax=- 5 
Ymin*# 

IF Range( I-Z) <=1000 THEN 
Ymax=3 
ELSE 

Ymax =4 
END IF 

Tit 1 e$="RANGE vs ACOUSTIC INDEX STRUCTURE PARAMETER" 

XI abel $*=" ACOUSTI C REFRACTIVE INDEX STRUCTURE PARAMETER- 

CALL Log_l og( Xmin , Xmax , Ymin , Ymax , CneZ( * ) , I- 1 , Range ( * ) f I - 1 , Tit 1 e$ , XI abel $ 
el $ ) 
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! PLOT RANGE VERSUS POUER RETURNED TO ANTENNA 

PAUSE 

GCLEAR 

IF Range* I -Z X-1000 THEN 
Xnin=-1 5 
Xmax=-5 
Ynax 5 ^ 

ELSE 

Xnin*- 17 

Xma 

Ymax 5 ^ 

END IF 
Ymin=0 

Title$="RANGE VERSUS POUER RETURNED TO ANTENNA'* 

XI abel$="POWER RETURNED TO ANTENNA" 

Ylabel$="RANGE (METERS)" 

CALL Log_log( Xnin,Xmax , Ymin, Ymax ,Pou_ret( * ) , 1-1 , Range* * > , I- 1 , Title$,Xlab 
Ylabel$> 

PAUSE 

GCLEAR 



LINPUT “WOULD YOU LIKE TO MAKE ANOTHER RUN (ENTER Y OR N)?",X$ 
Again=FNYes( X$> 

Mess_up= 1 
WHILE Mess_up= 1 
SELECT Again 
CASE 1 

LINPUT "WOULD YOU LIKE TO CHANGE A VARIABLE (ENTER Y OR N>?",)($ 
Neu_va=FNYes( X$) 

Ano_change= 1 
WHI LE Ano_change= 1 
SELECT Neu va 



CASE 1 
PRINT 


» 


VARIABLE CURRENT 


VALUE" 


PRINT 


USING 


"K" ;"1 


TEMPERATURE 


, Temp, "CELSIUS" 


PRINT 


USING 


"K'V'Z 


ATOMOSPHERIC PRESSURE " , Atom_pres , " mb 


PRINT 


USING 


"K" ; " 3 


WATER VAPOR PRESSURE 


" , Wat_pres , "mb" 


PRINT 

PRINT 


• l ii 

USING 


" K " ; "4 


FREQUENCY OF ECHOSOUNDER " ,Freq," Hz" 


PRINT 


USING 




ANTENNA DIAMETER 


" , Ant_diam , " m . " 


PRINT 


USING 


"K" ;"G 


PULSE LENGTH 


" , Pulse , " ms" 


PRINT 


USING 


"K" ;"7 


POWER TRANSMITTED ", 


Pou_trans," WATTS 


PRINT 


USING 


"K" ; "8 


ATMOSPHERIC PROFILE 


" , Prof lie 


PRINT 

PRINT 


n •• 

"ENTER 


THE NUMBER OF THE VARIABLE 


YOU WISH TO" 


PRINT 


"CHANGE" 
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INPUT Var 
SELECT Var 

CASE 1 
INPUT 
CASE Z 
INPUT 
CASE 3 
INPUT 
CASE 4 
INPUT 
CASE 5 

INPUT “ANTENNA DIAMETER IN m. 

CASE G 

’PULSE LENGTH IN ms=", Pulse 



* TEMPERATURE 55 " .Temp 



"ATMOSPHERIC PRESSURE IN mb* 



•UATER VAPOR PRESSURE IN mb= 



'FREQUENCY IN Hz=" f Freq 



, Atom_pr es 
, Uat_pres 



Ant diam 



INPUT 
CASE 7 
INPUT 
CASE 8 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 



•POWER TRANSMITTED IN UATTS=" , Poii_tr ans 



1 



=> Ct A Z PROFILE OF Z A (-I . 16)" 

FROM UALTERS/KUNDEL 1981" 

Z ==> SAME AS ONE BUT WITH FACTOR" 

OF TWO TO APPROXIMATE LOOKING- 
UP A CONVECTIVE PLUME" 

3 ==> Ct A Z PROFILE OF EXP( - . 001 *Z ) " 
FROM UALTERS/KUNDEL 1981." 

4 ==> CT A Z PROFILE OF Z A <-4/3)" 

ENTER NUMBER OF DESIRED PROFILE" , Pr of i 1 e 

IF Prof i 1 e=4 THEN 

INPUT "HEIGHT OF INVERSION IN METERS=" , I nver 
END IF 
CASE ELSE 

PRINT Var, "IS NOT ONE OF THE OPTIONS" 

END SELECT 

LINPUT "MADE ANOTHER CHANGE< ENTER Y OR N)?",X$ 
Neu_va*FNYes( X$) 

Mess_up=Z 
Z 

Mess_up s Z 
Ano_change=Z 
CASE ELSE 

PRINT "YOUR CHOICES WHERE Y OR N MM!" 

Ano__change = 1 

LINPUT "CHANGE A VARIABLE (ENTER Y OR N)?",X$ 
Neu_va=FNYes( X$ ) 

END SELECT 
END WHILE 
CASE Z 

Mess_up=2 
CASE ELSE 

PRINT "YOUR CHOICES WHERE Y OR N 
Mess_up= 1 



CASE 



i I u i i 
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LINPUT “MAKE ANOTHER RUN (ENTER Y OR N ) ?“ , X$ 



Again=FNYes< )($) 
END SELECT 
END WHILE 
END WHILE 

PRINT “THAT’S ALL, FOLKS*' 



end 



! 

! CALCULTE THE ATTENUATION 

sub Attenuat ion( Atom_pres , At ten , Freq , Tenp , Wat_pres ) 

! THIS SUBPROGRAM CALCULATES THE ATTENUATION OF SOUND 
! IN AIR BASED UPON EQUATIONS IN NEFF 1975 



! INPUT ATMOSPHERIC PRESSURE IN MILLIBARS 

! FREQUENCY OF SOUND WAVE IN HERTZ 

! TEMPERATURE IN DEGREES CELCIUS 

! WATER-VAPOR PRESSURE IN MILLIBARS 

! OUTPUT ATTENUATION IN 1 /METERS 

! VARIABLES 

! Aton_pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 

! Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 

! SUBPROGRAM ATTENUATION. 

! Att_max VARIABLE IN SUBPROGRAM ATTENUATION. IT IS THE 
! ATTENUATION AT THE FREQUENCY OF THE MAXIMUM 

! ATTENUATION FOR THE INPUT CONDITIONS. 

! F VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 

! RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 

! Froax FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 

! ATTENUATION. 

! Freq INPUT FREQUENCY OF ECHOSOUNDER. 

I H VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 

! Pstar VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 

! CALCULATION. 

! Temp INPUT TEMPERATURE IN DEGREES CELSIUS. 

! Tstar INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 

! IN SUBPROGRAM ATTENUATION. 

! Wat_pres ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 

! OPERATOR. 

H= 1 00* Wat_pres/Aton_pres 
Tstar=( 1 . 8*Temp+49Z )/51 9 
Pstar= Aton_pres/ 1014 

Fnax=( 10+6600*H+44400*H*H)*Pstar/Tstar A . 8 

Att_nax= . 0078*Fmax*Tstar A ( -Z . 5 )*EXP( 7. 77* ( 1-1/Tstar ) ) 

F=Freq/Fmax 

Atten=( Att_max/304.8>*< ( . 1 8*F ) A Z + ( 2*F*F/( 1 +F*F > > A Z) A .5 

Atten=( Atten4 1 . 74E- 1 0* Freq* Freq ) /4 . 35 

SUBEND 
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def FNYes(XS) 

! THI S FUNCTION INTERPRETS THE OPERATORS RESPONSES TO YES NO QUESTIONS 



! INPUT X$ 

! OUTPUT FNYes 

! VARIABLES 

! Tenp$ VARIABLE STRING USED IN FUNCTION YES. 

I X$ STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 

! RESPONSE TO YES OR NO QUESTION. 

DIM Tenp$[ 1 3 
Temp$[ 1 , 11=TRIM$( X$ ) 

SELECT TempS 
CASE M Y M , M y M 
RETURN I 
CASE M N M i ,, n M 
RETURN 2 
CASE " ” 

RETURN t 
CASE ELSE 

RETURN -2 
END SELECT 

FNEND 

! 

sub Att_f req( Atom_pres ,Freq,Temp) 

| 

! THIS SUBPROGRAM PLOTS THE ATTENUATION OF SOUND VERSUS WATER- 

! VAPOR PRESSURE FOR FIVE DIFFERENT FREQUENCIES AT 1/3 OCTIVE 
! INTERVALS AROUND THE INPUT FREQUENCY. 

! 

I INPUT ATOMOSPHERIC PRESSURE IN MILLIBARS 

! FREQUENCY IN HERTZ 

! TEMPERATURE IN CELCIUS 



OUTPUT PLOT OF ATTENUATION VERSUS WATER-VAPOR PRESSURE 



VARIABLES 

A*tom_pres 

Atten 

Att_f req 



! Es 
! Freq 
! Freq_con 



INPUT OF ATMOSPHERIC PRESSURE IN mb. 

ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 
SUBPROGRAM ATTENUATION. 

VARIABLE USED TO DETERMINE IF OPERATOR 
WANTS TO PLOT ATTENUATION VERSUS WATER_VAPOR 
PRESSURE FOR VARIOUS FREQUENCIES AROUND THE INPUT 
FREQUENCY. IF SO THEN PLOT IS DONE IN SUBPROGRAM 
Att_f req. 

SATURATION VAPOR PRESSURE AT GIVEN TEMPERATURE. 

INPUT FREQUENCY OF ECHOSOUNDER. 

USED TO INSURE VALUE OF FREQUENCY PASSED TO SUBPROGRAM 
Att_f req WAS NOT CHANGED. 
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J 

L 

Temp 

Temp_con 
Wat _ pres 

X 

X$ 

Ynax 



FIRST ORDER INDEX FOR ASSORTED LOOPS. 

SECOND ORDER INDEX FOR ASSORTED LOOPS. 

INPUT TEMPERATURE IN DEGREES CELSIUS. 

USED TO INSURE VALUE OF TEMPERATURE PASSED TO SUBPROGRAM 
ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 
OPERATOR. 

THIRD ORDER INDEX USED IN VARIOUS LOOPS. 

STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 
RESPONSE TO YES OR NO QUESTION. 

LARGEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 



Freq_con=Freq 

Tenp__con=Tenp 

Ynax-0 

Freq=DROUND(Freq*2"<2/3) , 3) 

FOR J=0 TO 4 STEP . I 
Uat_pres=J 

CALL Attenuat ion( Aton_pres , Atten , Freq , Tenp , Uat _pres ) 
IF Atten>Ynax THEN 
Ynax=Atten 
END IF 
NEXT J 

Ynax s= PROUND< Ynax + .005 , -2 ) 

GRAPHICS ON 
VIEWPORT 15,120,10,70 
WINDOW 0,15,0, Ynax 
AXES .5, .005,0,0,2 ,2 
CLIP OFF 
! LABEL PLOT 
CSIZE 4, .5 
LDIR 0 

MOVE 7.5, Ynax* 1 . 1 
LORG 5 

LABEL “ATTENUATION VERSUS WATER- VAPOR PRESSURE 11 
CSIZE 4 , . G 

! LABEL HORIZONTAL AXES 
LDIR 0 
LORG G 

FOR J=0 TO 15 
MOVE J,0 
LABEL J 
NEXT J 

MOVE 7.5, -Ynax* . 1 
CSIZE 4 , . G 

LABEL “WATER-VAPOR PRESSURE mb“ 

(LABEL VERTICAL AXES 
LORG 8 
CSIZE 4 , . G 

FOR J=0 TO Ynax STEP .01 
MOVE . 25 , J 
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LABEL J 
NEXT J 
LDIR PI/2 
LORG B 

MOVE -2 , Ymax* . 5 
CSIZE 4, .6 

LABEL M ATTENUATION 1/m" 

! PRINT INPUT CONDITIONS 
LDIR 0 
LORG 1 

MOVE 5. 5, Ymax 
CSIZE 4,.6 

LABEL USING "K" ; "ATMOSPHERIC PRESSURE 32 " , Atom__pres , " mb" 

MOVE 5.5, Ymax* . 95 
CSIZE 4,. 6 

LABEL USING "K" ; "TEMPERATURE-" , Temp , " CELSIUS" 

J=0 

FOR L= 1 TO 5 

FOR X=0. TO G.5+J STEP .05 
Uat_pres=X 

CALL Attenuat ion( Atom_j3res , At ten , Freq , Temp ,Uat_pres ) 
IF I NT ( L/Z )=L/Z THEN " 

LINE TYPE 2 

ELSE 

LINE TYPE 1 
END IF 

PLOT X , Atten 
NEXT X 
LDIR 0 
LORG 2 
CSIZE 3, .6 
MOVE X , Atten 
LINE TYPE 1 

LABEL USING "K" ; Freq , "Hz" 

Freq^DROUNDf Freq* 2* ( — 1/3) ,3) 

J= J+ 1 . 5 
NEXT L 
! 

PRINT "HIT CONTINUE TO CONTINUE" 

PAUSE 

GCLEAR 

FOR J=0 TO 15 
PRINT " " 

NEXT J 

PRINT " UOULD YOU LIKE TO SEE THE SAME PLOT ONLY- 
PRINT "UITH ATTENUATION VERSUS RELATIVE HUMIDITY?" 

LINPUT "IF YES ENTER Y, ELSE ENTER N" , X$ 

Att_f r eq=FNYes( X$) 

IF Att_f req= 1 THEN 
OUTPUT KBD; "K" ; 
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WINDOW 0, 100,0, Ymax 
CLIP ON 

AXES 5,-005,0,0,2,2 
CLIP OFF 
! LABEL PLOT 
CSIZE 4, .5 
LDIR 0 

MOVE 50, Ymax* 1 . 1 
LORG 5 

LABEL "ATTENUATION VERSUS RELATIVE HUMIDITY" 

CSIZE 4,-5 

! LABEL HORIZONTAL AXES 
LDIR 0 
LORG G 

FOR J=0 TO 100 STEP 10 
MOVE J,0 
LABEL J 
NEXT J 

MOVE 50 , - Ymax* . 1 
CSIZE 4,-6 

LABEL "RELATIVE HUMIDITY %" 

! LABEL VERTICAL AXES 
LORG 8 
CSIZE 4 , . G 

FOR J-0 TO Ymax STEP .01 
MOVE 1 ,J 
LABEL J 
NEXT J 
LDIR PI /Z 
LORG G 

MOVE -15, Ymax* . 5 
CSIZE 4, .6 

LABEL "ATTENUATION 1/m" 

IPRINT INPUT CONDITIONS 
LDIR 0 
LORG 1 

MOVE 40, Ymax 
CSIZE 4, .5 

LABEL USING "K" ; "ATMOSPHERIC PRESSURE^’ , Atom_pres , " mb" 

MOVE 40, Ymax*. 95 
CSIZE 4, .5 

LABEL USING "K" ; " TEMPERATURE=" , Temp , " CELSIUS" 

CLIP ON 

T=Temp+Z73 

Es= 1 0 A ( 9 . 4-2 353/ T ) 

J=0 

Freq=DROUND<Freq_con*Z A (-Z/3> ,3) 

FOR L= 1 TO 5 

FOR X=0. TO 40+J STEP .5 

Wa-t_pres = Es* Atom_pres*X/ ( ( Atom__pres-Es* ( 1 +X/ 100) )* 100) 
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CALL At tenuat ion( Aton_pres , At t en , Freq , Temp , Uat_pres ) 
IF INK L/Z )*L/Z THEN ~ 

LINE TYPE Z 

ELSE 

LINE TYPE 1 
END IF 

PLOT X.Atten 
NEXT X 
LDIR 0 
LORG Z 
CSIZE 3, .6 
MOVE X.Atten 
LINE TYPE 1 

LABEL USING "K" ; Freq , "Hz" 

Freq=DROUND(Freq*ZK 1/3) ,3) 

J=J+1 1 
NEXT L 

PRINT "HIT CONTINUE TO CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT KBD; "K " ; 

ELSE 

GCLEAR 

OUTPUT KBD; " K " ; 

END IF 

Temp=Tenp_con 

Freq=Freq_con 

SUBEND 

! 

! 

sub Att_temp( Aton_pres , Freq , Tenp ) 



THIS SUBPROGRAM PLOTS THE ATTENUATION OF SOUND VERSUS UATER- 
VAPOR PRESSURE FOR FIVE DIFFERENT TEMPERATURES AT 10 DEGREE 
INTERVALS AROUND THE INPUT TEMPERATURE. 



INPUT ATOMOSPHERIC PRESSURE IN MILLIBARS 

FREQUENCY OF ACOUSTIC ENERGY IN HERTZ 
TEMPERATURE IN CELSIUS 

OUTPUT PLOTS OF ATTENUATION VERSUS UATER- VAPOR PRESSURE 



VARIABLES 

At on_pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 

Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 

SUBPROGRAM ATTENUATION. 

At t_t emp VARIABLE USED TO DETERMINE IF OPERATOR UANTS 
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TO PLOT ATTENUATION VERSUS WATER_VAPOR PRESSURE 
FOR VARIOUS TEMPERATURE AROUND THE INPUT TEMPERATURE. 
THE PLOT IS DONE IN THE SUBPROGRAM Att_temp. 

Es SATURATION VAPOR PRESSURE AT GIVEN TEMPERATURE. 

Freq INPUT FREQUENCY OF ECHOSOUNDER. 

I MAIN INDEX FOR LOOP TO CHECK RETURNED POWER UNTIL 

J FIRST ORDER INDEX FOR ASSORTED LOOPS. 

L SECOND ORDER INDEX FOR ASSORTED LOOPS. 

Temp INPUT TEMPERATURE IN DEGREES CELSIUS. 

T emp_con USED TO INSURE VALUE OF TEMPERATURE PASSED TO SUBPROGRAM 

Att_temp WAS NOT CHANGED. 

Wat_pres ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 
OPERATOR. 

X THIRD ORDER INDEX USED IN VARIOUS LOOPS. 

X$ STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 

RESPONSE TO YES OR NO QUESTION. 

Ymax LARGEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 



Ymax=0 

T emp_con=Temp 
Temp=Tenp+20 
FOR J =0 TO 4 STEP .1 
Wat_pres=J 

CALL Attenuat ion( Atom^pree , At ten , Freq , Temp , Wat_pres ) 
IF Atten>Ymax THEN 
Ymax=Atten 
END IF 
NEXT J 

Ymax=PROUND( Ymax*+ . 00E , -2 ) 

GRAPHICS ON 
VIEWPORT 15,120, 10,70 
WINDOW 0 , 15,0, Ymax 
AXES .5,-005,0,0,2,2 
CLIP OFF 
! LABEL PLOT 
CSIZE 4 , . B 
LDIR 0 

MOVE 7.5, Ymax* 1 . 1 
LORG 5 

LABEL "ATTENUATION VERSUS WATER- VAPOR PRESSURE" 

CSIZE 4 , . B 

! LABEL HORIZONTAL AXES 
LDIR 0 
LORG B 

FOR J=0 TO 15 
MOVE J,0 
LABEL J 
NEXT J 

MOVE 7.5 , - Ymax* . 1 



8o 



CSIZE 4,. 6 

LABEL "WATER- VAPOR PRESSURE nb" 

! LABEL VERTICAL AXES 
LORG 8 
CSIZE 4 , . G 

FOR J=0 TO Ynax STEP .01 
MOVE .25, J 
LABEL J 
NEXT J 
LDIR PI/2 
LORG B 

MOVE -2, Ynax*. 5 
CSIZE 4 , . 6 

LABEL "ATTENUATION 1/n" 

! PRINT INPUT CONDITIONS 

LDIR 0 

LORG 1 

MOVE G , Ynax 

CSIZE 4 , . B 

LABEL USING "K" ; "ATMOSPHERIC PRESSURE =" , Aton_pres , " nb" 

MOVE G, Ynax*. 95 
CSIZE 4 , . G 

LABEL USING "K" ; "FREQUENCY*" f Freq , " HZ" 

CLIP ON 
J-0 

FOR L= 1 TO 5 

FOR X=0. TO G+J STEP .05 
Uat_pres=X 

IF INT<L/Z)=L/Z THEN 
LINE TYPE 3 

ELSE 

LINE TYPE 1 
END IF 

CALL Attenuat ion( Aton_pres , At ten , Freq , Tenp , Uat_pres ) 
PLOT X , At ten 
NEXT X 
LDIR 0 
LORG 2 
CSIZE 3, .5 
MOVE X.Atten 
LINE TYPE 1 

LABEL USING "K" ; Tenp , " C" 

Tenp=Tenp- 1 0 
J=J+ 1 . G 
NEXT L 

PRINT "HIT CONTINUE TO CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT KBD; "K" ; 

FOR J-0 TO 15 
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PRINT " M 
NEXT J 
! 

PRINT "WOULD YOU ARE LIKE TO SEE THE SAME PLOT ONLY" 

PRINT "WITH ATTENUATION VERSUS RELATIVE HUMIDITY" 

LINPUT "IF SO ENTER Y, IF NOT ENTER N",X$ 

Att_tenp=FNYes( X$ ) 

IF Att_temp= 1 THEN 
OUTPUT KBD; "K" ; 

Temp=Tenp_con 
WINDOW 0J00,0,Ynax 
AXES 10, .005,0,0,Z,Z 
CLIP OFF 
! LABEL PLOT 
CSIZE 4, .6 
LDIR 0 

MOVE 50, Ynax* 1 . 1 
LORG 5 

LABEL "ATTENUATION VERSUS RELATIVE HUMIDITY" 

! LABEL HORIZONTAL AXIS 
LDIR 0 
LORG G 

FOR J=0 TO 100 STEP 10 
MOVE J,0 
LABEL J 
NEXT J 

MOVE 50, -Ynax*. 1 
CSIZE 4 , . G 

LABEL "RELATIVE HUMIDITY" 

! LABEL VERTICAL AXES 
LORG 8 
CSIZE 4, .6 

FOR J=0 TO Ynax STEP .01 
MOVE Z,J 
LABEL J 
NEXT J 
LDIR PI/Z 
LORG G 

MOVE -10, Ynax*. 5 
CSIZE 4, .6 

LABEL "ATTENUATION 1/n" 

IPRINT INPUT CONDITIONS 

LDIR 0 

LORG G 

MOVE 50, Ynax 

CSIZE 4 , . G 

LABEL USING "K" ; "ATMOSPHERIC PRESSURE^ 1 , Aton_pres , " nb" 
MOVE 50, Ynax*. 95 
CSIZE 4 , . G 

LABEL USING "K" ; "FREQUENCY^" , Freq , " Hz" 
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CLIP ON 
J=0 

Temp“Temp-Z0 
FOR L=1 TO 5 
T=Temp+Z73 
Es=10 A ( 9 . 4-2 353/ T ) 

FOR X=0 TO 90- J STEP .5 

L)at_pres aB Es* Aton_pres*X/( ( At om_pres-Es* ( 1 +X/ 1 00) )* 100) 
IF INT<L/2>-L/2 THEN 
LINE TYPE 3 

ELSE 

LINE TYPE 1 
END IF 

CALL Attenuation* Atom_pres, Atten,Freq,Temp , Uat__pres) 
PLOT X , Atten 
NEXT X 
LDIR 0 
LORG 2 
CSIZE 3, .6 
MOVE X, Atten 
LINE TYPE 1 

LABEL USING M K M ; Temp, M C” 

Temp=Tenp+ 1 0 
J= J+ 1 5 
NEXT L 

PRINT M HIT CONTINUE TO CONTINUE” 

PAUSE 

GCLEAR 

OUTPUT KBD; "K" ; 

ELSE 

OUTPUT KBD; "K" ; 

END IF 

Temp=Temp_con 

SUBEND 



sub Log_log* Xmin, Xnax , Ynin , Yroax , Xvar< * ) , J , Range* * ) .L.TitleS.XlabelS.Ylabel 



THIS SUBROUTINE MAKES A LOG-LOG PLOT OF DATA PASSED FROM THE 
MAIN PROGRAM. 

INPUT MINIMUM VALUES OF )* AND Y FOR PLOT 
MAXIMUM VALUES OF X AND Y FOR PLOT 
X AND Y VALUES TO BE PLOTTED 
TITLE OF PLOT 
LABELS FOR X AND Y AXIS 
OUTPUT LOG-LOG PLOT 
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VARIABLES 

Dx STEP SIZE FOR X AXIS FOR VARIOUS PLOTS. 

Dy STEP SIZE FOR Y AXIS FOR VARIOUS PLOTS. 

J FIRST ORDER INDEX FOR ASSORTED LOOPS. 

L SECOND ORDER INDEX FOR ASSORTED LOOPS. 

TitleS STRING PASSED TO SUBPROGRAM PT FOR TITLE OF PLOT. 

XlabelS LABEL ON X AXIS PASSED TO SUBPROGRAM Pt FOR PLOTTING. 

Xnax VALUE OF LARGEST VALUE OF X FOR VARIOUS PLOTS. 

Xmin VALUE OF SMALLEST VALUE OF X FOR VARIOUS PLOTS. 

Xrange VALUE OF RANGE OF X VALUES FOR VARIOUS PLOTS. 

Xvar< * ) VALUE OF ARRAY OF X VALUES TO BE PLOTTED IN SUBPROGRAM Pt . 

YlabelS LABEL ON Y AXIS PASSED TO SUBPROGRAM Pt FOR PLOTTING. 

Ynax LARGEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

Ynin SMALLEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

Yrange RANGE OF Y VARIABLES FOR VARIOUS PLOTS. 



Dx=t 
Dy= 1 

Xrange=ABS< Xnax -Xmin ) 

Yrange=ABS( Ymax-Ymin) 

WINDOW Xmin , Xnax , Ynin , Ynax 
AXES Dx ,Dy, Xmin, Ynin, 1 , 1 
CLIP OFF 

l 

! LABEL PLOT 
CSIZE 4 , . G 
LDIR 0 
LORG 5 

MOVE Xnin+ . 5*Xrange , 1 . 1 * Yrange 
T itle$=TRIM$( Title®) 

LABEL Titles 
(LABEL HORIZONTAL AXES 
LDIR 0 
LORG 5 

FOR J=Xmin TO Xnax STEP Dx 
CSIZE 4 , . 6 

MOVE J- . 01 3*Xrange 067* Yrange 
LABEL M 10 M 

MOVE J4 . 01 3*Xrange 033* Yrange 
CSIZE Z 
LABEL J 
NEXT J 

MOVE Xmin+ . 5*Xrange 1 2* Yrange 
LORG 5 
CSIZE 4 , . 6 

XI abel $=TRI M$< XI abel $ ) 

LABEL XlabelS 
(LABEL VERTICAL AXES 
LORG 8 

FOR J=0 TO Ynax STEP Dy 
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CSIZE 4, .6 

MOVE Xnin- .0Z5*Xrange , J 
LABEL M 10” 

CSIZE Z 

MOVE Xnin- .00Z5*Xrange , J+ .03* Yrange 
LABEL J 
NEXT J 
LDIR PI/Z 
LORG G 

MOVE Xnin- . 1 *Xrange 5* Yrange 
CSIZE 4, .6 

Ylabel$=TRIM$( Ylabel$) 

LABEL YlabelS 
CLIP ON 
FOR J=t TO L 

PLOT LGT< Xvar( J ) > ,LGT( Range( J ) ) 

NEXT J 
SUBEND 
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APPENDIX B 
COMPUTER PROGRAM 2 



FULLER, ROBERT PROG_Z FREQUENCY 

9 SEP 85 



♦♦♦♦♦♦♦PURPOSE******* 



THIS PROGRAM WILL TAKE PARAMETERS OF AN ECHOSOUNDER AND 
ESTIMATE THE RANGE AS A FUNCTION OF FREQUENCY. THE FOLLOWING INPUTS 
ARE REQUIRED: 

ATMOSPHERIC DATA 

1 )ATOMOSPHERIC PRESSURE in millibars 

Z ) SELECT ONE OF FOUR ATMOSPHERIC PROFILES FOR THE TEMPERATURE 
STRUCTURE PROFILE. 

a)FOR ONE PROFILE THE HEIGHT OF THE INVERSION LAYER IS ALSO INPUT 

3 ) TEMPERATURE IN DEGREES CELSIUS 

4 ) UATER VAPOR PRESSURE IN millibars 

ECHOSOUNDER DATA 

5 ) ANTENNA DIAMETER IN METERS 

6 ) POUER TRANSMITTED BY ECHOSOUNDER IN WATTS 

7 ) PULSE LENGTH OF THE TRANSMITTED ACOUSTIC ENERGY 

THE PROGRAM OUTPUTS THE FOLLOWING GRAPHS TO AN EXTERNAL PLOTTER: 

1 ) RANGE ( m. ) VERSUS FREQUENCY OF ECHOSOUNDER 

Z ) RANGE ( m. ) VERSUS EXCESS ATTENUATION FOR VARIOUS FREQUENCUIES 



***** VARIABLES***** 



Again 

Ano_change 



Ant_area 

Ant_diam 
Atom_pres 
Att en 

Bn 



C 

Cl 

C3 

CneZ < * ) 



TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 
WISHES TO MAKE ANOUTHER RUN. 

TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 
WISHES TO CHANGE THE VALUE OF INPUT DATA BEFORE MAKING 
ANOUTHER RUN. 

ECHOSOUNDER ANTENNA AREA CALCULATED FROM INPUT OF 
ANTENNA DIAMETER. 

INPUT OF ECHOSOUNDER ANTENNA DIAMETER IN METERS. 

INPUT OF ATMOSPHERIC PRESSURE IN mb. 

ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 
SUBPROGRAM ATTENUATION. 

BANDWIDTH OF RECEIVER NEEDED TO RECEIVE SIGNALS 
DOPPLER SHIFFED BY THREE METER PER SECOND 
VERTICAL VELOCITIES. 

VELOCITY OF SOUND CALCULATION AT ZERO CELSIUS. 

VARIABLE USED TO SCALE INVERSION HEIGHT IN 
TEMPERATURE STRUCTURE PROFILE THREE. 

VARIABLE USED IN CALCULATING THE FOURTH CteZ PROFILE. 

ARRAY OF VALUES OF THE ACOUSTIC REFRACTIVE INDEX 
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! N 



PARAMETER. CALCULATED BASED ON SELECTION OF 
PROFILE FOR TEMPERATURE STRUCTURE PARAMETER 
AND ASSUMED PROFILE FOR THE VELOCITY STRUCTURE 
PARAMETER. 

CteZ VALUE OF THE TEMPERATURE STRUCTURE PARAMETER. 

THE OPERATOR SELECTS A PROFILE FROM SEVERAL PROVIDED. 

CveZ VALUE OF THE VELOCITY STRUCTURE PARAMETER. 

VALUES BASED ON CALCULATION USING ASSUMED DISSAPATION 
RATE. 

Dx STEP SIZE FOR X AXIS FOR VARIOUS PLOTS. 

Dy STEP SIZE FOR Y AXIS FOR VARIOUS PLOTS. 

Epsilon DISSAPATION RATE USED IN CALCULATION OF CveZ. 

Et TRANSMISSION EFFICIENCY OF ECHOSOUNDER. 

Er EFFICIENCY OF ECHOSOUNDER WHEN RECEIVING ACOUSTIC 

BACKSCATTER. 

Exc_att( * , * ) EXCESS "ATTENUATION*' AT GIVEN RANGE. 

F VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 

RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 

Fmax FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 

ATTENUATION. 

Freq INPUT FREQUENCY OF ECHOSOUNDER. 

G ANTENNA EFFECTIVE APERATURE FACTOR. 

H VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 

I MAIN INDEX FOR LOOP TO CHECK RETURNED POWER UNTIL 

POWER WAS LESS THAN MINIMUM DETECTABLE SIGNAL. 

IND INDEX USED FOR LOOP FOR DIFFERENT FREQUENCIES 

Interval MAXIMUM EFFECTIVE- SCATTERING VOLUME THICKNESS. 

Inver HEIGHT OF INVERSION LAYER. 

J FIRST ORDER INDEX FOR ASSORTED LOOPS. 

K WAVENUMBER 

L SECOND ORDER INDEX FOR ASSORTED LOOPS. 

Mess_up TEST VALUE TO CHECK RESPONSES WHEN OPERATOR IS INPUTTING 
RESPONSES. 

TEST VALUE IN CALCULATION OF COHERENCE LENTH IN CALCULATION 
OF EXCESS ATTENUATION. 

TEST VALUE IN CHECKING IF OPERATOR WISHES TO CHANGE 
A VARIABLE BEFORE A NEW RUN. 

ASSUMED MINIMUM DETECTABLE SIGNAL. 

*) POWER BACKSCATTERED FROM GIVEN RANGE. 



Neu_va 

Noise 
Pou_back ( 
Pou_ret ( * 
Pou__trans 
Prof ile 
Pst ar 

Pulse 

R 

R1 

Range ( * ) 
Remainder 



) POWER BACKSCATTERED TO ECHOSOUNDER FROM GIVEN RANGE. 

INPUT OF POWER SUPPLIED TO TRANSDUCER OF ECHOSOUNDER. 
OPERATOR INPUT OF CteZ PROFILE FROM AVAILABLE PROFILES. 
VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 

TRANSMITTED PULSE LENGTH IN MILLISECONDS. INPUT BY 
OPERATOR. 

RANGE FOR INPUT TO INTEGRAL OF EXCESS ATTENUATION. 

VARIABLE USED IN SCALING INVERSION HEIGHT FOR PROFILE THREE. 
ARRAY OF RANGE VALUES. 

REMAINDER OF MODULO FUNCTION USED TO DECREES THE 
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! 

! 

! Rge 
! Rho 
I 


NUMBER OF PASSES THROUGH THE INTEGRATION FOR EXCESS 
ATTENUATION. 

VALUE OF RANGE IN INTEGRAL FOR EXCESS ATTENUATION. 
CORRELATION LENGTH USED IN CALCULATION OF EXCESS 
ATTENUATION. 


! Sigma 


FRACTION OF ENERGY BACKSCATTERED FROM GIVEN RANGE. 



Speed_sound SPEED OF SOUND IN AIR AT INPUT TEMPERATURE 
Sumpou_back SUM OF BACKSCATTERED ENERGY 



! T 
1 Temp 
! Temp$ 

I 


INPUT TEMPERATURE IN DEGREES KELVIN. 
INPUT TEMPERATURE IN DEGREES CELSIUS. 
VARIABLE STRING USED IN FUNCTION YES. 


! Title$ 
! Tstar 

j 


STRING PASSED TO SUBPROGRAM Semi_log FOR TITLE OF PLOT. 
INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 
IN SUBPROGRAM ATTENUATION. 


! Var 
! 


USED TO SELECT THE VARIABLE THE OPERATOR WISHED TO 
CHANGE BEFORE MAKING ANOTHER RUN. 


! Uat_pres 
! 


ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 
OPERATOR. 


! X 
! X$ 
! 


THIRD ORDER INDEX USED IN VARIOUS LOOPS. 

STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 
RESPONSE TO YES OR NO QUESTION. 


! XI abel $ 
! Xmax 
! Xmi n 
! Xrange 
! Xvar ( * ) 


LABEL ON X AXIS PASSED TO SUBPROGRAM Semi_log FOR PLOTTING. 
VALUE OF LARGEST VALUE OF X FOR VARIOUS PLOTS. 

VALUE OF SMALLEST VALUE OF X FOR VARIOUS PLOTS. 

VALUE OF RANGE OF X VALUES FOR VARIOUS PLOTS. 

VALUE OF ARRAY OF X VALUES TO BE PLOTTED IN SUBPROGRAM Semi 


! Y1 abel $ 
! Ymax 
! Ymin 
! Yrange 
! Ze 


LABEL ON Y AXIS PASSED TO SUBPROGRAM Semi_log FOR PLOTTING. 
LARGEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

SMALLEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

RANGE OF Y VARIABLES FOR VARIOUS PLOTS. 

INTERMEDIATE VALUE IN CALCULATION OF EXCESS ATTENUATION. 



! DECLARE VARIABLES 
INTEGER I 

DIM Pou_back( 1 500) , Pou_ret( 1 500) ,Range( 1 500) 

DIM CneZ ( 1500) ,Exc_att< 1500,30) ,Xvar( 1500) ,Ran< 30) 

DIM Tit 1 e$[ 503 ,Xlabel$[503 ,Ylabel$[ 1G3 
! 

PLOTTER IS 705 , " HPGL" 

LINE TYPE 1 

! 

! INPUT ATMOSPHERIC DATA 

INPUT "ENTER TEMPERATURE IN DEGREES CELCI US" , Temp 
INPUT "ENTER ATMOSPHERIC PRESSURE IN MILLIBARS" ,Atom_pres 
INPUT "ENTER WATER VAPOR PRESSURE IN MI LLI BARS" , Wat_pres 
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! INPUT ECHOSOUNDER DATA 

INPUT ’’ENTER ANTENNA DIAM IN METERS" ,Ant_dian 

INPUT "ENTER ECHOSOUNDER PULSE LENGTH IN MILLISECONDS-" .Pulse 

INPUT "ENTER POWER TO TRANSMITTER IN WATTS" ,Pou_trans 

Et*. Z5 {TRANSMIT EFFICIENCY 

Er=.Z5 ! RECEIVER EFFICIENCY 

G=.40 {ANTENNA EFFECTIVE APERATURE FACTOR 



{SELECT THE PROFILE FOR THE TEMPERATURE STRUCTURE 
PRINT "YOU HAVE A CHOICE OF TEMPERATURE STRUCTURE PROFILES" 

PRINT "PROFILE TO BE USED." 

PRINT "YOUR SELECTIONS ARE- 
PRINT " 1 

PRINT " 

PRINT " 

PRINT ” " 

PRINT " Z 

PRINT " 

PRINT " 

PRINT " " 

PRINT ’’ 3 

PRINT " 

PRINT " 

PRINT " 

PRINT " 

PRINT " " 

PRINT " 4 

PRINT " 

PRINT " " 

PRINT ” " 

Mess_up= 1 
WHILE Mess_up=1 

INPUT "ENTER THE DESIRED PROFILE (1 OR Z OR 3 OR 4)", Profile 

IF Prof ile=4 THEN 

INPUT "ENTER THE HEIGHT OF THE INVERSION IN METERS" . I nver 

END IF 

IF Profiled THEN 

Mess_up=0 

ELSE 

IF Prof i 1 e=Z THEN 

Mess_up=0 

ELSE 

IF Prof i 1 e=3 THEN 

Me5s_up=0 

ELSE 

IF Prof i 1 e=4 THEN 
Mess_up=0 

ELSE 



A TEMPERARURE STRUCTURE PROFILE BASED ON DATA AS 
PRESENTED IN WALTERS/KUNDEL 1381 PAGE 338 WHICH " 
GIVES A HEIGHT TO THE -1.16 PROFILE" 

THE TEMPERATURE STRUCTURE PROFILE ABOVE BUT" 

WITH A FACTOR OF TWO TO APPROXIMATE LOOKING" 

UP A THERMAL PLOOM" 

A TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS- 
PRESENTED IN WALTERS/KUNDEL 1381 PAGE 338 FOR" 
NIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT " 

IS EXP( - .001 ^HEIGHT ABOVE 66 METERS AND HEIGHT" 

TO THE -1.46 BELOW 65 METERS" 

A TEMPERATURE STRUCTURE PROFILE BASED ON A" 

HEIGHT TO THE -4/3 " 
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PRINT Profile," WAS NOT ONE OF THE OPTIONS MM" 
Mess_up= I 
END IF 
END IF 
END IF 
END IF 
END WHILE 
OUTPUT KBD;"K"s 
I 
! 

Again= 1 
WHILE Again=1 
FOR J=i TO I 

Pou_back ( J ) =0 
Pou_ret ( J ) =0 
Range( J ) =0 
CneZ ( J ) =0 
FOR L= 1 TO Ind 

Exc_att( J.L )=0 
NEXT L 
NEXT J 
( 

! CONVERT TEMPERATURE TO KELVIN 
T=Temp+Z73 
( 

(CALCULATE SPEED OF SOUND BASED ON INPUT TEMPERATURE( CELCIUS ) 
Speed_sound=Z0. 05* ( T ) A . 5 

(CALCULATE THE SPEED OF SOUND AT 0 DEGREES CELCIUS 
C=Z0. 05*Z73 A . 5 
! 

1=1 (INDEX FOR LOOP 
I nd= 1 

Pou_back<0>=0 (INITIALIZE VARIABLE FOR POWER BACKSCATERRED 
Ant_area=PI*< Ant_diam/Z ) A Z (ANTENNA AREA 
Intervals Speed_sound*Pul se* 1 .E-3)/Z 
( 

(INITIALIZE PLOTTING PARAMETERS SO CAN COMPARE DURING RUN TO SCALE 

(PLOTS 

Ymi n= 1 000 

Ymax=0 

(CALCULATE RETURNED SIGNAL POWER UNTIL IT IS LESS THAN BACKGROUND 
(FOR EACH FREQUENCY. 

Range( 0 )=0 

FOR Freq= 1 00 TO 3000 STEP 100 

Bn=Z*Freq* ( 1 - 1 /( 6/ Speed_sound+ 1 ) ) ! BANDWIDTH FOR 3M/S VERTICAL VELOCITY 
Noise= 1 . 38E-Z3*Bn* < T )+2 . E- 1 4 ! MINIMUM DETECTABLE SIGNAL AS 
! JOHNSON NOISE ( NEGLIGABLE )+ ESTIMATED BACKGROUND 

(CALCULATE THE ATTENUATION COEFFICIENT BASED ON ATMOSPHERIC DATA 
CALL Attenuat ion( Atom_pres , At ten , Freq .Temp , Wat_pres ) 

K=Z*PI*Freq/ Speedy sound ! WAVENUMBER 
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1-1 

Sumpou_back =0 
REPEAT 

Range( I )=Range( 1-1 ) + Z 
SELECT Profile 
CASE 1 

(THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
! AS PRESENTED IN WALTERS/ KUNDEL 1981 PAGE 398. 

! WAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
! TULAROSA BASIN, NEW MEXICO. 

IAN ADDITONAL FACTOR OF 1 . 8 WAS INCLUDED AS PER DIFFERENT 
! AVERAGING TIME. 

Cte2 = 2 . 1 2 *Range( I ) A <-1 . 1G) 

CASE 2 

(THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
IAS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398. 

! WAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
! TULAROSA BASIN, NEW MEXICO. 

IAN ADDITONAL FACTOR OF 1 . 8 WAS INCLUDED AS PER DIFFERENT 
! AVERAGING TIME AND AN ADDITIONAL FACTOR OF TWO WAS 
! INCLUDED TO APPROXIMATE LOOKING UP A THERMAL PLOOM 
CteZ = Z*2 . 1Z*Range(ir<-1. 16) 

CASE 3 

! THI S TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS” 

PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 FOR” 

! NIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT ” 

IIS EXP( - . 001 * HEIGHT ABOVE G5 METERS AND HEIGHT” 

! TO THE -1.4G BELOW 65 METERS” 

IF Range( I ) <65 THEN 

CteZ=75. 5*Range< I ) A <-1 .46) 

ELSE 

Cte2=3 . 66E-2 *EXP( - . 001 *Range( I ) ) 

END IF 
CASE 4 

! CALCULATE THE TEMPERATURE STRUCTURE FACTOR BASED ON ASSUMPTION 
! THAT IT’S PROPORTIONAL TO HEIGHT TO THE -4/3 POWER 

! THAT THE SURFACE VERTICAL HEAT FLUX IS .095 C.m./sec 

ITHE DEPENDENCE WITH HEIGHT ABOVE THE INVERSION I PULLED FROM 
! MY LEFT EAR 

! EQUATION FROM NEFF, 1975 
SELECT Range< I )/Inver 
CASE <.9 

C3=< ( . 024 ) * ( T ) * ( .667) ) 

Cte2 = C3*( Range( I ) ) A < - 1 .33) 

CASE .9 TO 1 
Ct eZ =Ct eZ 
C 1 =Cte2 
R1 =Range( I ) 

CASE 1 TO 1 .3 

CteZ= 1 0 A < (LGT( C3*Range( I ) )-LGT< C3*R1 ) ) + LGT<C1 ) ) 



R2=Range( I ) 

CZ=CteZ 
CASE ELSE 

CteZ=C3*< Range* I ) * ( - 1 . 33 )-RZ * ( - I . 33 ) >+CZ 
END SELECT 
END SELECT 

! CALCULATE VELOCITY STRUCTURE FACTOR. FORMULA FROM 6AYN0R 77 
Epsilon=< ,Z8GG/Range( I ) )*( 1 + . 07* ( Range* I ) ) * < .6))** 1 .5) 

CveZ=Z* Epsilon* ( . GG7 ) 

! CALCULATE ACOUSTIC REFRACTIVE INDEX STRUCTURE FACTOR. 

! FORMULA FROM TETARSKI 

Cne2( I ) = < CteZ/( 2-38EH5) ) + ( CveZ/< C*C) ) 

! CALCULATE THE FRACTION OF POUER BACKSCATTERED FROM INTERVAL. 

! EQUATION FROM NEFF, 1975 

Signa=* .0039** K** 1/3) )*CteZ )/( T)*Z 

j 

! CALCULATE THE EXCESS ATTENUATION. EQUATION FROM CLIFFORD, 1980 
! THE EXCESS ATTENUATION IS Exc_att< I ,L , X ) 

IF I <50 THEN 

Renainder=0 

ELSE 

IF I < 1 50 THEN 

Remainder=Range* I ) MODULO 10 

ELSE 

Renainder=Range( I ) MODULO Z0 
END IF 
END IF 

IF Renainder=0 THEN 
Rho=0 
L=0 
H=0 

Rge=Range( I ) ! CONSTANT IN INTEGRAL 

R=0 

FOR J=0 TO Z*I 

F=CneZ ( INT< J/Z+1 ) ) 

F=F* ( 1 -R/< Rge)*( 1 . G7 )+( R/ ( Rge ) ) * ( 1 .67)) 

IF J>0 THEN 

IF J<Z*I THEN 

IF INT( J/Z )= J/Z THEN 
L=L+F 
F=0 

ELSE 

H=H+F 

F=0 

END IF 
END IF 
END IF 
Rho=Rho+F 
R=R4 1 

NEXT J 
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Rho=Rho+4*L+2*H 

Rho=( ( ( Rho* . 33 ) *K*K* 1 .46) A (-.6) ) 

N®( Ant_dian/Rho) A Z 
IF NOI THEN 

Ze®1/< 1+N) 

ELSE 

Ze= 1 . 5/ ( 1+N) 

!STEP OF 1 . 5==>SEE CLIFFORD 1980 
END IF 

Exc_att( I ,Ind)=Ze*Ze 

ELSE 

Exc_att ( I , I nd )=Exc_att ( 1-1 , Ind) 

END IF 
{ 

! CALCULATE THE POUER BACKSCATTERED 

Pou_back ( I ) = ( Pou_trans*Et-Sumpou_back ) *EXP( -At ten* Range( I ) ) 

Pou_back( I )=Pou_back( I )* Interval* Ex c_att ( I , Ind ) * Sigma 
Sumpou_back=5umpou_back+Pou_back ( I ) 

! CALCULATE THE POUER BACKSCATTERED TO THE ANTENNA 

Pou_ret( I ) *P outback ( I )*EXP( - At ten* Range ( I ) )* Ant _ area* G* Er/ Range ( I ) A 2 
PRINT "RANGE®" ,Range( I ) 

PRINT "POUER RETURNED®" ,Pou_ret( I ) 

1=1 + 1 

UNTIL Pou_ret( 1-1 ) >0 AND Pou_ret ( I - 1 ) <=Noise 
Ran( I nd )® Range ( I -2 ) 

IF RangeC I -2 ) > Ymax THEN 
Ymax=Range( 1-2 ) 

END IF 

IF Range ( I-ZKYmin THEN 
Ymin=Range( I -2 ) 

END IF 
Ind=Ind+1 
Again® 1 
FOR J®1 TO I 

Pou_back ( J )=0 
Pou_ret ( J ) =0 
Range( J )=0 
Cne2( J ) =0 
NEXT J 

PRINT "FREQUENCY®" ,Freq 
NEXT Freq 
I 
| 

OUTPUT KBD; "K" ; 

PRINT "INPUT CONDITIONS" 

I 

PRINT " " 

PRINT USING "K"; "TEMPERATURE® ", Temp , "CELCI US- 

PRINT USING "K" ; "ATMOSPHERIC PRESSURE® " , At on_pres , " nb" 

PRINT USING "K" ; "UATER VAPOR PRESSURE® " , Uat_pres , " nb" 
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PRINT USING H K n ; "PULSE LENGTH^ M , Pulse f "ms" 

PRINT USING "K"; "ANTENNA DI RMETER= " f Ant_diam, " m." 

PRINT USING "K" ; "POWER TRANSMI TTED= " f Pou.tr ans , M WATTS" 

PRINT "TEMPERATURE STRUCTURE PROFILE USED ", Profile 
IF Profile“4 THEN 

PRINT "INVERSION HEIGHT ", Inver 
END IF 

SELECT Profile 
CASE 1 

PRINT " 1 

PRINT " 

PRINT M 
CASE 2 

PRINT " 2 

PRINT " 

PRINT " 

PRINT " 

PRINT " 

CASE 3 

PRINT " 3 

PRINT " 

PRINT " 

PRINT " 

PRINT " 

CASE 4 

PRINT " 4 

PRINT " 

END SELECT 
! 

PRINT USING "K" ; "THE MINIMUM DETECTABLE SIGNAL WAS SET AT", Noise , "WATTS" 
| 

PRINT " " 

! 

I 

PRINT "HIT CONTINUE TO CONTINUE" 

PAUSE 

OUTPUT KBD; "K" ; 

{ 

GRAPHICS ON 
VIEWPORT 15,120, 10,70 



! PLOT RANGE VERSUS FREQUENCY 

Xmi n-0 

Xmax=3000 

Ymin=PROUND< Ymin-50,2 ) 
Ymax=PROUND< Ymax+50,2 ) 

Dx=500 

Dy=100 

Xrange^Xmax-Xmi n 



A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 
AS PRESENTED IN UALTERS/KUNDEL 1981 PAGE 398 " 
WHICH GIVES A HEIGHT TO THE -1.16 PROFILE" 

A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 
AS PRESENTED IN UALTERS/KUNDEL 1981 PAGE 398 " 
WHICH GIVES A HEIGHT TO THE -1.16 PROFILE" 

BUT WITH A FACTOR OF TWO TO APPROXIMATE" 
LOOKING UP A THERMAL PLOOM" 

A TEMPERATURE STRUCTURE PROFILE BASED ON DATA" 
AS PRESENTED IN UALTERS/KUNDEL 1981 PAGE 398" 
FOR NIGHT CONDITIONS. THE DEPENDENCE WITH" 
HEIGHT IS EXP< - . 001 * HEIGHT ABOVE 65 METERS" 

AND HEIGHT TO THE -1.46 BELOW 65 METERS" 

A TEMPERATURE STRUCTURE PROFILE BASED ON A" 
HEIGHT TO THE -4/3 " 
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Yrange=Ynax-Ynin 
WINDOW Xnin,Xnax , Ynin, Ynax 
AXES Dx ,Dy ,Xnin,Ynin,Z , 1 
CLIP OFF 
! 

! LABEL PLOT 
CSIZE 4,. 6 
LDIR 0 
LORG 5 

MOVE . 5*Xrange , Ynax* 1 . 05 
LABEL “RANGE VERSUS FREQUENCY” 

! LABEL HORIZONTAL AXES 
LDIR 0 
LORG 5 

FOR J*0 TO Xnax STEP Dx 
CSIZE 4, .6 

MOVE J , Ynin- . 05* Yrange 
LABEL J 
NEXT J 

MOVE „5*Xrange,Ynin- . 1*Yrange 

CSIZE 4 , . G 

LABEL “FREQUENCY" 

! LABEL VERTICAL AXES 
LORG 8 

FOR J=Ynin TO Ynax STEP Dy 
CSIZE 4,. 6 

MOVE Xnin- . 00Z5*Xrange , J 
LABEL USING "K" ;J 
NEXT J 
LDIR PI/Z 
LORG G 

MOVE Xnin- « 1 5*Xrange , . 5* Yrange + Ynin 
CSIZE 4 , . 6 

LABEL "RANGE (METERS)" 

CLIP ON 
Ind-1 

FOR Freq= 1 00 TO 3000 STEP 100 
PLOT Freq,Ran( Ind) 

Ind=Ind+l 
NEXT Freq 
! 

PRINT "HIT CONTINUE TO CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT KBD; "K" ; 



! PLOT RANGE VERSUS EXCESS ATTENUATION 
! 

Xnin=0 
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Ymax^ 

Xmax = 1 
Ymin=0 
Range( 0) =0 

Title$="RANGE vs EXCESS ATTENUATION" 

XI abelS 3 " EXCESS ATTENUATION" 
Ylabel$="RANGE (METERS)" 

Dx= . 1 
Dy= 1 

Xrange=ABS( Xmax-Xmin) 

Yrange=ABS( Ymax- Ymi n ) 

UINDOU Xmin , Xmax , Ymi n , Ymax 
AXES Dx , Dy , Xmin , Ymin , 1 , 1 
CLIP OFF 

I 

! LABEL PLOT 
CSIZE 4, .6 
LDIR 0 
LORG 5 

MOVE Xmin+ . 5*Xrange , 1 . 1 * Yrange 
Titl e$=TRIM$< Title® ) 

LABEL Title$ 

! LABEL HORIZONTAL AXES 
LDIR 0 
LORG 5 

FOR J=Xmin TO Xmax STEP Dx 
MOVE J 033* Yrange 
CSIZE 4 , . G 
LABEL J 
NEXT J 

MOVE Xmin+.5*Xrange,-. 12*Yrange 
LORG 5 
CSIZE 4, .6 

Xlabel$=TRIM$( XlabelS) 

LABEL XlabelS 
[LABEL VERTICAL AXES 
LORG 8 

FOR J=0 TO Ymax STEP Dy 
CSIZE 4, .6 

MOVE Xmin-.0Z5*Xrange, J 
LABEL "10" 

CSIZE Z 

MOVE Xmin- . 00Z5*Xrange , J+ . 03* Yrange 
LABEL J 
NEXT J 
LDIR PI/Z 
LORG G 

MOVE Xmin- . 1 *Xrange 5* Yrange 
CSIZE 4, .6 

Ylabel $=TRIM$( YlabelS) 
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LABEL Ylabel$ 

FOR L=5 TO Ind-1 STEP 5 
X-1 



WHILE Exc_att(X,L>>0 

Xvar( X>=Exc_att( X,L ) 
Range( X)=Range( X- 1 )+2 
X=X+ 1 
END WHILE 

IF INT(L/2)«L/Z THEN 
LINE TYPE 1 

ELSE 



LINE TYPE 3 
END IF 
Freq=L* 1 00 

CALL Semi_log( Freq,Xvar< * ) ,X-1 , Range( * ) , X- I ) 

NEXT L 

PRINT M HIT CONTINUE TO CONTINUE" 

PAUSE 

! 

GCLEAR 

OUTPUT KBD; "K" ; 

LINPUT "WOULD YOU LIKE TO MAKE ANOTHER RUN (ENTER Y OR N)?" ,X$ 
Again=FNYes( X$> 



Mess_up= 1 
WHILE Mess_up= 1 
SELECT Again 
CASE 1 

LINPUT "WOULD YOU LIKE TO CHANGE A VARIABLE (ENTER Y OR N >?" , X$ 
Neu_va=FNYes( X$> 

Ano_change= 1 
WHILE Ano_change= 1 
SELECT Neu_va 
CASE 1 

PRINT " VARIABLE 

PRINT 
PRINT 
PRINT 
PRINT 



CURRENT VALUE- 



USING "K* 


VI 


TEMPERATURE 


Temp. "CELSIUS" 


USING "K' 


’ ; "Z 


ATOMOSPHERIC PRESSURE 


" , At on_pres , " nb" 


USING "K* 

M it 


' ; ”3 


WATER VAPOR PRESSURE 


" , Wat_pres , "mb" 


USING "K* 


‘ ; "4 


ANTENNA DIAMETER 


" , Ant_diam , " m. " 


USING "K 1 


' ; "5 


PULSE LENGTH 


" , Pul se , " ms" 


USING "K* 


' 5 ”G 


POWER TRANSMITTED " , Poujrans , " WATTS" 


USING " K ' 


‘ 5 "7 


ATMOSPHERIC PROFILE 


" , Prof ile 



PRINT " " 

PRINT "ENTER THE NUMBER OF THE VARIABLE YOU WISH TO" 
PRINT "CHANGE- 
INPUT Var 
SELECT Var 

CASE 1 

INPUT "TEMPERATURE=" , Temp 
CASE 2 
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INPUT 
CRSE 3 
INPUT 
CRSE 4 
INPUT 
CRSE 5 
INPUT 
CRSE G 
INPUT 
CRSE 7 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 
IF Prof 



"RTMOSPHERIC PRESSURE IN mb=" ,Rton _pres 

"WRTER VRPOR PRESSURE IN mb=",Wat jires 

"RNTENNR DIRMETER IN n. =” , Rnt_dian 

"PULSE LENGTH IN ns=‘\Pulse 

•‘POWER TRRNSMITTED IN WRTTS=” f Pow_trans 

1 ==> Ct A Z PROFILE OF Z*<-1.1G>" 

FROM URLTERS/KUNDEL 1981 n 
Z ==> SRME RS ONE BUT UITH FRCTOR” 

OF TUO TO RPPROXIMRTE LOOKING*’ 

UP R THERMRL PLOOM*' 

3 ==> Ct A Z PROFILE OF EXP< - . 001 *2 ) ” 

FROM URLTERS/KUNDEL 1981.” 

4 ==> CT A Z PROFILE OF Z A (-4/3) M 

ENTER NUMBER OF DESIRED PROFILE” , Prof i 1 e 
ile-4 THEN 



INPUT '‘HEIGHT OF INVERSION IN METERS=” , Inver 
END IF 
CRSE ELSE 

PRINT Var , ” I S NOT ONE OF THE OPTIONS” 

END SELECT 

LINPUT "MRDE RNOTHER CHRNGE( ENTER Y OR N)?”,X$ 
Neu_va=FNYes( X$) 

Mess_up=Z 
CRSE Z 



Mess_up=2 
Rno_change=Z 
CRSE ELSE 

PRINT "YOUR CHOICES WHERE Y OR N MM!” 
Rno_change= 1 

LINPUT "CHRNGE R VRRIRBLE (ENTER Y OR N)?”,X$ 
Neu_va=FNYes( X$) 

END SELECT 
END WHILE 
CRSE Z 

Mess_up=Z 
CRSE ELSE 

PRINT "YOUR CHOICES WHERE Y OR N ! M ! ! ! ” 

Mess_up= 1 

LINPUT "MRKE RNOTHER RUN (ENTER Y OR N)?”,X$ 

Rgain=FNYes( X$) 

END SELECT 
END WHILE 
END WHILE 

PRINT “THRT’S RLL , FOLKS” 
end 
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ICALCULTE THE ATTENUATION 

sub Attenuat ion( Atom_pres f At ten , Freq , Temp , Uat_pres ) 
ITHIS SUBPROGRAM CALCULATES THE ATTENUATION OF SOUND 
! I N AIR BASED UPON EQUATIONS IN NEFF 1975 



INPUT ATMOSPHERIC PRESSURE IN MILLIBARS 
FREQUENCY OF SOUND WAVE IN HERTZ 
TEMPERATURE IN DEGREES CELCIUS 
WATER- VAPOR PRESSURE IN MILLIBARS 



OUTPUT ATTENUATION IN 1 /METERS 



VARIABLES 

Atom_pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 

Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 

SUBPROGRAM ATTENUATION. 

Att_max VARIABLE IN SUBPROGRAM ATTENUATION. IT IS THE 
ATTENUATION AT THE FREQUENCY OF THE MAXIMUM 
ATTENUATION FOR THE INPUT CONDITIONS. 

F VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 

RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 

Fmax FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 

ATTENUATION. 

Freq INPUT FREQUENCY OF ECHOSOUNDER. 

H VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 

Pstar VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 

CALCULATION. 

Temp INPUT TEMPERATURE IN DEGREES CELSIUS. 

Tstar INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 

IN SUBPROGRAM ATTENUATION. 

Wat_pres ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 
OPERATOR. 



H= 1 00* Wat_pres/Atom_pres 
Tstar=< 1 . 8*Temp+492 >/51 9 
Pstar=Atom_pres/ 1 01 4 

Fmax=( 1 0+6600* H+44400*H*H ) *Pstar /Tstar A . 8 
Att_max=.0078*Fmax*Tstar A < - Z . 5 > *EXP< 7 . 77* < 1-1/Tstar ) ) 

F=Freq/Fmax 

Atten=( Att_max/304. 8 )*( ( . 1 8*F ) A Z+< 2*F*F/ ( 1+F*F) ) A Z) A .5 
Atten=( Atten+ 1 . 74E- 1 0* Freq* Freq )/ 4. 35 
SUBEND 
! 

def FNYes( X$ ) 

! THIS FUNCTION INTERPRETS THE OPERATORS RESPONSES TO YES NO QUESTIONS 
! 

! INPUT X$ 

! 

! OUTPUT FNYes 
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! VARIABLES 

! Temp$ VRRIRBLE STRING USED IN FUNCTION YES. 

! X$ STRING PRSSED TO FUNCTION YES TO DETERMINE OPERRTOR 

! RESPONSE TO YES OR NO QUESTION. 

DIM Temp$[ 13 
T emp$[ 1 , 13=TRIM$(X$) 

SELECT Tenp$ 

CRSE 'T/'y" 

RETURN 1 
CRSE M N ,, , ,, n M 
RETURN Z 
CRSE M M 

RETURN 1 
CRSE ELSE 

RETURN -Z 
END SELECT 

FNEND 

! 

! 

sub Semi_log( Freq,Xvar< * ) , L , Range< * ) , J ) 
j 

! THIS SUBROUTINE MRKES R SEMI-LOG PLOT OF DRTR PRSSED FROM THE 

IMRIN PROGRRM. 

I 

! INPUT 

! X RND Y VRLUES TO BE PLOTTED 

! OUTPUT SEMI-LOG PLOT 

i 

! VRRIRBLES 

l J FIRST ORDER INDEX FOR RSSORTED LOOPS. 

! L SECOND ORDER INDEX FOR RSSORTED LOOPS. 

! Xvar<*) VRLUE OF RRRRY OF X VRLUES TO BE PLOTTED IN SUBPROGRRM Se 
! 

CLIP ON 

FOR J»1 TO L STEP 1 

PLOT Xvar ( J ) , LGT ( Range< J ) ) 

NEXT J 
LDIR PI/4 
LORG Z 
CSIZE 3, .6 
LINE TYPE 1 

MOVE Xvar( J- 1 ) , LGT ( Range ( J- 1 ) ) 

LRBEL USING n K" ; Freq , M Hz M 
SUBEND 
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APPENDIX 



COMPUTER PROGRAM 3 

! FULLER, ROBERT PR0G_3 EFFICIENCY 

! 

! 1 0 SEP 85 



j *******PURP0SE**» »**» 

! THIS PROGRAM WILL TAKE PARAMETERS OF AN ECHOSOUNDER AND 

! ESTIMATE THE RANGE AS A FUNCTION OF EFFICIENCY OF THE THE TRANSDUCER. 
ITHE FOLLOWING INPUTS ARE REQUIRED: 

! ATMOSPHERIC DATA 

l 1 )ATOMOSPHERIC PRESSURE in millibars 

! 2 ) SELECT ONE OF FOUR ATMOSPHERIC PROFILES FOR THE TEMPERATURE 

! STRUCTURE PROFILE. 

! a)FOR ONE PROFILE THE HEIGHT OF THE INVERSION LAYER IS ALSO INPUT 

! 3 ) TEMPERATURE IN DEGREES CELSIUS 

! 4 ) WATER VAPOR PRESSURE IN millibars 

! ECHOSOUNDER DATA 

! 5 ) ANTENNA DIAMETER IN METERS 

! G) FREQUENCY 

! 7 )POWER TRANSMITTED BY ECHOSOUNDER IN WATTS 

! 8 )PULSE LENGTH OF THE TRANSMITTED ACOUSTIC ENERGY 

ITHE PROGRAM OUTPUTS THE FOLLOWING GRAPH TO AN EXTERNAL PLOTTER: 

! 1 ) RANGE ( m. ) VERSUS EFFICIENCY OF ECHOSOUNDER 



* **** VARI ABLES**** * 



Again TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 

WISHES TO MAKE ANOUTHER RUN. 

Ano_change TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 

WISHES TO CHANGE THE VALUE OF INPUT DATA BEFORE MAKING 
ANOUTHER RUN. 

Ant_area ECHOSOUNDER ANTENNA AREA CALCULATED FROM INPUT OF 
ANTENNA DIAMETER. 

Ant_diam INPUT OF ECHOSOUNDER ANTENNA DIAMETER IN METERS. 

AtomjDres INPUT OF ATMOSPHERIC PRESSURE IN mb. 

Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 

SUBPROGRAM ATTENUATION. 

Bn BANDWIDTH OF RECEIVER NEEDED TO RECEIVE SIGNALS 

DOPPLER SHIFFED BY THREE METER PER SECOND 
VERTICAL VELOCITIES. 

C VELOCITY OF SOUND CALCULATION AT ZERO CELSIUS. 

Cl VARIABLE USED TO SCALE INVERSION HEIGHT IN 

TEMPERATURE STRUCTURE PROFILE THREE. 

C3 VARIABLE USED IN CALCULATING THE FOURTH CteZ PROFILE. 

Cne2 ( * ) ARRAY OF VALUES OF THE ACOUSTIC REFRACTIVE INDEX 



CteZ 


PARAMETER. CALCULATED BASED ON SELECTION OF 
PROFILE FOR TEMPERATURE STRUCTURE PARAMETER 
AND ASSUMED PROFILE FOR THE VELOCITY STRUCTURE 
PARAMETER. 

VALUE OF THE TEMPERATURE STRUCTURE PARAMETER. 

THE OPERATOR SELECTS A PROFILE FROM SEVERAL PROVIDED. 


CveZ 


VALUE OF THE VELOCITY STRUCTURE PARAMETER. 

VALUES BASED ON CALCULATION USING ASSUMED DISSAPATION 
RATE. 


Dx 

□y 

Epsilon 

Eff 

Eff 


STEP SIZE FOR X AXIS FOR VARIOUS PLOTS. 

STEP SIZE FOR Y AXIS FOR VARIOUS PLOTS. 
DISSAPATION RATE USED IN CALCULATION OF CveZ. 
TRANSMISSION EFFICIENCY OF ECHOSOUNDER. 

EFFICIENCY OF ECHOSOUNDER WHEN RECEIVING ACOUSTIC 
BACKSCATTER. 


Exc_att 

F 


EXCESS “ATTENUATION” AT GIVEN RANGE. 

VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 


Fnax 


FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
ATTENUATION. 


Freq 

G 

H 

I 


INPUT FREQUENCY OF ECHOSOUNDER. 

ANTENNA EFFECTIVE APERATURE FACTOR. 

VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 
MAIN INDEX FOR LOOP TO CHECK RETURNED POWER UNTIL 
POWER WAS LESS THAN MINIMUM DETECTABLE SIGNAL. 


IND 

I nterval 

Inver 

J 

K 

L 

Mess_up 


INDEX USED FOR LOOP FOR DIFFERENT EFFICIENCIES 
MAXIMUM EFFECTIVE-SCATTERING VOLUME THICKNESS. 

HEIGHT OF INVERSION LAYER. 

FIRST ORDER INDEX FOR ASSORTED LOOPS. 

WAVENUMBER 

SECOND ORDER INDEX FOR ASSORTED LOOPS. 

TEST VALUE TO CHECK RESPONSES WHEN OPERATOR IS INPUTTING 
RESPONSES. 


N 


TEST VALUE IN CALCULATION OF COHERENCE LENTH IN CALCULATION 
OF EXCESS ATTENUATION. 


Neu__va 


TEST VALUE IN CHECKING IF OPERATOR WISHES TO CHANGE 
A VARIABLE BEFORE A NEW RUN. 



Noise ASSUMED MINIMUM DETECTABLE SIGNAL. 

Pow_back ( * ) POWER BACKSCATTERED FROM GIVEN RANGE. 

Pou_ret( * ) POWER BACKSCATTERED TO ECHOSOUNDER FROM GIVEN RANGE. 



Pou_trans 

Profile 

Pstar 


INPUT OF POWER SUPPLIED TO TRANSDUCER OF ECHOSOUNDER. 
OPERATOR INPUT OF CteZ PROFILE FROM AVAILABLE PROFILES. 
VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 


Pulse 


TRANSMITTED PULSE LENGTH IN MILLISECONDS. INPUT BY 
OPERATOR. 


R 

R1 

Range( * ) 
Remainder 


RANGE FOR INPUT TO INTEGRAL OF EXCESS ATTENUATION. 

VARIABLE USED IN SCALING INVERSION HEIGHT FOR PROFILE THREE 
ARRAY OF RANGE VALUES. 

REMAINDER OF MODULO FUNCTION USED TO DECREES THE 
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! NUMBER OF PASSES THROUGH THE INTEGRATION FOR EXCESS 

! ATTENUATION. 

! Rge VALUE OF RANGE IN INTEGRAL FOR EXCESS ATTENUATION. 

! Rho CORRELATION LENGTH USED IN CALCULATION OF EXCESS 

! ATTENUATION. 

! Sigma FRACTION OF ENERGY BACKSCATTERED FROM GIVEN RANGE. 

! Speed_sound SPEED OF SOUND IN AIR AT INPUT TEMPERATURE. 

! Sumpou_back 

! T INPUT TEMPERATURE IN DEGREES KELVIN. 

! Temp INPUT TEMPERATURE IN DEGREES CELSIUS. 

! Temp$ VARIABLE STRING USED IN FUNCTION YES. 



Tstar INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 

IN SUBPROGRAM ATTENUATION. 

Var USED TO SELECT THE VARIABLE THE OPERATOR WISHED TO 

CHANGE BEFORE MAKING ANOTHER RUN. 

Wat_pres ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 
OPERATOR. 

X THIRD ORDER INDEX USED IN VARIOUS LOOPS. 

X$ STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 

RESPONSE TO YES OR NO QUESTION. 

Xmax VALUE OF LARGEST VALUE OF X FOR VARIOUS PLOTS. 

Xmin VALUE OF SMALLEST VALUE OF X FOR VARIOUS PLOTS. 

Xrange VALUE OF RANGE OF X VALUES FOR VARIOUS PLOTS. 

YIabel$ LABEL ON Y AXIS PASSED TO SUBPROGRAM Semi_log FOR PLOTTING. 
Ymax LARGEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

Ymin SMALLEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

Yrange RANGE OF Y VARIABLES FOR VARIOUS PLOTS. 

Ze INTERMEDIATE VALUE IN CALCULATION OF EXCESS ATTENUATION. 



(DECLARE VARIABLES 
INTEGER I 

DIM Pou_back ( 150®) ,Pou_ret( 1500) ,Range( 1500) 

DIM CneZ< 1500) ,Ran< 30) 

! 

PLOTTER IS 705 , "HPGL" 

LINE TYPE 1 
! 

(INPUT ATMOSPHERIC DATA 

INPUT "ENTER TEMPERATURE IN DEGREES CELCI US" , Temp 
INPUT "ENTER ATMOSPHERIC PRESSURE IN MILLIBARS" , Atom_pres 
INPUT "ENTER WATER VAPOR PRESSURE IN MI LLI BARS" , Wat_pres 



(INPUT ECHOSOUNDER DATA 

INPUT "ENTER ANTENNA DIAM IN METERS" ,Ant_diam 
INPUT "ENTER ECHOSOUNDER FREQUENCY IN HERTZ", Freq 
INPUT "ENTER ECHOSOUNDER PULSE LENGTH IN MILLISECONDS=" , Pul se 
INPUT "ENTER POWER TO TRANSMITTER IN WATTS" ,Pou_trans 
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G-.40 ! ANTENNA EFFECTIVE APERATURE FACTOR 

I 
! 

! SELECT THE PROFILE FOR THE TEMPERATURE STRUCTURE 
PRINT “YOU HAVE A CHOICE OF TEMPERATURE STRUCTURE PROFILES** 

PRINT "PROFILE TO BE USED." 

PRINT "YOUR SELECTIONS ARE" 

PRINT " 1 

PRINT " 

PRINT " 

PRINT " " 

PRINT " Z 

PRINT " 

PRINT " 

PRINT " " 

PRINT " 3 

PRINT " 

PRINT " 

PRINT " 

PRINT " 

PRINT " " 

PRINT " 4 

PRINT " 

| 

PRINT " M 
PRINT " '* 

Mess_up= 1 
WHILE Mess_up= 1 

INPUT "ENTER THE DESIRED PROFILE (1 OR Z OR 3 OR 4>" .Profile 

IF Profile-4 THEN 

INPUT "ENTER THE HEIGHT OF THE INVERSION IN METERS" , Inver 

END IF 

IF Profile-1 THEN 
Mess_up=0 

ELSE 

IF Prof i 1 e=Z THEN 
Mess_up=0 

ELSE 

IF Profile-3 THEN 
Mess__up=0 

ELSE 

IF Profile-4 THEN 

Mess__up=0 

ELSE 

PRINT Profile," WAS NOT ONE OF THE OPTIONSM!!" 

Mess_up= 1 
END IF 
END IF 
END IF 
END IF 



A TEMPERARURE STRUCTURE PROFILE BASED ON DATA AS " 
PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 WHICH " 
GIVES A HEIGHT TO THE -1.16 PROFILE" 

THE TEMPERATURE STRUCTURE PROFILE ABOVE BUT" 

WITH A FACTOR OF TWO TO APPROXIMATE LOOKING" 

UP A THERMAL PLOOM" 

A TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS" 
PRESENTED IN WALTERS/KUNDEL 1381 PAGE 398 FOR" 
NIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT " 

IS EXP( - . 001 * HEIGHT ABOVE 65 METERS AND HEIGHT" 

TO THE -1.46 BELOW 65 METERS" 

A TEMPERATURE STRUCTURE PROFILE BASED ON A" 

HEIGHT TO THE -4/3 " 
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END WHILE 
OUTPUT KBD; M K " ; 

! 

! 

Again=1 
WHILE Again” 1 
FOR J-1 TO I 

P outback ( J ) ”0 
Pou_ret( J)=0 
Range( J )=0 
CneZ( J ) e 0 
NEXT J 
I 

! CONVERT TEMPERATURE TO KELVIN 
T=Temp+273 



! CALCULATE SPEED OF SOUND BASED ON INPUT TEMPERATURE( CELCIUS ) 
Speed_sound=Z0. 05-* ( T ) * . 5 

I CALCULATE THE SPEED OF SOUND AT 0 DEGREES CELCIUS 
C = Z0. 05*Z73 A . 5 
! 

1=1 ! INDEX FOR INNER LOOP (POWER RETURNED TO NOISE) 

I nd= 1 ! INDEX FOR OUTER LOOP (EFFICIENCY) 

Pou_back ( 0 )=0 UNITIALIZE VARIABLE FOR POWER BACKSC ATERRED 
Ant_area=PI*( Ant_dian/Z ) A Z ! ANTENNA AREA 
I nterval=( Speed_sound* Pulse* 1 . E-3 )/Z 
! 

UNITIALIZE PLOTTING PARAMETERS SO CAN COMPARE DURING RUN TO SCALE 
! PLOTS 
Ymin= 1 000 
Ymax=0 

! CALCULATE RETURNED SIGNAL POWER UNTIL IT IS LESS THAN BACKGROUND 
IFOR EACH FREQUENCY. 

Range( 0 )=0 

FOR Ef f = . 05 TO .5 STEP .05 

Bn=Z*Freq*( 1 - 1 /( G/Speed_sound4 1 ) ) ! BANDWIDTH FOR 3M/S VERTICAL VELOCITY 
Noise* 1 . 38E-Z3*Bn*( T )+Z . E- 1 4 ! MINIMUM DETECTABLE SIGNAL AS 
! JOHNSON NOISE( NEGLI GABLE )+ ESTIMATED BACKGROUND 

! CALCULATE THE ATTENUATION COEFFICIENT BASED ON ATMOSPHERIC DATA 
CALL Attenuat ion( At on_pres , Att en , Freq , Tenp , Wat_pres ) 

K=Z*PI *Freq/Speed_sound ! WAVENUMBER 
1 = 1 

Sunpou_back =0 
REPEAT 

Range( I ) = Range( I- 1 ) + Z 
SELECT Profile 
CASE I 

ITHIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
! AS PRESENTED IN UALTERS/KUNDEL 1981 PAGE 398. 
iwns FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
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! TULAROSA BASIN , NEU MEXICO. 

! AN ADDITONAL FACTOR OF 1 . 8 UAS INCLUDED AS PER DIFFERENT 
! AVERAGING TIME. 

Cte2=2. 1 2*Range( I ) A <-1 . 1G) 

CASE 2 

! THI S TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
! AS PRESENTED IN UALTERS/KUNDEL 1981 PAGE 388. 

! UAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
! TULAROSA BASIN, NEU MEXICO. 

! AN ADDITONAL FACTOR OF 1 . 8 UAS INCLUDED AS PER DIFFERENT 
! AVERAGING TIME AND AN ADDITIONAL FACTOR OF TUO UAS 
! INCLUDED TO APPROXIMATE LOOKING UP A THERMAL PLOOM 
Cte2=2*2 . 12*Range(ir<-1 . 16) 

CASE 3 

ITHIS TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS" 

PRESENTED IN UALTERS/KUNDEL 1981 PAGE 398 FOR M 
! NIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT " 

IIS EXP( - . 001 * HEIGHT ABOVE G5 METERS AND HEIGHT" 

! TO THE -1.46 BELOU 65 METERS" 

IF Range( I ) <65 THEN 

Ct e2=75 . 5*Range( I ) A < -1 .46) 

ELSE 

Ct e2 = 3. 66E-2*EXP( - . 001 * Range ( I ) ) 

END IF 
CASE 4 

! CALCULATE THE TEMPERATURE STRUCTURE FACTOR BASED ON ASSUMPTION 
! THAT IT’S PROPORTIONAL TO HEIGHT TO THE -4/3 POWER 

! THAT THE SURFACE VERTICAL HEAT FLUX IS .095 C.m./sec 

! THE DEPENDENCE WITH HEIGHT ABOVE THE INVERSION I PULLED FROM 
! MY LEFT EAR 

! EQUATION FROM NEFF, 1975 
SELECT Range< I )/I nver 
CASE <.9 

C3=< ( . 024 ) * ( T ) * ( .667)) 

Ct e2= s C3*( Range ( I ) ) A ( - 1 . 33 ) 

CASE .9 TO 1 
Cte2=Cte2 
C 1 =Cte2 
R1 =Range( I ) 

CASE 1 TO 1.3 

Ct e2= 1 0 A ( ( LGT ( C3*Range< I ) )-LGT( C3*R1 ) )+LGT< Cl ) ) 
R2=Range( I ) 

C2=Ct e2 
CASE ELSE 

Cte2=C3*( Range( I ) A < - 1 . 33 )-R2 A < - 1 . 33 ) ) + C2 
END SELECT 
END SELECT 

! CALCULATE VELOCITY STRUCTURE FACTOR. FORMULA FROM GAYNOR 77 
Epsllon=( . 2866/Range( I ) )*( 1 + . 07* < Range( I ) ) A < ,6)) A ( 1.5) 
Cve2=2*Epsilon A < .667) 
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! CALCULATE ACOUSTIC REFRACTIVE INDEX STRUCTURE FACTOR. 

•FORMULA FROM TETARSKI 

Cne2( I )*( CteZ/ ( 2 . 98E+5 ) ) + ( CveZ/ ( C*C ) ) 

(CALCULATE THE FRACTION OF POUER BACKSCATTERED FROM INTERVAL. 
(EQUATION FROM NEFF, 1375 
Sigma=( .00 39*<K A < 1 /3 ) )*CteZ )/< T) A Z 
! 

•CALCULATE THE EXCESS ATTENUATION. EQUATION FROM CLIFFORD, 1980 
(THE EXCESS ATTENUATION IS Exc_att 
IF I <50 THEN 

Remainder*© 

ELSE 

IF I < 1 50 THEN 

Remainder*Range( I ) MODULO 10 

ELSE 

Remainder=Range( I ) MODULO Z0 
END IF 
END IF 

IF Remainder*© THEN 
Rho*0 
L=0 
H=0 

Rge=Range( I ) (CONSTANT IN INTEGRAL 
R=0 

FOR J=0 TO Z*I 

F=CneZ< INT( J/Z+1 ) ) 

F=F*( 1-R/(Rge) A ( 1 . G7 ) + ( R/ < Rge ) ) A < 1 .57)) 

IF J>0 THEN 

IF J<Z«I THEN 

IF INT ( J/Z )= J/Z THEN 
L=L + F 
F=0 

ELSE 

H=H+F 

F=0 

END IF 
END IF 
END IF 
Rho=Rho+F 
R=R+ 1 

NEXT J 

Rho=Rho+4*L+Z*H 

Rho=( ( (Rho*.33)*K*K*1 .4G) A (-.G> > 

N=( Ant_diam/ Rho ) A Z 
IF N<= 1 THEN 

Ze= 1 / ( 1+N) 

ELSE 

Ze= 1 . 5/ ( 1+N) 

•STEP OF 1 . 5*=>SEE CLIFFORD 1980 
END IF 
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Exc_att=Ze*Ze 

ELSE 

Exc_att=Exc_att 
END IF 
| 

! CALCULATE THE POWER BACKSCATTERED 

Pou_back ( I ) -( Pou_trans*Ef f-Sumpou_back )*EXP< - Att en*Range( I ) ) 
Pou_back( I )=PotJ_back( I ) ^Interval *Exc_ att*Sigma 
Sumpou_back = SumpoiJ_back+Pou_back ( I ) 

! CALCULATE THE POWER BACKSCATTERED TO THE ANTENNA 
Pou_ret( I )=Pou_back( I )*EXP( -At*ten* Range ( I ) )* Ant_ar ea*G*Ef f /Range ( I ) A Z 
PRINT M RANGE=" , Range( I ) 

PRINT "POWER RETURNED^" ,Pou_ret( I ) 

1 = 1 + 1 

UNTIL Pou_ret( I- 1 ) >0 AND Pou_ret( I- 1 ) <=Noise 
Ran( Ind)=Range( I-Z ) 

IF Range( I - Z ) >Ynax THEN 
Ynax=Range( I-Z ) 

END IF 

IF Range( I-Z KYmin THEN 
Ymin=Range( I-Z ) 

END IF 
I nd=I nd+ 1 
Again= 1 
FOR J=1 TO I 

Pou_back ( J ) =0 
Pou_ret ( J )=0 
Range( J ) -<d 
CneZ ( J )=0 
NEXT J 
NEXT Eff 
! 
l 

OUTPUT KBD; "K" ; 

PRINT "INPUT CONDITIONS" 

{ 

PRINT " " 

PRINT USING "K" ; "TEMPERATURE= 

PRINT USING "K" ; "ATMOSPHERIC PRESSURE= 

PRINT USING "K"; "WATER VAPOR PRESSURE = 

PRINT USING "K" j "PULSE LENGTH= 

PRINT USING "K"; "ANTENNA DI AMETER= 

PRINT USING "K" ; "ECHOSOUNDER FREOUENCY= 

PRINT USING "K" ; "POWER TRANSMITTED^ 

PRINT "TEMPERATURE STRUCTURE PROFILE USED " .Profile 
IF Prof i 1 e=4 THEN 

PRINT "INVERSION HEIGHT ", Inver 
END IF 

SELECT Profile 
CASE 1 



" ,Tenp,"CELCIUS" 

" , Atom_pres , "mb" 

" , Wat_pres , "mb" 

" ,Pul se , "ms" 

" , Ant_diam , " m. " 

" ,Freq,"Hz" 

" ,Pou_trans," WATTS" 
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PRINT " 1 

PRINT M 
PRINT M 
CASE 2 

PRINT " 2 

PRINT M 
PRINT M 
PRINT " 

PRINT M 
CASE 3 

PRINT M 3 

PRINT ” 

PRINT " 

PRINT " 

PRINT M 
CASE 4 

PRINT " 4 

PRINT M 
END SELECT 
! 

PRINT USING "K'VTHE MINIMUM DETECTABLE SIGNAL UAS SET AT" . Noise , "WATTS" 
! 

PRINT M M 

j 

j 

PRINT "HIT CONTINUE TO CONTINUE 11 
PAUSE 

OUTPUT KBD; M K" ; 

! 

GRAPHICS ON 
VIEWPORT 15, 120, 10,70 
! 

! 

! PLOT RANGE VERSUS FREQUENCY 

Xnin=0 

Xnax = . 5 

Ynin=PROUND< Ynin-50,2 ) 

Ynax=PROUND( Ynax+50 , 2 ) 

Dx s . 05 
Dy- 100 

Xrange=Xnax-Xnin 
Yr ange=Ynax- Ynin 
WINDOW Xnin , Xnax , Ynin , Ynax 
AXES Dx , Dy f Xnin , Ymin , 2 , 1 
CLIP OFF 
! 

! LABEL PLOT 
CSIZE 4 , . 6 
LDIR 0 
LORG 5 



A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 
AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 " 
WHICH GIVES A HEIGHT TO THE -1.16 PROFILE" 

A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 
AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 " 
WHICH GIVES A HEIGHT TO THE -1.16 PROFILE" 

BUT WITH A FACTOR OF TWO TO APPROXIMATE" 
LOOKING UP A THERMAL PLOOM" 

A TEMPERATURE STRUCTURE PROFILE BASED ON DATA" 
AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398" 
FOR NIGHT CONDITIONS. THE DEPENDENCE WITH" 
HEIGHT IS EXP< - . 001 * HEIGHT ABOVE G5 METERS" 

AND HEIGHT TO THE -1.46 BELOW G5 METERS" 

A TEMPERATURE STRUCTURE PROFILE BASED ON A" 
HEIGHT TO THE -4/3 " 
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MOVE .5*Xrange,Ymax* 1 .05 

LABEL "RANGE VERSUS TRANSDUCER EFFICIENCY" 

! LABEL HORIZONTAL AXES 
LDIR 0 
LORG 5 

FOR J-0 TO Xnax STEP Dx 
CSIZE 4 , . G 

MOVE J, Ynin- . 05* Yrange 
LABEL J 
NEXT J 

MOVE . 5*Xrange , Ynin- . 1 *Yrange 

CSIZE 4 , . G 

LABEL "EFFICIENCY" 

! LABEL VERTICAL AXES 
LORG 8 

FOR J=Ynin TO Ynax STEP Dy 
CSIZE 4 , . G 

MOVE Xnin- .00Z5*Xrange , J 
LABEL USING "K";J 
NEXT J 
LDIR PI/Z 
LORG 6 

MOVE Xnin- . 1 5*Xrange , . 5* Yrange+Ymin 
CSIZE 4, .6 

LABEL "RANGE (METERS)" 

CLIP ON 
I nd= 1 

FOR Ef f = . 05 TO .5 STEP .05 
PLOT Eff.Randnd) 

I nd=I nd+ 1 
NEXT Eff 

l 

PRINT "HIT CONTINUE TO CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT KBD; "K" ; 

! 

I 

LINPUT "WOULD YOU LIKE TO MAKE ANOTHER RUN (ENTER Y OR N)?",X$ 
Again=FNYes( X$) 

Mess^up- 1 
WHILE Mess_up= 1 
SELECT Again 
CASE 1 

LINPUT "WOULD YOU LIKE TO CHANGE A VARIABLE (ENTER Y OR N)?",X$ 
Neu_va=FNYes( X$) 

Ano_change= 1 
WHILE Ano_change=1 
SELECT Neu_va 
CASE 1 



VARIABLE 



CURRENT VALUE* 



USING "K"; 


"1 


TEMPERATURE 


", Temp, "CELSIUS" 


USING "K"; 


"Z 


ATOMOSPHERIC PRESSURE " ,Aton jDres/'nb" 


USING "K"; 

M M 


"3 


WATER VAPOR PRESSURE " , Wat_pres . "mb" 


USING "K"; 


"4 


ANTENNA DIAMETER 


" , Ant_dian , " m. " 


USING "K"; 


"5 


PULSE LENGTH 


".Pulse," ms" 


USING "K"; 


"G 


POWER TRANSMITTED " 


,Pou_trans," WATTS" 


USING M K"; 


"7 


ATMOSPHERIC PROFILE 


" , Prof lie 



PRINT 



PRINT 



PRINT M M 

PRINT "ENTER THE NUMBER OF THE VARIABLE YOU WISH TO" 

PRINT "CHANGE" 

INPUT Var 
SELECT Var 

CASE 1 

INPUT "TEMPERATURE=" , Temp 
CASE Z 

INPUT "ATMOSPHERIC PRESSURE IN mb=" , Atom_pres 
CASE 3 

INPUT "WATER VAPOR PRESSURE IN mb=" , Wat_pres 
CASE 4 

INPUT "ANTENNA DIAMETER IN m . =" , Ant_diam 
CASE 5 

INPUT "PULSE LENGTH IN ns*", Pulse 



CASE S 
INPUT 
CASE 7 
PRINT " 
PRINT * 
PRINT * 
PRINT 1 
PRINT ' 
PRINT " 
PRINT ‘ 
PRINT ‘ 



‘POWER TRANSMITTED IN WATTS=" , Pou_trans 
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— > Ct A Z PROFILE OF Z A <-1.1G>" 

FROM WALTERS/KUNDEL 1981" 

Z ==> SAME AS ONE BUT WITH FACTOR" 

OF TWO TO APPROXIMATE LOOKING- 
UP A THERMAL PLOOM" 

3 ==> Ct A Z PROFILE OF EXP< - . 001 *Z >" 
FROM WALTERS/KUNDEL 1981." 

4 ==> CT A Z PROFILE OF Z A <-4/3)" 

INPUT "ENTER NUMBER OF DESIRED PROFILE" , Prof i 1 e 
IF Prof i 1 e=4 THEN 

INPUT "HEIGHT OF INVERSION IN METERS**" , I nver 
END IF 
CASE ELSE 

PRINT Var, "IS NOT ONE OF THE OPTIONS" 

END SELECT 

LINPUT "MADE ANOTHER CHANGE( ENTER Y OR N)?",X$ 
Neu_va=FNYes(X$) 

Mess_up=Z 
CASE Z 

Mes5_up=Z 
Ano_change c= Z 
CASE ELSE 

PRINT "YOUR CHOICES WHERE Y OR N MM!" 



Ano_change= 1 

LINPUT "CHANGE A VARIABLE (ENTER Y OR N)?",X$ 
Neu^va^FNYesC X$) 

END SELECT 
END WHILE 
CASE 2 

Mess_up=2 

CASE ELSE 

PRINT "YOUR CHOICES UHERE Y OR N !!!!!!" 

Mess_up= 1 

LINPUT "MAKE ANOTHER RUN (ENTER Y OR N)?",X$ 

Again=FNYes( X$ ) 

END SELECT 
END WHILE 
END WHILE 

PRINT "THAT’S ALL, FOLKS" 
end 



! CALCULTE THE ATTENUATION 

sub Attenuat ion( Atom_pres , At ten ,Freq , Temp , Wat_pres ) 

! THI S SUBPROGRAM CALCULATES THE ATTENUATION OF SOUND 
! IN AIR BASED UPON EQUATIONS IN NEFF 1975 



INPUT ATMOSPHERIC PRESSURE IN MILLIBARS 
FREQUENCY OF SOUND WAVE IN HERTZ 
TEMPERATURE IN DEGREES CELCIUS 
WATER-VAPOR PRESSURE IN MILLIBARS 

OUTPUT ATTENUATION IN 1 /METERS 
VARIABLES 

Atom_pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 

Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 

SUBPROGRAM ATTENUATION. 

At-t_max VARIABLE IN SUBPROGRAM ATTENUATION. IT IS THE 
ATTENUATION AT THE FREQUENCY OF THE MAXIMUM 
ATTENUATION FOR THE INPUT CONDITIONS. 

F VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 

RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 

Fmax FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 

ATTENUATION. 

Freq INPUT FREQUENCY OF ECHOSOUNDER. 

H VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 

Ps-tar VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 

CALCULATION. 

Temp INPUT TEMPERATURE IN DEGREES CELSIUS. 

Tstar INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 

IN SUBPROGRAM ATTENUATION. 

Wat_pres ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 



! OPERATOR. 

H= 1 00*Uat_pres/Atom_pres 
Tstar*( 1 .8*Temp+49Z)/519 
Pstar-Atom_pres/ 1014 

Fmax*( 10+GB00*H+44400*H*H)*Pstar/Tstar" .8 

A*tt_max= . 0078*Fmax*Tstar* ( -2 . 5)*EXP( 7.77*( 1-1/Tstar ) ) 

F=Freq/Fmax 

Attend Att_max/ 304.8 >*(< . 1 8*F ) "2 + ( Z*F*F/( 14F*F> >"Zr .5 

Attend Atten+1 . 74E- 1 0*Freq*Freq )/4 . 35 

SUBEND 

j 

def FNYes( X$ ) 

{THIS FUNCTION INTERPRETS THE OPERATORS RESPONSES TO YES NO QUESTIONS 
! 

! INPUT X$ 

! 

! OUTPUT FNYes 

! VARIABLES 

! Temp$ VARIABLE STRING USED IN FUNCTION YES. 

! X$ STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 

! RESPONSE TO YES OR NO QUESTION. 

DIM Temp$[ 1 ] 

TempSt 1 ,1]=TRIM$(X$) 

SELECT TempS 
CASE " Y" , “y" 

RETURN 1 
CASE "N" , "n" 

RETURN 2 
CASE M " 

RETURN 1 
CASE ELSE 

RETURN -2 
END SELECT 

FNEND 

| 
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APPENDIX D 
COMPUTER PROGRAM 4 

! FULLER, ROBERT PR0G_4 ANTENNA SIZE 

! 

! 10 SEP 85 
! 

! 

| *******puRPOSE******* 

! 

! THIS PROGRAM WILL TAKE PARAMETERS OF AN ECHOSOUNDER AND 

! ESTIMATE THE RANGE AS A FUNCTION OF ANTENNA SIZE OF THE THE ECHOSOUNDER. 

! THE FOLLOWING INPUTS ARE REQUIRED: 

! ATMOSPHERIC DATA 

! 1 )ATOMOSPHERIC PRESSURE in millibars 

! 2 )SELECT ONE OF FOUR ATMOSPHERIC PROFILES FOR THE TEMPERATURE 

! STRUCTURE PROFILE. 

! a)FOR ONE PROFILE THE HEIGHT OF THE INVERSION LAYER IS ALSO INPUT 

! 3 )TEMPCRATURE IN DEGREES CELSIUS 

! 4) WATER VAPOR PRESSURE IN millibars 

! ECHOSOUNDER DATA 

I 5)FREQUENCY 

! G)POWER TRANSMITTED BY ECHOSOUNDER IN WATTS 

! 7 )PULSE LENGTH OF THE TRANSMITTED ACOUSTIC ENERGY 

! THE PROGRAM OUTPUTS THE FOLLOWING GRAPH TO AN EXTERNAL PLOTTER: 

! 1 )RANGE( m. ) VERSUS ANTENNA DIAMETER FOR VARIOUS FREQUENCIES 



***** VARIABLES***** 



Again TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 

WISHES TO MAKE ANOUTHER RUN. 

Ano_change TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 

WISHES TO CHANGE THE VALUE OF INPUT DATA BEFORE MAKING 
ANOUTHER RUN. 

Ant_area ECHOSOUNDER ANTENNA AREA CALCULATED FROM INPUT OF 
ANTENNA DIAMETER. 

Ant_diam INPUT OF ECHOSOUNDER ANTENNA DIAMETER IN METERS. 

Atom_pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 

Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 

SUBPROGRAM ATTENUATION. 

Bn BANDWIDTH OF RECEIVER NEEDED TO RECEIVE SIGNALS 

DOPPLER SHIFFED BY THREE METER PER SECOND 
VERTICAL VELOCITIES. 

C VELOCITY OF SOUND CALCULATION AT ZERO CELSIUS. 

Cl VARIABLE USED TO SCALE INVERSION HEIGHT IN 

TEMPERATURE STRUCTURE PROFILE THREE. 

C3 VARIABLE USED IN CALCULATING THE FOURTH CteZ PROFILE. 

CneZ < * ) ARRAY OF VALUES OF THE ACOUSTIC REFRACTIVE INDEX 
PARAMETER. CALCULATED BASED ON SELECTION OF 



PROFILE FOR TEMPERATURE STRUCTURE PARAMETER 
AND ASSUMED PROFILE FOR THE VELOCITY STRUCTURE 
PARAMETER. 

CteZ VALUE OF THE TEMPERATURE STRUCTURE PARAMETER. 

THE OPERATOR SELECTS A PROFILE FROM SEVERAL PROVIDED. 

CveZ VALUE OF THE VELOCITY STRUCTURE PARAMETER. 

VALUES BASED ON CALCULATION USING ASSUMED DISSAPATION 
RATE. 

Dx STEP SIZE FOR X AXIS FOR VARIOUS PLOTS. 

Dy STEP SIZE FOR Y AXIS FOR VARIOUS PLOTS. 

Epsilon DISSAPATION RATE USED IN CALCULATION OF CveZ. 

Ef f TRANSMISSION EFFICIENCY OF ECHOSOUNDER. 

Eff EFFICIENCY OF ECHOSOUNDER WHEN RECEIVING ACOUSTIC 

BACKSCATTER. 

Exc_att EXCESS "ATTENUATION" AT GIVEN RANGE. 

F VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 

RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 

Fmax FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 

ATTENUATION. 

Freq INPUT FREQUENCY OF ECHOSOUNDER. 

G ANTENNA EFFECTIVE APERATURE FACTOR. 

H VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 

I MAIN INDEX FOR LOOP TO CHECK RETURNED POWER UNTIL 

POWER WAS LESS THAN MINIMUM DETECTABLE SICNAL. 

IND INDEX USED FOR LOOP FOR DIFFERENT EFFICIENCIES 

Interval MAXIMUM EFFECTIVE-SCATTERING VOLUME THICKNESS. 

Inver HEIGHT OF INVERSION LAYER. 

J FIRST ORDER INDEX FOR ASSORTED LOOPS. 

K WAVENUMBER 

L SECOND ORDER INDEX FOR ASSORTED LOOPS. 

Mess_up TEST VALUE TO CHECK RESPONSES WHEN OPERATOR IS INPUTTING 
RESPONSES. 

N TEST VALUE IN CALCULATION OF COHERENCE LENTH IN CALCULATION 

OF EXCESS ATTENUATION. 

Neu_va TEST VALUE IN CHECKING IF OPERATOR WISHES TO CHANGE 
A VARIABLE BEFORE A NEW RUN. 

Noise ASSUMED MINIMUM DETECTABLE SIGNAL. 

Pou_back ( * ) POWER BACKSCATTERED FROM GIVEN RANGE. 

Pou_ret< * ) POWER BACKSCATTERED TO ECHOSOUNDER FROM GIVEN RANGE. 

Pou_trans INPUT OF POWER SUPPLIED TO TRANSDUCER OF ECHOSOUNDER. 

Profile OPERATOR INPUT OF CteZ PROFILE FROM AVAILABLE PROFILES. 

Pstar VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 

CALCULATION. 

Pulse TRANSMITTED PULSE LENGTH IN MILLISECONDS. INPUT BY 

OPERATOR. 

R RANGE FOR INPUT TO INTEGRAL OF EXCESS ATTENUATION. 

R1 VARIABLE USED IN SCALING INVERSION HEIGHT FOR PROFILE THREE. 

Range( * ) ARRAY OF RANGE VALUES. 

Remainder REMAINDER OF MODULO FUNCTION USED TO DECREES THE 

NUMBER OF PASSES THROUGH THE INTEGRATION FOR EXCESS 



ATTENUATION. 

Rge VALUE Or RANGE IN INTEGRAL FOR EXCESS ATTENUATION. 

Rho CORRELATION LENGTH USED IN CALCULATION OF EXCESS 

ATTENUATION. 

Signa FRACTION OF ENERGY BACKSCATTERED FROM GIVEN RANGE. 

Speed_sound SPEED OF SOUND IN AIR AT INPUT TEMPERATURE. 

Sunpou_back 

T INPUT TEMPERATURE IN DEGREES KELVIN. 

Temp INPUT TEMPERATURE IN DEGREES CELSIUS. 

Tenp$ VARIABLE STRING USED IN FUNCTION YES. 

Tstar' INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 

IN SUBPROGRAM ATTENUATION. 

Var USED TO SELECT THE VARIABLE THE OPERATOR UISHED TO 

CHANGE BEFORE MAKING ANOTHER RUN. 

Wat_pres ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 
OPERATOR. 

X THIRD ORDER INDEX USED IN VARIOUS LOOPS. 

X$ STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 

RESPONSE TO YES OR NO QUESTION. 

Xmax VALUE OF LARGEST VALUE OF X FOR VARIOUS PLOTS. 

Xmin VALUE OF SMALLEST VALUE OF X FOR VARIOUS PLOTS. 

Xrange VALUE OF RANGE OF X VALUES FOR VARIOUS PLOTS. 

YlabelS LABEL ON Y AXIS PASSED TO SUBPROGRAM Semi_log FOR PLOTTING. 
Ynax LARGEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

Ymn SMALLEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

Yrange RANGE OF Y VARIABLES FOR VARIOUS PLOTS. 

Ze INTERMEDIATE VALUE IN CALCULATION OF EXCESS ATTENUATION. 



(DECLARE VARIABLES 
INTEGER I 

DIM Pou_back< 1500) ,Pou_ret< 1500) ,Range< 1500) 

DIM Cne2< 1500) ,Ran<50,6) 

! 

PLOTTER IS 705 , " HPGL" 

LINE TYPE 1 

i 

(INPUT ATMOSPHERIC DATA 

INPUT "ENTER TEMPERATURE IN DEGREES CELCIUS" , Tenp 
INPUT "ENTER ATMOSPHERIC PRESSURE IN MILLIBARS" ,Atom_pres 
INPUT "ENTER WATER VAPOR PRESSURE IN MILLIBARS" ,Wat_pres 



! 

(INPUT ECHOSOUNDER DATA 

INPUT "ENTER ECHOSOUNDER PULSE LENGTH IN MILLI SECONDS^' , Pulse 
INPUT "ENTER POWER TO TRANSMITTER IN WATTS" ,Pou_1;rans 
G= . 40 (ANTENNA EFFECTIVE APERATURE FACTOR 

Ef f = . 25 



! 

! SELECT THE PROFILE FOR THE TEMPERRTURE STRUCTURE 
PRINT "YOU HRVE R CHOICE OF TEMPERRTURE STRUCTURE PROFILES'* 

PRINT "PROFILE TO BE USED." 

PRINT "YOUR SELECTIONS RRE" 

PRINT M 
PRINT " 

PRINT " 

PRINT " " 

PRINT " 

PRINT " 

PRINT " 

PRINT " " 

PRINT " 

PRINT " 

PRINT " 

PRINT " 

PRINT " 

PRINT " " 

PRINT " 

PRINT " 

PRINT " " 

PRINT' " " 

Mess_up*= 1 
UHILE Mess 

INPUT "ENTER THE DESIRED PROFILE (1 OR Z OR 3 OR 4)", Prof lie 

IF Prof i 1 e=4 THEN 

INPUT "ENTER THE HEIGHT OF THE INVERSION IN METERS" , I nver 

END IF 

IF Profiled THEN 
Mess_up=0 

ELSE 

IF Profiled THEN 
Mess_up=0 

ELSE 

IF Profiled THEN 
Mess_up=0 

ELSE 

IF Prof il e=4 THEN 

Mess_up=0 

ELSE 

PRINT Profile," UR S NOT ONE OF THE OPTIONS!!"" 

Mess_up= 1 
END IF 
END IF 
END IF 
END IF 
END UHILE 
OUTPUT KBD; "K " ; 



1 R TEMPERRRURE STRUCTURE PROFILE BRSED ON DRTR RS " 

PRESENTED IN URLTERS/KUNDEL 1981 PRGE 398 UHICH " 
GIVES R HEIGHT TO THE -1.16 PROFILE" 

Z THE TEMPERRTURE STRUCTURE PROFILE RBOVE BUT" 

WITH R FRCTOR OF TWO TO RPPROXIMRTE LOOKING" 

UP R THERMRL PLOOM" 

3 R TEMPERRTURE STRUCTURE PROFILE BRSED ON DRTR RS" 
PRESENTED IN URLTERS/KUNDEL 1981 PRGE 398 FOR" 
NIGHT CONDITIONS. THE DEPENDENCE UITH HEIGHT " 

IS EXP( - . 001 * HEIGHT RBOVE 65 METERS RND HEIGHT" 

TO THE -1.46 BELOU 65 METERS" 

4 R TEMPERRTURE STRUCTURE PROFILE BRSED ON R" 

HEIGHT TO THE -4/3 " 



up* 1 
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Again= 1 
WHILE Again* 1 1 
FOR J»1 TO I 

Pou_back ( J )=0 
Pou_ret( J ) =0 
Range( J )=0 
CneZ( J ) =0 
NEXT J 

I 

! CONVERT TEMPERATURE TO KELVIN 
T=Tenp+Z73 

I 

! CALCULATE SPEED OF SOUND BASED ON INPUT TEMPERATURE( CELSIUS) 
Speed_sound=Z0. 05* ( T ) A . 5 

! CALCULATE THE SPEED OF SOUND AT 0 DEGREES CELCIUS 
C=Z0. 05*Z73 A . S 

1=1 (INDEX FOR INNER LOOP (POWER RETURNED TO NOISE) 

I nd= 1 (INDEX FOR SECOND LOOP (EFFICIENCY) 

I f req= 1 (INDEX FOR THIRD LOOP (FREQUENCY) 

Pou_back ( 0 )=0 (INITIALIZE VARIABLE FOR POWER BACKSCATERRED 
Intervals Speed_sound*Pul se+ 1 . E-3)/Z 

(INITIALIZE PLOTTING PARAMETERS SO CAN COMPARE DURING RUN TO SCALE 

(PLOTS 

Ynin= 1 000 

Ynax=0 

(CALCULATE RETURNED SIGNAL POWER UNTIL IT IS LESS THAN BACKGROUND 
(FOR EACH FREQUENCY. 

Range( 0)=0 

FOR Freq=500 TO 1500 STEP 500 
Sumpou_back=0 

FOR Ant_diam=.5 TO G.Z5-Hreq STEP .Z5 
Ant_area=PI*( Ant_diam/Z ) A Z 

Bn=Z*Freq*( 1 - 1 /( 6/Speed_sound+ 1 ) ) ( BANDWIDTH FOR 3M/S VERTICAL VELOCITY 
Noise= 1 . 38E-Z3*Bn* ( T HZ . E- 1 4 ! MINI MUM DETECTABLE SIGNAL AS 
! JOHNSON NOISE( NEGLIGABLE H ESTIMATED BACKGROUND 

(CALCULATE THE ATTENUATION COEFFICIENT BASED ON ATMOSPHERIC DATA 
CALL Attenuat ion( Atom_pres ,Atten,Freq, Tenp , Wat_pres ) 

K=Z*PI *Freq/Speed_sound ! WAVENUMBER 
1=1 

REPEAT 

Range( I )=Range( I- 1 HZ 
SELECT Profile 
CASE 1 

(THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
IAS PRESENTED IN UALTERS/KUNDEL 1981 PAGE 398. 

(WAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 



! TULAROSA BASIN , NEW MEXICO. 

! AN ADDITONAL FACTOR OF 1 . 8 UAS INCLUDED AS PER DIFFERENT 
! AVERAGING TIME. 

CteZ=Z . 1 Z*Range( I ) A (-1 . 1G) 

CASE Z 

•THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
•AS PRESENTED IN UALTERS/KUNDEL 1981 PAGE 398. 

! UAS FOR MIDDAY AND CLEAR UEATHER ABOVE THE 
•TULAROSA BASIN, NEW MEXICO. 

•AN ADDITONAL FACTOR OF 1 . 8 UAS INCLUDED AS PER DIFFERENT 
j AVERAGING TIME AND AN ADDITIONAL FACTOR OF TUO UAS 
! INCLUDED TO APPROXIMATE LOOKING UP A THERMAL PLOOM 
CteZ=Z*Z . 1 Z*Range( I ) A ( -1 . 1G) 

CASE 3 

! THIS TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS M 
• PRESENTED IN UALTERS/KUNDEL 1981 PAGE 398 FOR’ 1 
! NIGHT CONDITIONS. THE DEPENDENCE UITH HEIGHT " 

IIS EXP( - . 001 * HEIGHT ABOVE G5 METERS AND HEIGHT" 

! TO THE -1.4G BELOU G5 METERS" 

IF Range< I ) <65 THEN 

CteZ=75 . 5* Range ( I ) A ( - 1 .4G) 

' ELSE 

CteZ=3. GGE-Z*EXP( - .001 *Range( I ) ) 

END IF 
CASE 4 

! CALCULATE THE TEMPERATURE STRUCTURE FACTOR BASED ON ASSUMPTION 
! THAT IT’S PROPORTIONAL TO HEIGHT TO THE -4/3 POUER 

! THAT THE SURFACE VERTICAL HEAT FLUX IS .095 C.m./sec 

I THE DEPENDENCE UITH HEIGHT ABOVE THE INVERSION I PULLED FROM 
•MY LEFT EAR 

•EQUATION FROM NEFF, 1975 
SELECT Range( I )/Inver 
CASE <.9 

C3-< ( . 0Z4 )*( T ) A ( . GG7 ) ) 

CteZ = C3*( Range( I ) ) A ( - 1 . 33) 

CASE .9 TO 1 
CteZ=Ct eZ 
C 1 =CteZ 
R1 =Range< I ) 

CASE 1 TO 1 .3 

Ct eZ= 1 0 A ( ( LGT( C3*Range( I ) )-LGT< C3*R1 ) )+LGT< Cl ) ) 
RZ=Range( I ) 

CZ=Ct eZ 
CASE ELSE 

Ct eZ = C3* ( Range ( I ) A ( - 1 . 33 ) -RZ A ( - 1 . 33 ) ) + CZ 
END SELECT 
END SELECT 

! CALCULATE VELOCITY STRUCTURE FACTOR. FORMULA FROM GAYNOR 77 
Epsil on=( . Z866/Range( I ) )*( 1 + . 07* ( Range( I ) ) A ( . G ) ) A ( 1 . 5 ) 

CveZ=Z*Epsi Ion A ( .GG7) 



! CALCULATE ACOUSTIC REFRACTIVE INDEX STRUCTURE FACTOR. 

! FORMULA FROM TETARSKI 

CneZ< I ) = ( CteZ/< Z.98E+5) ) + ( CveZ/< C*C) ) 

! CALCULATE THE FRACTION OF POUER BACKSCATTERED FROM INTERVAL. 

! EQUATION FROM NEFF, 1975 
Signal .00 39*<K A < 1 /3) )*CteZ )/( T) A Z 

I 

(CALCULATE THE EXCESS ATTENUATION. EQUATION FROM CLIFFORD, 1980 
(THE EXCESS ATTENUATION IS Exc_att 
IF I <50 THEN 

Remainder=0 

ELSE 

IF I < 1 50 THEN 

Reroainder=Range< I ) MODULO 10 

ELSE 

Remainder=Range< I ) MODULO Z0 
END IF 
END IF 

IF Remainder^© THEN 
Rho=0 
L=0 
H=0 

Rge=Range( I ) (CONSTANT IN INTEGRAL 
R=0 

FOR J=0 TO Z*I 

F=CneZ ( INT< J/Z+1 ) ) 

F=F* < 1-R/<Rge) A < 1 . 57 ) + ( R/< Rge ) ) A < I .57) ) 

IF J>0 THEN 

IF J<Z*I THEN 

IF I NT ( J/Z )= J/Z THEN 
L=L+F 
F=0 

ELSE 

H=H+F 

F=0 

END IF 
END IF 
END IF 
Rho=Rho+F 
R=R+ 1 

NEXT J 

Rho=Rho+4*L+Z*H 

Rho=( < (Rho* .33>*K*K* 1 . 45 ) A ( - . G ) ) 

N~< Ant_diari/Rho ) A Z 
IF N<= 1 THEN 

Ze=1 /( 1+N) 

ELSE 

Ze- 1 . 5/ < 1+N) 

(STEP OF 1 . 5==>SEE CLIFFORD 1980 
END IF 
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Exc att*Ze»Ze 



ELSE 

Exc_att=Exc_att 
END IF ' 

I 

! CALCULATE THE POUER BACKSCATTERED 

Pou_back ( I ) = < Pou_trans*Ef f -Sumpou_back ) *EXP( - Atten*Range< I ) ) 
Pou_back( I )=Pou_back< I )*Interval*Exc_att*Signa 
Sunpou_back*Sunpou_back+Pou_back ( I > 

! CALCULATE THE POUER BACKSCATTERED TO THE ANTENNA 
Pou_ret< I )=Pou_back ( I ) *EXP< - Atten*Range< I ) ) * Ant_ar ea*G*Ef f / Range ( I ) A Z 
PRINT "RANGE* 3 " , Range( I ) 

PRINT "POUER RETURNED*" ,Pou_ret( I ) 

PRINT "FREQUENCY*", Freq 
1 = 1 + 1 

UNTIL Pou_ret< 1-1 ) >0 AND Pou_ret< I- 1 ) <=Noise 
Ran< Ind, If req )*Range< I-Z ) 

IF Ran( I nd , If req ) >Ynax THEN 
Ymax=Ran< Ind, If req) 

END IF 

IF Ran< Ind , If req > <Ymin THEN 
Ynin=Ran( Ind , If req ) 

END IF 
Ind=Ind+ 1 
Again* 1 
FOR J*0 TO I 

Pou_back < J ) =0 
Pou_ret( J >=0 
Range< J)=0 
CneZ< J>=0 
NEXT J 

NEXT Ant_dian 
If req=If req+ 1 
I nd= 1 
NEXT Freq 

I 
| 

OUTPUT KBD; "K" ; 

PRINT "INPUT CONDITIONS" 

| 

PRINT " " 

PRINT USING "K"; "TEMPERATURE* 

PRINT USING "K";" ATMOSPHERIC PRESSURE* 

PRINT USING "K" ; "WATER VAPOR PRESSURE* 

PRINT USING "K" ; "PULSE LENGTH* 

PRINT USING "K"; "TRANSDUCER EFFICIENCY* 

PRINT USING "K"; "POUER TRANSMITTED* 

PRINT "TEMPERATURE STRUCTURE PROFILE USED ", Profile 
IF Profile=4 THEN 

PRINT "INVERSION HEIGHT ".Inver 



" ,Tenp,"CELCIUS" 

" , Aton_pr es , "mb" 

" , Uat_pres , "nb" 

" ,Pul se , "ns" 

" ,Ef f 

" ,Pou_trans," UATTS" 
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END IF 

SELECT Profile 
CASE 1 

PRINT " 

PRINT M 
PRINT " 

CASE 2 

PRINT " 

PRINT " 

PRINT M 
PRINT M 
PRINT M 
CASE 3 

PRINT " 

PRINT " 

PRINT " 

PRINT M 
PRINT " 

CASE 4 

PRINT M 
PRINT M 
END SELECT 
! 

PRINT USING "K" ; "THE MINIMUM DETECTABLE SIGNAL UAS SET AT" , Noise , "WATTS" 
! 

PRINT M M 
! 

PRINT "HIT CONTINUE TO CONTINUE" 

PAUSE 

OUTPUT KBD? "K" ; 

! 

GRAPHICS ON 
VIEWPORT 15,120,10,70 
! 
l 

! PLOT RANGE VERSUS ANTENNA DIAMETER FOR VARIOUS FREQUENCIES 

Xroin=0 

Xmax=G 

Ymin=PROUND( Ymin-50,2 ) 

Ymax=PROUND( Ynax+50 , 2 ) 

Dx = 1 
Dy= 1 00 

Xrange=Xnax-Xmin 
Yrange=Ynax- Ymin 
WINDOW Xnin , Xnax , Ynin , Ymax 
AXES Dx ,Dy, Xnin, Ynin, 1 , 1 
CLIP OFF 
! 

! LABEL PLOT 



1 A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 
AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 M 
WHICH GIVES A HEIGHT TO THE -1.16 PROFILE" 

2 A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 
AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 " 
WHICH GIVES A HEIGHT TO THE -1.16 PROFILE" 

BUT WITH A FACTOR OF TWO TO APPROXIMATE" 
LOOKING UP A THERMAL PLOOM" 

3 A TEMPERATURE STRUCTURE PROFILE BASED ON DATA" 
AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398" 
FOR NIGHT CONDITIONS. THE DEPENDENCE WITH" 
HEIGHT IS EXP( - . 001 * HEIGHT ABOVE G5 METERS" 

AND HEIGHT TO THE -1.46 BELOW G5 METERS" 

4 A TEMPERATURE STRUCTURE PROFILE BASED ON A" 

HEIGHT TO THE -4/3 " 
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CSIZE 4, .6 
LDIR 0 
LORG 5 

MOVE . 5*Xrange , Ymax* 1 . 0S 

LABEL "RANGE VERSUS ANTENNA DIAMETER" 

! LABEL HORIZONTAL AXES 
LDIR 0 
LORG 5 

FOR J=0 TO Xmax STEP Dx 
CSIZE 4, .6 

MOVE J , Ymin- . 05* Yrange 
LABEL J 
NEXT J 

MOVE . 5*Xrange , Ymin- . 1 * Yrange 
CSIZE 4, .6 

LABEL "ANTENNA DIAMETER <m.)" 

! LABEL VERTICAL AXES 
LORG 8 

FOR J=Ymin TO Ynax STEP Dy 
CSIZE 4 , . G 

MOVE Xmin- . 00Z5*Xrange , J 
LABEL USING "K" ;J 
NEXT J 
LDIR PI/Z 
LORG G 

MOVE Xnin- - 1 S*Xrange , . S* Yrange + Ymin 
CSIZE 4 , . 6 

LABEL "RANGE (METERS)" 

CLIP ON 
I nd= 1 
If req 3 1 

FOR Freq=500 TO 1 G00 STEP B00 

FOR Ant^diam- . 5 TO 6.Z5-Ifreq STEP .Z5 
IF INT( Ifreq/Z )=Ifreq/Z THEN 
LINE TYPE 3 
ELSE 

LINE TYPE 1 
END IF 

PLOT Ant_diam , Ran< Ind , I f req ) 

I nd=I nd+ 1 
NEXT Airt_diam 
LDIR 0 
LORG Z 
CSIZE 3 , . G 

MOVE Ant_dian- . 1 , Ran( Ind- 1 , If req ) 

LINE TYPE I 

LABEL USING "K" ; Freq , "Hz " 

If req=If req+ 1 
Ind= 1 
NEXT Freq 
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! 

PRINT "HIT CONTINUE TO CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT KBD; "K" ; 

! 

! 

LINPUT "WOULD YOU LIKE TO MAKE ANOTHER RUN (ENTER Y OR N)?",X$ 
Again=FNYes( X$ ) 

Mess_up= 1 
WHILE Mess_up= 1 
SELECT Again 
CASE 1 

LINPUT "WOULD YOU LIKE TO CHANGE A VARIABLE (ENTER Y OR N>?“ ,X$ 
Neu_va=FNYes( X$ ) 

Ano_change= 1 
WHI LE Ano_change= I 
SELECT Neu_va 
CASE 1 

PRINT M VARIABLE CURRENT VALUE" 



PRINT 


USING 


"K" ; ' 


"1 


TEMPERATURE 


Tenp , "CELSIUS" 


PRINT 


USING 


"K'V 


•z 


ATOMOSPHERIC PRESSURE 


" , Atom _pres , " mb' 


PRINT 

PRINT 


USING 

«i ii 






WATER VAPOR PRESSURE 


" , Wat_pres , "mb" 


PRINT 


USING 




*4 


ANTENNA DIAMETER 


" , Ant_diam , " m, " 


PRINT 


USING 


"K'V 


"5 


PULSE LENGTH 


", Pulse," ms" 


PRINT 


USING 


"K'V 


"B 


POUER TRANSMITTED " .Poujrans," WATTS' 


PRINT 


USING 


"K" 5 


"7 


ATMOSPHERIC PROFILE 


" , Prof lie 



PRINT " " 

PRINT "ENTER THE NUMBER OF THE VARIABLE YOU WISH TO" 

PRINT "CHANGE" 

INPUT Var 
SELECT Var 

CASE 1 

INPUT "TEMPERATURE*" ,Tenp 
CASE Z 

INPUT "ATMOSPHERIC PRESSURE IN mb*" , Atom_pres 
CASE 3 

INPUT "WATER VAPOR PRESSURE IN mb=" , Wat_pres 
CASE 4 

INPUT "ANTENNA DIAMETER IN m. =" , Ant_diam 
CASE B 

INPUT "PULSE LENGTH IN ms=", Pulse 
CASE B 

INPUT "POWER TRANSMITTED IN WATTS=" , Pou_trans 
CASE 7 

PRINT " 1 ==> Ct A Z PROFILE OF Z A (-1.1B)" 

PRINT " FROM WALTERS/KUNDEL 1981" 

PRINT " Z ==> SAME AS ONE BUT WITH FACTOR" 
PRINT " OF TWO TO APPROXIMATE LOOKING" 
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PRINT “ UP A THERMAL PLOOM** 

PRINT M 3 = «> C-TZ PROFILE OF EXP( - . <d ®\ * Z ) M 

PRINT M FROM WALTERS/KUNDEL 1981.” 

PRINT “ 4 =“> CT* 2 PROFILE OF Z*(-4/3>“ 

INPUT “ENTER NUMBER OF DESIRED PROFILE” .Prof i 1 e 
IF Prof i 1 e=4 THEN 

INPUT “HEIGHT OF INVERSION IN METERS-** , I nver 
END IF 
CASE ELSE 

PRINT Var .“IS NOT ONE OF THE OPTIONS*' 

END SELECT 

LINPUT “MADE ANOTHER CHANGE( ENTER Y OR N)?“,X$ 
Neu_va=FNYes( X$) 

Mess_up= a Z 
CASE Z 

Mess_up e 2 
Ano^change-Z 
CASE ELSE 

PRINT “YOUR CHOICES UHERE Y OR N !!!!!“ 

Ano_change*= 1 

LINPUT “CHANGE A VARIABLE (ENTER Y OR N>?“ ,XS 
Neu_va=FNYes( X$ ) 

END SELECT 
END WHILE 
CASE 2 

Mess_up=2 
CASE ELSE 

PRINT “YOUR CHOICES UHERE Y OR N !!!!!!“ 

Mess_up= 1 

LINPUT “MAKE ANOTHER RUN (ENTER Y OR N)?“ ,X$ 

Again=FNYes( X$ ) 

END SELECT 
END WHILE 
END WHILE 

PRINT “THAT’S ALL. FOLKS" 
end 
! 

! 

ICALCULTE THE ATTENUATION 

sub Attenuat ion( At om_pres , At ten , Freq , Temp , Uat_pres ) 

! THI S SUBPROGRAM CALCULATES THE ATTENUATION OF SOUND 
! IN AIR BASED UPON EQUATIONS IN NEFF 1 97S 
I 

! INPUT ATMOSPHERIC PRESSURE IN MILLIBARS 
! FREQUENCY OF SOUND WAVE IN HERTZ 

! TEMPERATURE IN DEGREES CELCIUS 

! WATER- VAPOR PRESSURE IN MILLIBARS 

! 

! OUTPUT ATTENUATION IN 1 /METERS 
! 
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! VARI ABLES 
! Atom_pres 
! Atten 

! Att max 



F 

Fmax 

Freq 

H 

Pstar 

Temp 
Tst ar 

Uat_pres 



INPUT OF ATMOSPHERIC PRESSURE IN mb. 

ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 
SUBPROGRAM ATTENUATION. 

VARIABLE IN SUBPROGRAM ATTENUATION. IT IS THE 
ATTENUATION AT THE FREQUENCY OF THE MAXIMUM 
ATTENUATION FOR THE INPUT CONDITIONS. 

VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 

FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
ATTENUATION. 

INPUT FREQUENCY OF ECHOSOUNDER. 

VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 

VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 

INPUT TEMPERATURE IN DEGREES CELSIUS. 

INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 
IN SUBPROGRAM ATTENUATION. 

ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 
OPERATOR. 



H= I 00* Uat_pres/ At om_pres 
Tstar=( 1 . 8*Temp+49Z )/S1 9 
Pst ar= At om_pr es/ 1014 

Fmax=( 1 0+6600*H444400*H*H )*Pstar/Tstar" .8 



Att_max= . 0078* Fmax* Tst ar'' (-Z.5)*EXP(7.77*( 1-1/Tstar ) ) 
F=Freq/Fmax 

Atten=( Att_max/304. 8)*( ( . 1 8* F ) “ 2 + ( Z * F* F/ ( 1 +F * F ) ) A 2 ) A . 5 

Att en=( Att en+ 1 . 74E- 10*Freq*Freq>/4. 35 

SUBEND 



def FNYes( )($ > 

! THIS FUNCTION INTERPRETS THE OPERATORS RESPONSES TO YES NO QUESTIONS 
! 

! INPUT X$ 

I 

! OUTPUT FNYes 

! VARIABLES 

! Temp$ VARIABLE STRING USED IN FUNCTION YES. 

! X$ STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 

! RESPONSE TO YES OR NO QUESTION. 

DIM Temp$l 1 1 
T emp$[ 1 , 1 ] =TRI M$( X$ ) 

SELECT Temp$ 

CASE "rVy" 

RETURN 1 
CASE "N" , “ n" 

RETURN 2 
CASE M M 

RETURN 1 
CASE ELSE 
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RETURN -2 
END SELECT 

FNEND 

I 
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APPENDIX E 
COMPUTER PROGRAM 5 



! FULLER, ROBERT 
| 

! 10 SEP 85 
! 



PROG 5 



♦♦♦♦♦♦♦PURPOSE******* 



POUER 



! THIS PROGRAM WILL TAKE PARAMETERS OF AN ECHOSOUNDER AND 

! ESTIMATE THE RANGE AS A FUNCTION OF POUER TRANSMITTED BY THE ECHOSOUNDER 
! FOR A RANGE OF FREQUENCIES 
! THE FOLLOWING INPUTS ARE REQUIRED: 

! ATMOSPHERIC DATA 

! 1 ) ATOMOSPHERIC PRESSURE in millibars 

! 2 ) SELECT ONE OF FOUR ATMOSPHERIC PROFILES FOR THE TEMPERATURE 

! STRUCTURE PROFILE. 

! a)FOR ONE PROFILE THE HEIGHT OF THE INVERSION LAYER IS ALSO INPUT 

! 3) TEMPERATURE IN DEGREES CELSIUS 

! 4)UATER VAPOR PRESSURE IN millibars 

! ECHOSOUNDER DATA 
! 5) FREQUENCY 

! G) ANTENNA DIAMETER IN METERS 

! 7 )PULSE LENGTH OF THE TRANSMITTED ACOUSTIC ENERGY 

! THE PROGRAM OUTPUTS THE FOLLOWING GRAPH TO AN EXTERNAL PLOTTER: 

! 1 )RANGE< m. ) VERSUS POUER TRANSMITTED IN WATTS FOR VARIOUS FREQUENCIES 

! 

! *****VARI ABLES***** 



I 

l 



Again TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 

WISHES TO MAKE ANOUTHER RUN. 

Ano_change TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 

WISHES TO CHANGE THE VALUE OF INPUT DATA BEFORE MAKING 
ANOUTHER RUN. 

Ant_ar ea ECHOSOUNDER ANTENNA AREA CALCULATED FROM INPUT OF 

ANTENNA DIAMETER. 

Ant_diam INPUT OF ECHOSOUNDER ANTENNA DIAMETER IN METERS. 

Atom_pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 

Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 

SUBPROGRAM ATTENUATION. 

Bn BANDWIDTH OF RECEIVER NEEDED TO RECEIVE SIGNALS 

DOPPLER SHIFFED BY THREE METER PER SECOND 
VERTICAL VELOCITIES. 

C VELOCITY OF SOUND CALCULATION AT ZERO CELSIUS. 

Cl VARIABLE USED TO SCALE INVERSION HEIGHT IN 

TEMPERATURE STRUCTURE PROFILE THREE. 

C3 VARIABLE USED IN CALCULATING THE FOURTH Cte2 PROFILE. 

Cne2 ( * ) ARRAY OF VALUES OF THE ACOUSTIC REFRACTIVE INDEX 
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CteZ 


PARAMETER. CALCULATED BASED ON SELECTION OF 
PROFILE FOR TEMPERATURE STRUCTURE PARAMETER 
AND ASSUMED PROFILE FOR THE VELOCITY STRUCTURE 
PARAMETER. 

VALUE OF THE TEMPERATURE STRUCTURE PARAMETER. 

THE OPERATOR SELECTS A PROFILE FROM SEVERAL PROVIDED. 


C veZ 


VALUE OF THE VELOCITY STRUCTURE PARAMETER. 

VALUES BASED ON CALCULATION USING ASSUMED DISSAPATION 
RATE. 


Dx 

Dy 

Epsilon 

Eff 

Eff 


STEP SIZE FOR X AXIS FOR VARIOUS PLOTS. 

STEP SIZE FOR Y AXIS FOR VARIOUS PLOTS. 
DISSAPATION RATE USED IN CALCULATION OF CveZ. 
TRANSMISSION EFFICIENCY OF ECHOSOUNDER. 

EFFICIENCY OF ECHOSOUNDER WHEN RECEIVING ACOUSTIC 
BACKSCATTER. 


Exc_att 

F 


EXCESS M ATTENUATI ON 1 ' AT GIVEN RANGE. 

VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 


Fnax 


FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
ATTENUATION. 


Freq 

G 

H 

I 


INPUT FREQUENCY OF ECHOSOUNDER. 

ANTENNA EFFECTIVE APERATURE FACTOR. 

VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 
MAIN INDEX FOR LOOP TO CHECK RETURNED POWER UNTIL 
POWER WAS LESS THAN MINIMUM DETECTABLE SIGNAL. 


If req 
IND 

I nt erval 
I nver 
J 
K 
L 

Mess_up 


INDEX FOR INCREMENTS OF FREQUENCY 

INDEX USED FOR LOOP FOR DIFFERENT TRANSMITTED POWER 
MAXIMUM EFFECTIVE-SCATTERING VOLUME THICKNESS. 

HEIGHT OF INVERSION LAYER. 

FIRST ORDER INDEX FOR ASSORTED LOOPS. 

WAVENUMBER 

SECOND ORDER INDEX FOR ASSORTED LOOPS. 

TEST VALUE TO CHECK RESPONSES WHEN OPERATOR IS INPUTTING 
RESPONSES. 


N 


TEST VALUE IN CALCULATION OF COHERENCE LENTH IN CALCULATION 
OF EXCESS ATTENUATION. 


Neu_va 


TEST VALUE IN CHECKING IF OPERATOR WISHES TO CHANGE 
A VARIABLE BEFORE A NEW RUN. 


Noise 


ASSUMED MINIMUM DETECTABLE SIGNAL. 



Pou_back<*> POWER BACKSCATTERED FROM GIVEN RANGE. 

Pou ret( * ) POWER BACKSCATTERED TO ECHOSOUNDER FROM GIVEN RANGE. 



Pou_trans 
Prof i 1 e 
Pstar 


INPUT OF POWER SUPPLIED TO TRANSDUCER OF ECHOSOUNDER. 
OPERATOR INPUT OF CteZ PROFILE FROM AVAILABLE PROFILES. 
VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 


Pul se 


TRANSMITTED PULSE LENGTH IN MILLISECONDS. INPUT BY 
OPERATOR. 


R 

R1 

Range( * ) 


RANGE FOR INPUT TO INTEGRAL OF EXCESS ATTENUATION. 

VARIABLE USED IN SCALING INVERSION HEIGHT FOR PROFILE THREE 
ARRAY OF RANGE VALUES. 
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I 

! 



Remainder 



Rge 

Rho 

Sigma 

Speed_sound 



REMAINDER OF MODULO FUNCTION USED TO DECREES THE 
NUMBER OF PASSES THROUGH THE INTEGRATION FOR EXCESS 
ATTENUATION. 

VALUE OF RANGE IN INTEGRAL FOR EXCESS ATTENUATION. 
CORRELATION LENGTH USED IN CALCULATION OF EXCESS 
ATTENUATION. 

FRACTION OF ENERGY BACKSCATTERED FROM GIVEN RANGE. 
SPEED OF SOUND IN AIR AT INPUT TEMPERATURE. 



Sumpou_back 

T INPUT TEMPERATURE IN DEGREES KELVIN. 

Temp INPUT TEMPERATURE IN DEGREES CELSIUS. 

Temp$ VARIABLE STRING USED IN FUNCTION YES. 



Tstar INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 

IN SUBPROGRAM ATTENUATION. 

Var USED TO SELECT THE VARIABLE THE OPERATOR UISHED TO 

CHANGE BEFORE MAKING ANOTHER RUN. 

Uat_pres ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 
OPERATOR. 

X THIRD ORDER INDEX USED IN VARIOUS LOOPS. 

X$ STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 

RESPONSE TO YES OR NO QUESTION. 

Xmax VALUE OF LARGEST VALUE OF X FOR VARIOUS PLOTS. 

Xmin VALUE OF SMALLEST VALUE OF X FOR VARIOUS PLOTS. 

Xrange VALUE OF RANGE OF X VALUES FOR VARIOUS PLOTS. 

YlabelS LABEL ON Y AXIS PASSED TO SUBPROGRAM SemiJLog FOR PLOTTING. 
Ymax LARGEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

Ymin SMALLEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

Yrange RANGE OF Y VARIABLES FOR VARIOUS PLOTS. 

Ze INTERMEDIATE VALUE IN CALCULATION OF EXCESS ATTENUATION. 



! DECLARE VARIABLES 
INTEGER I 

DI M Pou_back ( 1 500) ,Pou_ret ( 1 500 ) , Range ( 1 500) 

DIM CneZ( 1500) ,Ran< 50,5) 

I 

PLOTTER IS 705, "HPGL" 

LINE TYPE 1 

l 

J INPUT ATMOSPHERIC DATA 

INPUT "ENTER TEMPERATURE IN DEGREES CELCI US" , Temp 
INPUT "ENTER ATMOSPHERIC PRESSURE IN MILLIBARS" , Atom_pres 
INPUT "ENTER WATER VAPOR PRESSURE IN MILLIBARS" , Wat_pres 



! INPUT ECHOSOUNDER DATA 

INPUT "ENTER ECHOSOUNDER PULSE LENGTH IN MILLI SECONDS=" , Pulse 
INPUT "ENTER THE ECHOSOUNDER ANTENNA DIAMETER IN METERS" , Ant_diam 
G= . 40 ! ANTENNA EFFECTIVE APERATURE FACTOR 
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Ef f = . Z5 
! 

! 

! SELECT THE PROFILE FOR THE TEMPERATURE STRUCTURE 
PRINT "YOU HAVE A CHOICE OF TEMPERATURE STRUCTURE PROFILES" 

PRINT "PROFILE TO BE USED." 

PRINT "YOUR SELECTIONS ARE" 

PRINT " 1 

PRINT " 

PRINT " 

PRINT " " 

PRINT " Z 

PRINT " 

PRINT " 

PRINT " " 

PRINT " 3 

PRINT " 

PRINT " 

PRINT " 

PRINT " 

PRINT " " 

PRINT " 4 

PRINT " 

! 

PRINT " " 

PRINT " " 

Mess_up= 1 
UHILE Mess_up=1 

INPUT "ENTER THE DESIRED PROFILE (1 OR Z OR 3 OR 4)" .Profile 

IF Profiled THEN 

INPUT "ENTER THE HEIGHT OF THE INVERSION IN METERS" , Inver 

END IF 

IF Profiled THEN 
Mess_up=0 

ELSE 

IF Profiled THEN 
Mess_up=0 

ELSE 

IF Prof i 1 e=3 THEN 
Mess__up=0 

ELSE 

IF Prof i 1 e=4 THEN 

Mess_up=0 

ELSE 

PRINT Profile," UAS NOT ONE OF THE OPTIONS!!!!" 

Mess_up= 1 
END IF 
END IF 
END IF 
END IF 



A TEMPERARURE STRUCTURE PROFILE BASED ON DATA AS 
PRESENTED IN UALTERS/KUNDEL 1381 PAGE 398 UHICH " 
GIVES A HEIGHT TO THE -1.16 PROFILE" 

THE TEMPERATURE STRUCTURE PROFILE ABOVE BUT" 

WITH A FACTOR OF TUO TO APPROXIMATE LOOKING" 

UP A THERMAL PLOOM" 

A TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS" 
PRESENTED IN UALTERS/KUNDEL 1981 PAGE 398 FOR" 
NIGHT CONDITIONS. THE DEPENDENCE UITH HEIGHT " 

IS EXP( -.001 * HEIGHT ABOVE 65 METERS AND HEIGHT" 

TO THE -1.46 BELOU 65 METERS" 

A TEMPERATURE STRUCTURE PROFILE BASED ON A" 

HEIGHT TO THE -4/3 " 




END UHILE 
OUTPUT KBD; ,, K n ; 

! 

! 

Again* 1 
UHILE Again* I 
FOR M TO I 

Pou_back ( J ) =0 
Pou_re*t( J ) =0 
Range< J ) =0 
CneZ ( J ) *0 
NEXT J 
! 

! CONVERT TEMPERATURE TO KELVIN 

T=Temp+Z73 

l 

! CALCULATE SPEED OF SOUND BASED ON INPUT TEMPERATURE( CELCIUS) 
Speed_sound=Z0. 05* ( T ) A . 5 

! CALCULATE THE SPEED OF SOUND AT 0 DEGREES CELCIUS 
C=Z0.05*Z73 A . 5 
I 

1-1 ! INDEX FOR INNER LOOP (POUER RETURNED TO NOISE) 

Ind-1 ! INDEX FOR SECOND LOOP (EFFICIENCY) 

Ifreq-1 ! INDEX FOR THIRD LOOP (FREQUENCY) 

Pou_back<0)=0 UNITIALIZE VARIABLE FOR POUER BACKSCATERRED 
Interval** Speed_sound*Pul se* 1 .E-3)/Z 
Ant_area=PI *( Ant_dian/Z ) A Z 
! 

(INITIALIZE PLOTTING PARAMETERS SO CAN COMPARE DURING RUN TO SCALE 

(PLOTS 

Yn i n= 1 000 

Ymax=0 

(CALCULATE RETURNED SIGNAL POUER UNTIL IT IS LESS THAN BACKGROUND 
(FOR EACH FREQUENCY. 

Range( 0)=0 

FOR Fr eq=500 TO 1500 STEP 500 
FOR Pou_trans=50 TO 500 STEP 50 

Bn=Z*Freq*( 1 - 1 / < 5/ Speed__sound+ 1 ) ) (BANDUIDTH FOR 3M/S VERTICAL VELOCITY 
Noise=1 . 38E-Z3*Bn*( T)+Z .E- 14! MINIMUM DETECTABLE SIGNAL AS 
! JOHNSON NOISE( NEGLIGABLEH ESTIMATED BACKGROUND 

(CALCULATE THE ATTENUATION COEFFICIENT BASED ON ATMOSPHERIC DATA 
CALL Attenuat ion( Atom_pres , At ten , Freq, Temp , Uat_pres ) 
K=Z*PI*Freq/Speed_sound ! UAVENUMBER 
1=1 

Sunpou_back =0 
REPEAT 

Range( I )=Range( I- 1 ) + Z 
SELECT Profile 
CASE 1 

(THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
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mS PRESENTED IN UALTERS/KUNDEL 1981 PAGE 398. 

!UAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
! TULAROSA BASIN, NEW MEXICO. 

IAN ADDITONAL FACTOR OF 1 . 8 WAS INCLUDED AS PER DIFFERENT 
‘AVERAGING TIME. 

Ct eZ**2 . 1 Z*Range* I ) A * - 1 . 1 6 ) 

CASE Z 

‘THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
I AS PRESENTED IN UALTERS/KUNDEL 1981 PAGE 398. 

! WAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
‘TULAROSA BASIN, NEW MEXICO. 

‘AN ADDITONAL FACTOR OF 1 . 8 WAS INCLUDED AS PER DIFFERENT 
‘AVERAGING TIME AND AN ADDITIONAL FACTOR OF TWO WAS 
‘INCLUDED TO APPROXIMATE LOOKING UP A THERMAL PLOOM 
Cte2“Z*2 . 1 Z*Range* I ) A (-1 . 16) 

CASE 3 

ITHIS TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS" 

‘PRESENTED IN UALTERS/KUNDEL 1981 PAGE 398 FOR" 

! NIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT " 

IIS EXP*-.001*HEIGHT ABOVE G5 METERS AND HEIGHT" 

‘TO THE -1.4G BELOW G5 METERS" 

IF Range* I) <65 THEN 

Ct eZ=75 . 5* Range* I ) A *-1 .46) 

ELSE 

CteZ=3.GGE-Z*EXP* -. 001 ♦ Range* I ) ) 

END IF 
CASE 4 

‘.CALCULATE THE TEMPERATURE STRUCTURE FACTOR BASED ON ASSUMPTION 
I THAT IT’S PROPORTIONAL TO HEIGHT TO THE -4/3 POWER 

! THAT THE SURFACE VERTICAL HEAT FLUX IS .095 C.n./sec 

‘THE DEPENDENCE WITH HEIGHT ABOVE THE INVERSION I PULLED FROM 
! MY LEFT EAR 

‘EQUATION FROM NEFF, 1975 
SELECT Range* I ) / I nver 
CASE <.9 

C3=* * .0Z4)*< T) A * .667) ) 

CteZ=C3* * Range* I ) ) A *-1 .33) 

CASE .9 TO 1 
C teZ«Ct eZ 
C 1 =C teZ 
R1 =Range* I ) 

CASE 1 TO 1 . 3 

Ct eZ= 1 0 A ( * LGTCC3* Range* I ) )-LGT(C3*R1 ) )+LGT( Cl ) ) 
R2=Range* I ) 

CZ=CteZ 
CASE ELSE 

CteZ = C3*( Range* I ) A ( - 1 . 33 ) - RZ A * - 1 . 33 ) ) + C2 
END SELECT 
END SELECT 

‘CALCULATE VELOCITY STRUCTURE FACTOR. FORMULA FROM GAYNOR 77 
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Epsilon=( . Z866/Range( I ) )*( 1+.07*( RangeC I ) ) A ( . B ) ) * ( 1 . 5 ) 
Cve2 e Z*Epsilon A ( .667) 

! CALCULATE ACOUSTIC REFRACTIVE INDEX STRUCTURE FACTOR. 

! FORMULA FROM TETARSKI 

CneZ( I )»( CteZ/( Z.98E+5) )+(CveZ/( C*C) ) 

! CALCULATE THE FRACTION OF POLIER BACKSCATTERED FROM INTERVAL. 

! EQUATION FROM NEFF , 1976 
Sigma=( .0039*(K A ( 1/3) )*CteZ )/( T) A Z 
! 

! CALCULATE THE EXCESS ATTENUATION. EQUATION FROM CLIFFORD, 1980 
! THE EXCESS ATTENUATION IS Exc_att 
IF I <50 THEN 

Remainder-0 

ELSE 

IF I < 1 50 THEN 

Renainder=Range( I ) MODULO 10 

ELSE 

Renai nder^RangeC I ) MODULO Z0 
END IF 
END IF 

IF Renainder-0 THEN 
Rho=0 
L*0 
H-0 

Rge=Range( I ) ! CONSTANT IN INTEGRAL 

R=0 

FOR J=0 TO 2*1 

F=CneZ ( INT( J/Z+1 ) ) 

F=F*( 1-R/(Rge) A ( 1 . 67 ) + ( R/ ( Rge ) ) A ( 1 .67) ) 

IF J>0 THEN 

IF J<Z*I THEN 

IF INT( J/Z )= J/Z THEN 
L=L+F 
F=0 

ELSE 

H=H+F 
F=0 
END IF 
END IF 
END IF 
Rho=Rho+F 
R=R+ 1 

NEXT J 

Rho=Rho+4* L+Z*H 

Rho=( < (Rho*.33)*K*K*1 .46> A (-.6> ) 

N=( Arrt_diam/Rho ) A Z 
IF N<= 1 THEN 

Ze= 1 /( 1+N) 

ELSE 

Ze= 1 . 5/( 1+N) 
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! STEP OF 1 . 5®=>SEE CLIFFORD 1980 
END IF 

Exc_att=Ze*Ze 

ELSE 

Exc_a-t-t-Exc_att 
END IF 
! 

! CALCULATE THE POWER BACKSCATTERED 

Pou_back( I ) =( Pou_trans*Ef f -Sumpou_back )*EXP( - Atten*Range( I ) ) 
Pou_back( I )=Pou_back( I )*I rrterval *Exc_atL* Sigma 
Sumpou_back=Sumpou_back+Pou_back( I ) 

! CALCULATE THE POWER BACKSCATTERED TO THE ANTENNA 
Pou_ret( I )=Pou_back( I ) *EXP< -At t en*Range< I ) ) * Ant_area*G*Ef f /Range( I ) * 2 
PRINT "RANGE®" , Range( I ) 

PRINT "POWER RETURNED®" ,Pou_ret( I ) 

PRINT "FREQUENCY®" ,Freq 
1 = 1 + 1 

UNTIL Pou_ret (1-1 ) >0 AND Pou_ret( 1-1 )<®Noise 
Ran( I nd , I f req ) =Range( I - Z ) 

IF Ran( I nd , If req ) > Ymax THEN 
Ymax=Ran( Ind , I f req ) 

END IF 

IF Ran< I nd , I f req ) <Ymin THEN 
Ymin=Ran( Ind,Ifreq) 

END IF 
I nd=I nd+ 1 
Again=1 
FOR J=0 TO I 

Pou_back ( J ) =0 
Pou_ret ( J )=0 
Range( J )®0 
CneZ ( J ) =0 
NEXT J 

NEXT Pou_trans 
I f r eq=I f req+ 1 
Ind® 1 
NEXT Freq 
! 

| 

OUTPUT KBD; "K" ; 

PRINT "INPUT CONDITIONS" 

! 

PRINT " " 

PRINT USING "K" ; "TEMPERATURE® 

PRINT USING "K" ; "ATMOSPHERIC PRESSURE® 

PRINT USING "K"; "WATER VAPOR PRESSURE® 

PRINT USING "K"; "PULSE LENGTH® 

PRINT USING "K"; "TRANSDUCER EFFICIENCY® 

PRINT USING "K M ; "POWER TRANSMITTED® 

PRINT "TEMPERATURE STRUCTURE PROFILE USED ", Prof lie 



" ,Temp,"CELCIUS" 

" , Atom_pres , " nb" 

" , Wat_pres , "mb" 

" ,Pul se , "ms" 

" ,Eff 

" ,Pou_trans," WATTS" 
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IF Profiled THEN 



PRINT "INVERSION HEIGHT 
END IF 

SELECT Profile 
CASE 1 

PRINT " 1 

PRINT " 

PRINT " 

CASE 2 

PRINT " 2 

PRINT " 

PRINT " 

PRINT M 
PRINT " 

CASE 3 

PRINT " 3 

PRINT " 

PRINT " 

PRINT " 

PRINT " 

CASE 4 

PRINT " 4 

PRINT " 

END SELECT 



" , I nver 



A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 
AS PRESENTED IN UALTERS/KUNDEL 198! PAGE 398 M 
WHICH GIVES A HEIGHT TO THE - 1 . 1 G PROFILE" 

A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 
AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 " 
WHICH GIVES A HEIGHT TO THE -1.16 PROFILE" 

BUT WITH A FACTOR OF TWO TO APPROXIMATE" 
LOOKING UP A THERMAL PLOOM" 

A TEMPERATURE STRUCTURE PROFILE BASED ON DATA" 
AS PRESENTED IN UALTERS/KUNDEL 1981 PAGE 398" 
FOR NIGHT CONDITIONS. THE DEPENDENCE WITH" 
HEIGHT IS EXP< - . 001 * HEIGHT ABOVE G5 METERS" 

AND HEIGHT TO THE -1.4G BELOW G5 METERS" 

A TEMPERATURE STRUCTURE PROFILE BASED ON A" 
HEIGHT TO THE -4/3 " 



j 

PRINT USING "K" ; "THE MINIMUM DETECTABLE SIGNAL WAS SET AT" , Noise , "WATTS" 
j 

PRINT " " 



! 



! 

PRINT "HIT CONTINUE TO CONTINUE" 
PAUSE 

OUTPUT KBD; "K" ; 



! 



GRAPHICS ON 
VIEWPORT 15,120, 10,70 



! PLOT RANGE VERSUS POWER TRANSMITTED FOR VARIOUS FREQUENCIES 

Xnin=0 

Xmax e =600 

Ymi n=PROUND( Ymin-50,2 ) 

Ymax=PROUND< Ymax+50, Z ) 

Dx= 100 
Dy= 1 00 

Xrange=Xmax-Xmin 
Yrange=Yroax-Ynin 
WINDOW Xmin , Xmax , Ymin , Ymax 
AXES Dx ,Dy ,Xmin, Ynin, 1 ,1 
CLIP OFF 
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! LABEL PLOT 
CSIZE 4, .6 
LDIR 0 
LORG 5 

MOVE . 5*Xrange , Ynax* 1 .05 

LABEL “RANGE VERSUS POWER TRANSMITTED" 

! LABEL HORIZONTAL AXES 

LDIR 0 

LORG 5 

FOR J-0 TO Xnax STEP Dx 
CSIZE 4, .6 

MOVE J # Ynin- . 05* Yrange 
LABEL J 
NEXT J 

MOVE . 5* Xrange , Ynin- . 1 * Yrange 
CSIZE 4, .6 

LABEL “POWER TRANSMITTED ( WATTS )“ 

! LABEL VERTICAL AXES 
LORG 8 

FOR J-Ymin TO Ynax STEP Dy 
CSIZE 4 , . G 

MOVE Xnin- . 0025* Xrange , J 
LABEL USING "K“;J 
NEXT J 
LDIR PI/2 
LORG B 

MOVE Xnin- . 1 5* Xrange , . 5* Yrange + Ymn 
CSIZE 4, .6 

LABEL “RANGE ( METERS )“ 

CLIP ON 
I nd= 1 
I f req* 1 

FOR Freq=500 TO 1500 STEP 500 
FOR Pou_trans=50 TO 500 STEP 50 
IF INT( Ifreq/2 ) = If req/ 2 THEN 
LINE TYPE 3 
ELSE 

LINE TYPE 1 
END IF 

PLOT Pou_irans , Ran( I nd , I f req ) 

I nd=I nd+ 1 
NEXT Pou_trans 
LDIR 0 
LORG 2 
CSIZE 3 , . 6 

MOVE Pou_trans-50, Ran( Ind- 1 , If req ) 
LINE TYPE 1 

LABEL USING "K“ ; Freq , “Hz “ 

I f req=I f req+ 1 
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I nd= 1 
NEXT Freq 

! 

PRINT "HIT CONTINUE TO CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT KBD; "K" ; 

| 

! 

LINPUT "WOULD YOU LIKE TO MAKE ANOTHER RUN (ENTER Y OR N)?",X$ 
Again=FNYes( X$) 

Mess_up= 1 
WHILE Mess_up= 1 
SELECT Again 
CASE 1 

LINPUT "WOULD YOU LIKE TO CHANGE A VARIABLE (ENTER Y OR N)?",X$ 
Neu_va=FNYes( X$) 

Ano_change= 1 
WHILE Ano_change=1 
SELECT Neu_va 
CASE 1 

PRINT " VARIABLE CURRENT VALUE" 

PRINT USING "K" ; " 1 TEMPERATURE " , Tenp , "CELSIUS" 

PRINT USING "K" ; " 2 ATOMOSPHERIC PRESSURE " , Aton_pres , " mb" 
PRINT USING "K " ; " 3 WATER VAPOR PRESSURE " , Wat_pres , " mb" 
PRINT " " 

PRINT USING "K " ; " 4 ANTENNA DIAMETER " , Ant_diam , " m." 

PRINT USING "K" ; "5 PULSE LENGTH ", Pulse," ms" 

PRINT USING "K " ; " B POWER TRANSMITTED " ,Pou_trans , " WATTS" 
PRINT USING "K" ; "7 ATMOSPHERIC PROFILE ".Profile 
PRINT " " 

PRINT "ENTER THE NUMBER OF THE VARIABLE YOU WISH TO" 

PRINT "CHANGE- 
INPUT Var 
SELECT Var 

CASE 1 

INPUT "TEMPERATURE=" ,Tenp 
CASE 1 

INPUT "ATMOSPHERIC PRESSURE IN mb=" , Aton_pr es 
CASE 3 

INPUT "WATER VAPOR PRESSURE IN mb=" , Wat_pres 
CASE 4 

INPUT "ANTENNA DIAMETER IN m. =" , Ant_diam 
CASE 5 

INPUT "PULSE LENGTH IN ns=" ,Pulse 
CASE G 

INPUT "POWER TRANSMITTED IN WATTS®" ,Pou_trans 
CASE 7 

PRINT " 1 ==> Ci A Z PROFILE OF Z A (-1.16)" 

PRINT " FROM WALTERS/KUNDEL 1981" 
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CASE 



PRINT M Z -«> SAME AS ONE BUT UITH FACTOR" 
PRINT " OF TUO TO APPROXIMATE LOOKING" 

PRINT " UP A THERMAL PLOOM" 

PRINT " 3 = = > Ct A Z PROFILE OF EXP< - . 00 1 * Z > " 

PRINT M FROM WALTERS/ KUNDEL 1981." 

PRINT " 4 — > CT A Z PROFILE OF Z A (-4/3>" 

INPUT "ENTER NUMBER OF DESIRED PROFILE" , Prof i 1 e 



IF Prof i 1 ©=4 THEN 

INPUT "HEIGHT OF INVERSION IN METERS'" , Inver 
END IF 
CASE ELSE 

PRINT Var,"IS NOT ONE OF THE OPTIONS" 

END SELECT 

LINPUT "MADE ANOTHER CHANGE( ENTER Y OR N)?",X$ 
Neu_va=FNYes( X$) 

Mess^up^Z 

Z 



Mess_up=Z 
Ano_change=Z 
CASE ELSE 

PRINT "YOUR CHOICES WHERE Y OR N MU!" 
Ano__change= 1 

LINPUT "CHANGE A VARIABLE (ENTER Y OR N)?",X$ 
Neu_va=FNYes( X$) 

END SELECT 
END WHILE 
CASE Z 



Mess__up=Z 
CASE ELSE 

PRINT "YOUR CHOICES WHERE Y OR N MUM" 
Mess__up= 1 

LINPUT "MAKE ANOTHER RUN (ENTER Y OR N)?",X$ 
Again=FNYes( X$) 

END SELECT 
END WHILE 
END WHILE 

PRINT "THAT’S ALL, FOLKS" 
end 



ICALCULTE THE ATTENUATION 

sub Att enuat ion( At om__pres , At ten , Freq , Temp , Wat _ pres ) 

! THIS SUBPROGRAM CALCULATES THE ATTENUATION OF SOUND 
! IN AIR BASED UPON EQUATIONS IN NEFF 1975 



INPUT ATMOSPHERIC PRESSURE IN MILLIBARS 
FREQUENCY OF SOUND WAVE IN HERTZ 
TEMPERATURE IN DEGREES CELCIUS 
WATER- VAPOR PRESSURE IN MILLIBARS 
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OUTPUT ATTENUATION IN 1/METERS 



VARIABLES 

Aton_pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 
Atten 

Att nax 



! F 

! Froax 

! Freq 
! H 

! Pstar 
! 

! T enp 
! Tstar 

! Uat_pres 

H= I 00*Uat_pres/ At on_pres 
Tstar = ( 1 . 8* Tenp + 49Z ) / 5 I 9 
Pstar=At on_pres/ 1014 

F max = ( 1 0+6600* H+44400*H*H)*Ps tar /Tstar* . 8 

At t_nax = . 0078* F max * Tst ar* ( -Z.5)*EXP(7.77*( 1-1 /Tstar) ) 

F=Freq/Fnax 

At t en= ( At t_max / 304 . 8 ) * ( ( . 1 8*F ) ( Z*F*F/( 1 + F*F ) ) * Z ) * . 5 

Atten=(Atten+l.74E-10*Freq*Freq)/4.35 

SUBEND 

i 

def FNYes(X$> 

! THIS FUNCTION INTERPRETS THE OPERATORS RESPONSES TO YES NO QUESTIONS 
! 

! INPUT X$ 

! 

! OUTPUT FNYes 

! VARIABLES 

! Tenp$ VARIABLE STRING USED IN FUNCTION YES. 

! X$ STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 

! RESPONSE TO YES OR NO QUESTION. 

DIM TempSI 1 ] 

T enp$[ 1 , 1 1 = TRI M$( X$) 

SELECT Tenp$ 

CASE #, Y M , M y M 
RETURN 1 
CASE 1" , M n M 
RETURN Z 
CASE M M 



ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 
SUBPROGRAM ATTENUATION. 

VARIABLE IN SUBPROGRAM ATTENUATION. IT IS THE 
ATTENUATION AT THE FREQUENCY OF THE MAXIMUM 
ATTENUATION FOR THE INPUT CONDITIONS. 

VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 

FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
ATTENUATION. 

INPUT FREQUENCY OF ECHOSOUNDER. 

VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 

VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 

INPUT TEMPERATURE IN DEGREES CELSIUS. 

INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 
IN SUBPROGRAM ATTENUATION. 

ATMOSPHERIC UATER PRESSURE IN MILLIBARS. INPUT BY 
OPERATOR. 



I 4 O 



RETURN 1 
CASE ELSE 

RETURN -2 
END SELECT 

FNEND 

! 
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APPENDIX F 
COMPUTER PROGRAM 6 



! FULLER, ROBERT PR0G_G NOISE 

! 

! 10 SEP 85 

j ******* PURPOSE* ****** 



! THIS PROGRAM WILL TAKE PARAMETERS OF AN ECHOSOUNDER AND 

! ESTIMATE THE RANGE AS A FUNCTION OF BACKGROUND NOISE 
! THE FOLLOWING INPUTS ARE REQUIRED: 

! ATMOSPHERIC DATA 

I 1 ) ATOMOSPHERIC PRESSURE in Millibars 

) 2 ) SELECT ONE OF FOUR ATMOSPHERIC PROFILES FOR THE TEMPERATURE 

! STRUCTURE PROFILE. 

! a)FOR ONE PROFILE THE HEIGHT OF THE INVERSION LAYER IS ALSO INPUT 

! 3 ) TEMPERATURE IN DEGREES CELSIUS 

! 4 > WATER VAPOR PRESSURE IN millibars 

J ECHOSOUNDER DATA 

! 5)FREQUENCY 

! B) ANTENNA DIAMETER IN METERS 

! 7 )PULSE LENGTH OF THE TRANSMITTED ACOUSTIC ENERGY 

! 8 )POWER TRANSMITTED 

! 

! THE PROGRAM OUTPUTS THE FOLLOWING GRAPH TO AN EXTERNAL PLOTTER: 

! 1)RANGE(m.) VERSUS BACKGROUND NOISE IN WATTS 



♦♦♦♦♦VARIABLES***** 



Again 

Ano_change 



Ant_area 

Ant_diam 
Atom __pres 
Att en 

C 

Cl 

C3 

CneZ< * > 



TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 
WISHES TO MAKE ANOUTHER RUN. 

TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 
WISHES TO CHANGE THE VALUE OF INPUT DATA BEFORE MAKING 
ANOUTHER RUN. 

ECHOSOUNDER ANTENNA AREA CALCULATED FROM INPUT OF 
ANTENNA DIAMETER. 

INPUT OF ECHOSOUNDER ANTENNA DIAMETER IN METERS. 

INPUT OF ATMOSPHERIC PRESSURE IN mb. 

ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 
SUBPROGRAM ATTENUATION. 

VELOCITY OF SOUND CALCULATION AT ZERO CELSIUS. 

VARIABLE USED TO SCALE INVERSION HEIGHT IN 
TEMPERATURE STRUCTURE PROFILE THREE. 

VARIABLE USED IN CALCULATING THE FOURTH CteZ PROFILE. 

ARRAY OF VALUES OF THE ACOUSTIC REFRACTIVE INDEX 
PARAMETER. CALCULATED BASED ON SELECTION OF 
PROFILE FOR TEMPERATURE STRUCTURE PARAMETER 
AND ASSUMED PROFILE FOR THE VELOCITY STRUCTURE 
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CteZ 


PARAMETER. 

VALUE OF THE TEMPERATURE STRUCTURE PARAMETER. 

THE OPERATOR SELECTS A PROFILE FROM SEVERAL PROVIDED. 


CveZ 


VALUE OF THE VELOCITY STRUCTURE PARAMETER. 

VALUES BASED ON CALCULATION USING ASSUMED DISSAPATION 
RATE. 


Dx 

Dy 

Epsilon 

Eff 

Eff 


STEP SIZE FOR X AXIS FOR VARIOUS PLOTS. 

STEP SIZE FOR Y AXIS FOR VARIOUS PLOTS. 
DISSAPATION RATE USED IN CALCULATION OF CveZ. 
TRANSMISSION EFFICIENCY OF ECHOSOUNDER. 

EFFICIENCY OF ECHOSOUNDER UHEN RECEIVING ACOUSTIC 
BACKSCATTER. 


Exc_att 

F 


EXCESS "ATTENUATION" AT GIVEN RANGE. 

VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 


Fnax 


FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
ATTENUATION. 


Freq 

G 

H 

I 


INPUT FREQUENCY OF ECHOSOUNDER. 

ANTENNA EFFECTIVE APERATURE FACTOR. 

VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 
MAIN INDEX FOR LOOP TO CHECK RETURNED POUER UNTIL 
POUER UAS LESS THAN MINIMUM DETECTABLE SIGNAL. 


IND 

Interval 

Inver 

J 

K 

L 

Mess_up 


INDEX USED FOR LOOP FOR DIFFERENT NOISE CUT OFF LEVELS 
MAXIMUM EFFECTIVE-SCATTERING VOLUME THICKNESS. 

HEIGHT OF INVERSION LAYER. 

FIRST ORDER INDEX FOR ASSORTED LOOPS. 

UAVENUMBER 

SECOND ORDER INDEX FOR ASSORTED LOOPS. 

TEST VALUE TO CHECK RESPONSES UHEN OPERATOR IS INPUTTING 
RESPONSES. 


N 


TEST VALUE IN CALCULATION OF COHERENCE LENTH IN CALCULATION 
OF EXCESS ATTENUATION. 


Neu_va 


TEST VALUE IN CHECKING IF OPERATOR UISHES TO CHANGE 
A VARIABLE BEFORE A NEU RUN. 


Noise 


ASSUMED MINIMUM DETECTABLE SIGNAL. 



Pou_back ( * ) POUER BACKSCATTERED FROM GIVEN RANGE. 

Pou_ret ( * ) POUER BACKSCATTERED TO ECHOSOUNDER FROM GIVEN RANGE. 
Pou trans INPUT OF POUER SUPPLIED TO TRANSDUCER OF ECHOSOUNDER 



Prof ile 
Pstar 


OPERATOR INPUT OF CteZ PROFILE FROM AVAILABLE PROFILES. 
VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 


Pul se 


TRANSMITTED PULSE LENGTH IN MILLISECONDS. INPUT BY 
OPERATOR. 


R 

R1 

Range ( * ) 
Remainder 


RANGE FOR INPUT TO INTEGRAL OF EXCESS ATTENUATION. 

VARIABLE USED IN SCALING INVERSION HEIGHT FOR PROFILE THREE 
ARRAY OF RANGE VALUES. 

REMAINDER OF MODULO FUNCTION USED TO DECREES THE 
NUMBER OF PASSES THROUGH THE INTEGRATION FOR EXCESS 
ATTENUATION. 


Rge 


VALUE OF RANGE IN INTEGRAL FOR EXCESS ATTENUATION. 
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! Rho CORRELATION LENGTH USED IN CALCULATION OF EXCESS 

I ATTENUATION. 

! Sigma FRACTION OF ENERGY BACKSCATTERED FROM GIVEN RANGE. 

I Speed_sound SPEED OF SOUND IN AIR AT INPUT TEMPERATURE. 

! T INPUT TEMPERATURE IN. DEGREES KELVIN. 

{ Temp INPUT TEMPERATURE IN DEGREES CELSIUS. 

! Sumpou_back 

! Temp$ VARIABLE STRING USED IN FUNCTION YES. 



I Tstar 

I Var 

! Wat_pres 
! 

! X 
! X$ 

! Xmax 
I Xmin 
! Xrange 
! Ylabel$ 

! Ymax 
! Y m i n 
! Yrange 
! Ze 



INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 
IN SUBPROGRAM ATTENUATION. 

USED TO SELECT THE VARIABLE THE OPERATOR WISHED TO 
CHANGE BEFORE MAKING ANOTHER RUN. 

ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 
OPERATOR. 

THIRD ORDER INDEX USED IN VARIOUS LOOPS. 

STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 
RESPONSE TO YES OR NO QUESTION. 

VALUE OF LARGEST VALUE OF X FOR VARIOUS PLOTS. 

VALUE OF SMALLEST VALUE OF X FOR VARIOUS PLOTS. 

VALUE OF RANGE OF X VALUES FOR VARIOUS PLOTS. 

LABEL ON Y AXIS PASSED TO SUBPROGRAM Semi_log FOR PLOTTING. 
LARGEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

SMALLEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

RANGE OF Y VARIABLES FOR VARIOUS PLOTS. 

INTERMEDIATE VALUE IN CALCULATION OF EXCESS ATTENUATION. 



! DECLARE VARIABLES 
INTEGER I 

DIM Pou_back( 5000) ,Pou_ret( 5000) ,Range( 5000) 

DIM CneZ( 5000) , Ran( 50) 

! 

PLOTTER IS 705 , "HPGL M 
LINE TYPE 1 
! 

! INPUT ATMOSPHERIC DATA 

INPUT “ENTER TEMPERATURE IN DEGREES CELCI US" , Temp 
INPUT “ENTER ATMOSPHERIC PRESSURE IN MILLIBARS" , Atom_pres 
INPUT “ENTER WATER VAPOR PRESSURE IN MILLIBARS" ,Wat_pres 



! INPUT ECHOSOUNDER DATA 

INPUT “ENTER ECHOSOUNDER FREQUENCY IN HERTZ", Freq 
INPUT “ENTER ECHOSOUNDER PULSE LENGTH IN MILLISECONDS^" , Pulse 
INPUT “ENTER THE ECHOSOUNDER ANTENNA DIAMETER IN METERS" , Ant_diam 
INPUT "ENTER THE POWER TRANSMITTED IN WATTS" ,Pou_trans 
G= . 40 ! ANTENNA EFFECTIVE APERATURE FACTOR 

Ef f = . 25 
! 
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! 

! SELECT THE PROFILE FOR THE TEMPERRTURE STRUCTURE 
PRINT “YOU HRVE R CHOICE OF TEMPERRTURE STRUCTURE PROFILES 11 
PRINT "PROFILE TO BE USED." 

PRINT "YOUR SELECTIONS R RE" 

PRINT " 1 

PRINT " 

PRINT " 

PRINT " " 

PRINT " Z 

PRINT " 

.PRINT " 

PRINT " " 

PRINT " 3 

PRINT " 

PRINT " 

PRINT " 

PRINT " 

PRINT " " 

PRINT " 4 

PRINT " 

I 

PRINT " " 

PRINT " " 

Mess_up= 1 
WHILE Mess_up=1 

INPUT "ENTER THE DESIRED PROFILE ( 1 OR Z OR 3 OR 4)" .Profile 
IF Prof i 1 e=4 THEN 

INPUT "ENTER THE HEIGHT OF THE INVERSION IN METERS" , I nver 
END IF 

IF Profile-! THEN 
Mess_up=0 

ELSE 

IF Prof il e=Z THEN 
Mess_up=0 

ELSE 

IF Profile-3 THEN 
Mess_up=0 

ELSE 

IF Profile-4 THEN 
Mess__up=0 

ELSE 

PRINT Profile," WRS NOT ONE OF THE OPTIONS! 
Mess_up= 1 
END IF 
END IF 
END IF 
END IF 
END WHILE 
OUTPUT KBD; "K" ; 



R TEMPERRRURE STRUCTURE PROFILE BRSED ON DRTR RS 
PRESENTED IN WRLTERS/KUNDEL 198! PRGE 398 WHICH 
GIVES R HEIGHT TO THE -1.16 PROFILE" 

THE TEMPERRTURE STRUCTURE PROFILE RBOVE BUT" 

WITH R FRCTOR OF TWO TO RPPROXIMRTE LOOKING" 

UP R THERMRL PLOOM" 

R TEMPERRTURE STRUCTURE PROFILE BRSED ON DRTR RS 
PRESENTED IN WRLTERS/KUNDEL 1981 PRGE 398 FOR" 
NIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT " 
IS EXP(-.001*HEIGHT RBOVE 65 METERS RND HEIGHT" 
TO THE -1.46 BELOW G5 METERS" 

R TEMPERRTURE STRUCTURE PROFILE BRSED ON fl" 
HEIGHT TO THE -4/3 " 



145 



! 

! 

Again= I 
WHILE Again= I 
FOR J=1 TO I 

Pou_back ( J ) =0 
Pou_ret ( J )=0 
Range< J )=0 
CneZ< J>«0 
NEXT J 

j 

! CONVERT TEMPERATURE TO KELVIN 
T=Temp+Z73 
! 

! CALCULATE SPEED OF SOUND BASED ON INPUT TEMPERATURE< CELCIUS ) 
Speed_sound=Z0. 05* ( T) A .S 

! CALCULATE THE SPEED OF SOUND AT 0 DEGREES CELCIUS 
C=Z0. 05*Z73 A . 5 
! 

1=1 ! INDEX FOR INNER LOOP (POWER RETURNED TO NOISE) 

I nd= 1 ! INDEX FOR SECOND LOOP (NOISE) 

Pou_back(0)=0 UNITIALIZE VARIABLE FOR POWER BACKSCATERRED 
I nterval =( Speed_sound*Pul se* 1 . E-3 ) /Z 
Ant__area=PI *( Arrt__dian/Z ) A Z 
( 

! CALCULATE RETURNED SIGNAL POWER UNTIL IT IS LESS THAN BACKGROUND 
IFOR EACH FREQUENCY. 

Range( 0 ) =0 
Noise= 1 .E-18 
FOR Ind= 1 TO 5 

CALL Attenuat ion( Atom_pres , Atten,Freq,Temp, Wat_pres ) 

K=Z*PI *Freq/Speed_sound ! WAVENUMBER 
1 = 1 

REPEAT 

Range( I )=Range( I- 1 ) + Z 
Sunpou_back =0 
SELECT Profile 
CASE 1 

ITHIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
IAS PRESENTED IN WALTERS/KUNDEL 1381 PAGE 398. 

! WAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
! TULAROSA BASIN, NEW MEXICO. 

! AN ADDITONAL FACTOR OF 1 . 8 WAS INCLUDED AS PER DIFFERENT 
! AVERAGING TIME. 

CteZ=Z . 1 Z*Range( I ) A <-1 . 1G) 

CASE Z 

ITHIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
IAS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398. 

(WAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
l TULAROSA BASIN, NEW MEXICO. 
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! AN ADDI TONAL FACTOR OF t . 8 WAS INCLUDED AS PER DIFFERENT 
I AVERAGING TIME AND AN ADDITIONAL FACTOR OF TWO UAS 
! INCLUDED TO APPROXIMATE LOOKING UP A THERMAL PLOOM 
C-teZ“Z*Z . 1 Z*Range( I ) A ( - 1 . 1 G) 

CASE 3 

ITHIS TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS" 

PRESENTED IN UALTERS/KUNDEL 1981 PAGE 398 FOR" 

! NIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT " 

IIS EXP( - . 001 * HEIGHT ABOVE G5 METERS AND HEIGHT" 

I TO THE -1.4G BELOU G5 METERS* 1 
IF Range(IXG5 THEN 

CteZ=75. 5*Range( I > * ( - 1 . 4G ) 

ELSE 

CteZ=3 . 66E-Z*EXP( - . 001 * Range ( I > > 

END IF 
CASE 4 

! CALCULATE THE TEMPERATURE STRUCTURE FACTOR BASED ON ASSUMPTION 
! THAT IT’S PROPORTIONAL TO HEIGHT TO THE -4/3 POWER 

! THAT THE SURFACE VERTICAL HEAT FLUX IS .095 C.n./sec 

! THE DEPENDENCE WITH HEIGHT ABOVE THE INVERSION I PULLED FROM 
* MY LEFT EAR 

! EQUATION FROM NEFF , 1975 
SELECT Range( I > / I nver 
CASE <.9 

C3=( ( . 0Z4 )*( T) A ( . G67 ) > 

CteZ=C3* ( Range ( I ))"(-! .33) 

CASE .9 TO 1 
Ct eZ = Ct eZ 
C 1 =CteZ 
R1 =Range( I ) 

CASE 1 TO 1.3 

Ct eZ= 1 0 A ( ( LGT ( C3*Range( I ) >-LGT(C3*Rl ) ) + LGT(Ct ) ) 
RZ=Range( I > 

CZ=CteZ 
CASE ELSE 

CteZ=C3*( Range( I > A ( - 1 . 33 >-RZ A ( - 1 . 33 > ) + CZ 
END SELECT 
END SELECT 

! CALCULATE VELOCITY STRUCTURE FACTOR. FORMULA FROM GAYNOR 77 
Eps il on=( . Z866/Range( I ) >*( 1 + . 07* ( Range( I)) A (.6)) A (1.5) 
CveZ=Z*EpsilorT( .GG7> 

! CALCULATE ACOUSTIC REFRACTIVE INDEX STRUCTURE FACTOR. 

! FORMULA FROM TETARSKI 

CneZ ( I > = (CteZ/(Z.98E+5> > + (CveZ/( C*C ) ) 

! CALCULATE THE FRACTION OF POWER BACKSCATTERED FROM INTERVAL, 
i EQUATION FROM NEFF, 1975 
S igna=( .0039*(K A < 1 /3 > >*CteZ )/( T) A Z 
| 

! CALCULATE THE EXCESS ATTENUATION. EQUATION FROM CLIFFORD, 1980 
ITHE EXCESS ATTENUATION IS Exc_att 
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IF I <5® THEN 

Remainder®® 

ELSE 

IF I<15® THEN 

Remainder®Range( I ) MODULO 1® 

ELSE 

Remainder=Range( I ) MODULO 2® 

END IF 
END IF 

IF Remainder®® THEN 
Rho=® 

L®® 

H®® 

Rge=Range( I ) ! CONSTANT IN INTEGRAL 

R®® 

FOR J®® TO 2*1 

F=Cne2< INT( J/2+1 ) ) 

F=F*( 1-R/(Rge) A ( 1 . 67 ) +< R/ < Rge ) ) A < 1 .67) ) 

IF J>® THEN 

IF J<2 * I THEN 

IF INT( J/2 ) = J/2 THEN 
L=L + F 
F=® 

ELSE 

H=H4F 

F=® 

END IF 
END IF 
END IF 
Rho=Rho+F 
R®R+ 1 

NEXT J 

Rho=Rho+4*L+2*H 

Rho®(<(Rho*.33>*K*K*1 .4G) A <-.G) > 

N=( Ant_diam/Rho) A 2 
IF N<= 1 THEN 

Ze= 1 / ( 1+N) 

ELSE 

2e= 1 . 5/ ( HN > 

ISTEP OF 1 .5==>SEE CLIFFORD 198® 

END IF 

Exc_att=Ze*Ze 

ELSE 

Exc_at t=Exc_att 
END IF “ 

l 

! CALCULATE THE POWER BACKSCATTERED 

Pou__back ( I ) = ( Pou__trans*Ef f-Sumpou__back ) *EXP( - At t en* Range ( I ) ) 
P ou_back ( I )=P outback ( I >*Interval*Exc_att*Sigma 
Sumpou_back = Sumpou_back +Pou_back ( I ) 
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! CALCULATE THE POWER BACKSCATTERED TO THE ANTENNA 
Pou_ret( I )=Pou_back( I )*EXP( -At ten* Range ( I ) )* Ant _ area* G*Ef f /Range ( I ) *2 
PRINT "RANGE*” ,Range( I ) 

PRINT "POWER RETURNED*" ,Pou_ret( I ) 

PRINT "NOISE*” ,Noise 
1*1 + 1 

UNTIL Pou_ret< I- 1 )>0 AND Pou_ret ( I- 1 ) <=Noise 
Ran( Ind)-Range( I-Z ) 

Again* 1 
FOR J*0 TO I 

P ou_back ( J )=0 
Pou_ret ( J ) =0 
Range( J )*0 
CneZ(J)=0 
NEXT J 

Noise*Noise* 1 0 
NEXT Ind 
! 

! 

OUTPUT KBD; "K" ; 

PRINT "INPUT CONDITIONS" 

! 

PRINT " " 

PRINT USING ”K" ; "TEMPERATURE* " , Temp f "CELCIUS” 

PRINT USING "K"; "ATMOSPHERIC PRESSURE- " , Atom_pres , " mb” 

PRINT USING "K" ; "WATER VAPOR PRESSURE* " , Wat_pres , " mb" 

PRINT USING "K" ; "FREQUENCY* ",Freq," Hz" 

PRINT USING "K" ? "PULSE LENGTH* " ,Pul se , " ns" 

PRINT USING "K" ; "TRANSDUCER EFFICIENCY* ",Eff 

PRINT USING "K"; "POWER TRANSMITTED* ".Poujrans," WATTS" 

PRINT "TEMPERATURE STRUCTURE PROFILE USED ", Prof lie 
IF Prof il e=4 THEN 

PRINT "INVERSION HEIGHT ", Inver 
END IF 

SELECT Profile 
CASE 1 

PRINT " 1 

PRINT " 

PRINT " 

CASE 2 

PRINT " 2 

PRINT " 

PRINT " 

PRINT " 

PRINT " 

CASE 3 

PRINT " 3 

PRINT " 

PRINT " 

PRINT " 



A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 
AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 " 
WHICH GIVES A HEIGHT TO THE -1.16 PROFILE" 

A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 
AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 " 
WHICH GIVES A HEIGHT TO THE -1.16 PROFILE" 

BUT WITH A FACTOR OF TWO TO APPROXIMATE" 
LOOKING UP A THERMAL PLOOM" 

A TEMPERATURE STRUCTURE PROFILE BASED ON DATA" 
AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398" 
FOR NIGHT CONDITIONS. THE DEPENDENCE WITH" 
HEIGHT IS EXP( - . 001 * HEIGHT ABOVE 65 METERS" 
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PRINT " FIND HEIGHT TO THE - 1 . 4G BELOU G5 METERS’* 

CASE 4 

PRINT 4 A TEMPERFITURE STRUCTURE PROFILE BFISED ON FI” 

PRINT ** HEIGHT TO THE -4/3 M 

END SELECT 
! 

! 

PRINT " M 
! 

! 

PRINT "HIT CONTINUE TO CONTINUE” 

PFIUSE 

OUTPUT KBD;"K"; 

! 

GRAPHICS ON 
VIEWPORT 15, 120,10,70 
(PLOT RANGE VERSUS NOISE 
! 

Xnin=- 1 8 
Xnax=- 1 3 
Ynax=4 
Yni n=Z 
Dx = 1 
Dy= 1 

Yrange=2 

Xrange=5 

WINDOW Xni n , Xmax , Ynin , Ynax 
AXES Dx , Dy , Xnin , Ynin , 1 , I 
CLIP OFF 
! 

! LABEL PLOT 
CSI2E 4 , . G 
LDIR 0 
LORG 4 

MOVE Xnin+ . 5* Xrange , Ynin+ 1 * Yrange 
LABEL "RANGE vs NOISE" 

! LABEL HORIZONTAL AXES 
LDIR 0 
LORG 5 

FOR J=Xnin TO Xnax STEP Dx 
CSIZE 4, .6 

MOVE J- . 01 3* Xrange , Ynin- . 0G7* Yrange 
LABEL "10" 

MOVE J+ . 01 3* Xrange , Ynin- . 033* Yrange 
CSIZE 2 
LABEL J 
NEXT J 

MOVE Xnin+ . 5*Xrange , Ynin- . 1Z*Yrange 
LORG 5 
CSIZE 4 , . G 
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LABEL "NOISE ( WATTS ) " 

! LABEL VERTICAL AXES 
LORG 8 

FOR J=Ynin TO Ynax STEP Dy 
CSIZE 4 , . G 

MOVE Xmi n- .01 *Xrange , J 
LABEL 10 
CSIZE Z 

MOVE Xroin- . 00Z5*Xrange , J+ . 03*Yrange 
LABEL J 
NEXT J 
LDIR PI/Z 
LORG G 

MOVE Xnin-. 1 *Xrange , Ynin+ . 5* Yrange 
CSIZE 4 , . 6 

LABEL "RANGE (METERS)" 

CLIP ON 
Ind= 1 

Noise= 1 . E- 1 8 
FOR I nd= 1 TO 5 

PLOT LGT( Noise > ,LGT(Ran< Ind) > 
Noise=Noise* 1 0 
NEXT Ind 
! 

PRINT "HIT CONTINUE TO CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT KBD; "K" ; 

! 

! 



LINPUT "WOULD YOU LIKE TO MAKE ANOTHER RUN (ENTER Y OR N>?" ,X$ 
Again=FNYes(X$> 

Mess_up= 1 
WHILE Mess_up*1 
SELECT Again 
CASE 1 

LINPUT "WOULD YOU LIKE TO CHANGE A VARIABLE (ENTER Y OR N)?",X$ 
Neu_va-FNYes( X$) 



Ano_change= 1 



WHILE Ano_change=1 
SELECT Neu va 



CASE 1 




PRINT 


•• 


PRINT 


USING 


PRINT 


USING 


PRINT 


USING 


PRINT 


II M 


PRINT 


USING 


PRINT 


USING 


PRINT 


USING 



VARIABLE CURRENT VALUE" 

"K" ; " 1 TEMPERATURE ", Temp , "CELSIUS" 

"K " ; "Z ATOMOSPHERIC PRESSURE " , At on_pres , " mb' 
"K" ; " 3 WATER VAPOR PRESSURE " , Wat_pres , " mb" 

"K" ; "4 ANTENNA DIAMETER " , Ant_diam , " m." 
"K " ; "5 PULSE LENGTH ".Pulse," ns" 

" K" ; "G POWER TRANSMITTED " , Poujrans f " WATTS' 



PRINT USING "K" ; "7 ATMOSPHERIC PROFILE ", Prof lie 
PRINT " M 

PRINT "ENTER THE NUMBER OF THE VARIABLE YOU WISH TO" 

PRINT "CHANGE" 

INPUT Var 
SELECT Var 

CASE 1 

INPUT "TEMPERATURE 52 " .Temp 
CASE 2 

INPUT "ATMOSPHERIC PRESSURE IN mb=" , Atom_pres 
CASE 3 

INPUT "WATER VAPOR PRESSURE IN mb=" , Wat_pres 
CASE 4 

INPUT "ANTENNA DIAMETER IN n. =" , Ant_dian 
CASE 5 

INPUT "PULSE LENGTH IN ms=", Pulse 
CASE G 

INPUT "POWER TRANSMITTED IN WATTS=" ,Pou_i;rans 
CASE 7 

PRINT " 1 ==> Ct A Z PROFILE OF Z A (-1.16)" 

PRINT " FROM WALTERS/KUNDEL 1981" 

PRINT " 2 ==> SAME AS ONE BUT WITH FACTOR- 

PRINT " OF TWO TO APPROXIMATE LOOKING- 

PRINT " UP A THERMAL PLOOM" 

PRINT " 3 = =*> Ct A Z PROFILE OF EXP( - . 001 * 1 ) " 

PRINT " FROM WALTERS/KUNDEL 1981." 

PRINT " 4 = = > CT A 2 PROFILE OF Z A <-4/3>" 

INPUT "ENTER NUMBER OF DESIRED PROFILE" .Prof i 1 e 
IF Prof i 1 e=4 THEN 

INPUT "HEIGHT OF INVERSION IN METERS=" , I nver 
END IF 
CASE ELSE 

PRINT Var, "IS NOT ONE OF THE OPTIONS" 

END SELECT 

LINPUT "MADE ANOTHER CHANGE( ENTER Y OR N)?",X$ 
Neu_va=FNYes( X$) 

Mess_up-2 
CASE 2 



Mess_up=2 
Ano_change=2 
CASE ELSE 

PRINT "YOUR CHOICES WHERE Y OR N !!!!!" 
Ano_change= 1 

LINPUT "CHANGE A VARIABLE (ENTER Y OR N)?",XS 
Neu_va=FNYes( X$) 

END SELECT 
END WHILE 
CASE 2 



Mess_up=Z 
CASE ELSE 



I 52 



PRINT "YOUR CHOICES UHERE Y OR N ! I ! ! i I ,# 
Mess_up= 1 

LINPUT "MAKE ANOTHER RUN (ENTER Y OR N)7",X$ 
Again=FNYes( X$) 

END SELECT 
END WHILE 
END WHILE 

PRINT "THAT’S ALL, FOLKS" 

end 

i 



ICALCULTE THE ATTENUATION 

sub Att enuat ion( At om_pres , At t en , Freq , T emp , Wat_pres ) 

! THI S SUBPROGRAM CALCULATES THE ATTENUATION OF SOUND 
I IN AIR BASED UPON EQUATIONS IN NEFF 1975 
j 

! INPUT ATMOSPHERIC PRESSURE IN MILLIBARS 
I FREQUENCY OF SOUND WAVE IN HERTZ 

! TEMPERATURE IN DEGREES CELCIUS 

! WATER- VAPOR PRESSURE IN MILLIBARS 

! 

! OUTPUT ATTENUATION IN I /METERS 



VARIABLES 
Atom _pres 
At ten 

At t_nax 



F 

Fmax 

Freq 

H 

Pstar 

Temp 

Tstar 

Wat_pres 



INPUT OF ATMOSPHERIC PRESSURE IN mb. 

ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 
SUBPROGRAM ATTENUATION. 

VARIABLE IN SUBPROGRAM ATTENUATION. IT IS THE 
ATTENUATION AT THE FREQUENCY OF THE MAXIMUM 
ATTENUATION FOR THE INPUT CONDITIONS. 

VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 

FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
ATTENUATION. 

INPUT FREQUENCY OF ECHOSOUNDER. 

VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 

VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 

INPUT TEMPERATURE IN DEGREES CELSIUS. 

INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 
IN SUBPROGRAM ATTENUATION. 

ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 
OPERATOR. 



H= 1 00* Wat.__pres/Aton_pres 
Tstar = < 1 . 8« T enpH 492 ) / 5 1 9 
Pst ar = At on^pres/ 1014 

Fmax=( 1 0+6600* H + 44400* H* H )*Pstar /Tstar* . 8 

Att_max=. 0078* Fmax * Tst ar A < -Z .5)*EXP( 7.77* ( 1- 1 /Tstar ) ) 

F=Freq/Fmax 

At t en=( Att max/ 304. 8)*( ( . 1 8*F ) *Z + ( 2*F*F/< 1 +F*F ) )* 2 ) A . 5 



Atten=< Atten-U . 74E- 1 0*Freq*Freq ) /4 . 35 
SUBEND 
| 

def FNYes(X$) 

! THIS FUNCTION INTERPRETS THE OPERATORS RESPONSES TO YES NO QUESTIONS 
! 

! INPUT X$ 

! 

! OUTPUT FNYes 

! VARIABLES 

! Tenp$ VARIABLE STRING USED IN FUNCTION YES. 

! X$ STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 

! RESPONSE TO YES OR NO QUESTION. 

DIM Tenp$C 1 1 

Tenp$C 1 , 1 3 =TRIM$( X$ ) 

SELECT Temp$ 

CASE " Y" , "y" 

RETURN 1 
CASE "N M , "n" 

RETURN 2 
CASE M 11 

RETURN 1 
CASE ELSE 

RETURN -2 
END SELECT 



FNEND 
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